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I.  Introduction  and  Suomaurv 


High  angle -of -attack  phenomena  have  been  of  interest  to  aero- 
dynamicists,  aircraft  designers,  pilots  and  control  system  analysts 
ever  since  the  advent  of  modern  high  performance  aircraft.  Due  to 
the  concentration  of  inertia  along  the  fuselage,  the  modern  jet 
fighters  are  highly  susceptible  to  post-stall  departures  and  spin. 
Extensive  wind-tunnel  testing  and  radio-controlled  flight  testing  has 
been  done  over  the  last  twenty  years  to  gain  better  understanding  of 
the  dynamic  Instabilities  at  high  angles-of-attack.  A basic  problem 
has  existed  in  interpreting  this  data  and  in  making  predictions  of 
aircraft  dynamic  behavior  so  as  to  achieve  close  agreement  with  flight 
test  data. 

Aircraft  dynamic  behavior  at  high  angles-of-attack  is  highly 
nonlinear  and,  in  the  past,  there  has  been  a lack  of  suited>le  techniques 
for  analyzing  the  global  behavior  of  nonlinear  systems.  Under  am  on-going 
project  with  the  Office  of  Naval  Research,  Scientific  Systems,  Inc.  has 
developed  a new  approach  based  on  Bifurcation  ^alysis  and  Catastrophe 
Theory  Methodology  (BACTM) . The  approach  has  been  applied  to  specific 
jump,  hysteresis  and  limit  cycle  phenomena  such  as  roll -coupling,  pitch- 
up,  wing  rock,  buffeting,  departure  ^md  divergence.  Three  different  air- 
craft have  been  considered  for  con^rison  purposes,  and  it  has  been  shown 
how  different  types  of  instabilities  and  families  of  limit  cycles  arise 


aircraft  equlllbrluB  and  bifurcation  surfaces  is  given  in  zm  eight 
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dinwnsional  space  consisting  of  roll  rate,  pitch  rate,  yaw  rate,  zmgle- 
of>attack,  sideslip  angle,  elevator,  aileron  amd  rudder  deflections. 

Two  dimensional  projections  of  the  equilibrium  emd  bifurcation  surfaces 
provide  pictorial  representations  of  the  aircraft  global  stability  zmd 
control  behavior  at  high  2mgles-of-attack.  The  use  of  BACTM  for  spin 
entry,  spin  prevention,  stability  augmentation  at  high  angles-of -attack 
and  nonlinear  system  identification  are  also  considered. 

This  chapter  presents  a discussion  of  the  high  angle -of -attack 
stability  problem  followed  by  a brief  survey  of  the  previous  work  in 
this  area.  The  new  methodology  called  BkCTM  is  Introduced  in  section 
1.3  and  the  significant  results  obtained  during  the  current  stxidy  are 
suomuurized  in  section  1.4.  The  last  section  in  this  chapter  discusses 
the  overall  organization  of  the  report. 

(1.1)  Problem  Description. 

Aircraft  motion  at  high  angles-of'-attack  (a)  is  chuacterlzed  by 
l2u:ge  deviations  in  both  state  and  control  variables.  Linearized  emalysis 
of  aircraft  stability  and  control  provides  only  limited  localized  informa- 
tion regarding  aircraft  behavior  at  high  a.  There  is  a oon^lcuous  lack 
of  techniques  for  analyzing  global  stability  and  large  maneuver 
response  of  aircraft.  Only  specialized  nonllne^u:  phenomena  such  as 
roll-coupling  and  spin  have  been  partially  analyzed  in  an  Isolated  fashion 
There  is  a clear  need  for  a unified  approach  to  analyze  global  aircraft 
behavior  at  high  a In  an  Integrated  and  systematic  meuuier. 


IJ 

li 


A methodology  for  global  stability  emalysis  and  control  is  also 
a first  step  towards  safer  piloting  procedures,  stability  augmentation, 
control  system  design,  aircraft  model  structure  determination  and  design. 

A suitable  methodology  must  contribute  towards  solutions  to  these  problems 
in  addition  to  simply  giving  a yes/no  type  of  information  regarding 
stability.  In  other  words,  the  global  stability  analysis  should  also 
provide  guidelines  for  stabilizing  an  unstable  nonlinear  system.  Since 
aircraft  behavior  at  high  a is  known  to  exhibit  discontinuous  and  limit 
cycle  phenomena,  the  methodology  should  also  be  able  to  predict  and  explain 
such  phenomena. 

(1.2)  Previous  Wbrk 

Early  work  on  high  a stability  and  control  during  the  forties  and 
fifties,  concentrated  mainly  on  qxialitative  description  and  understanding 
Qf  the  phenomena  [1,2, 3, 4, 5].  The  study  of  specific  phencxnena  was 
motivated  by  observed  flight  behavior.  Attempts  were  made  to  understand 
these  pheomena  and  to  correct  them  via  changes  in  aircraft  design. 
Conventional  analysis  techniques  based  on  steady  state  assumptions  and 
linearization  were  used  to  obtain  quantitative  information.  Itowever, 
a lack  of  suitable  experimental  facilities,  analysis  techniques  and 
meager  c<»nputing  capabilities  hindered  progress  in  this  area.  The  major 
emphasis  during  the  sixties  was  placed  on  the  establishment  of  experimental 

facilities  such  as  radio *<;on trolled  flights  and  winc^^tunnels  for  spin 

If 

testing.  At  the  same  time,  simulation  techniques  were  used  to  reproduce 
with  six  degree-of-freedom  models  behavior  typical  of  high  a flights 
[6,7,8].  Again,  due  to  lack  of  analytical  techniques,  only  limited  pro- 
gress has  been  made  using  simulation  technlAies. 
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Recent  work  in  high  a areas  covers  a wide  spectrtsn  of  subjects 
reuiging  frctn  basic  reseeurch  in  aerodynamics,  wind  tunnel  and  flight 
testing  [9]  to  advemced  mathematical  techniques  for  stability  analysis 
[10]  euid  peurameter  identification  [11].  The  highly  nonlinear  zmd  complex 
nature  of  aircraft  dynamic  behavior  at  high  a coupled  with  poorly  under- 
stood aerodynamics  has  presented  formidable  difficulties  to  significant 
progress  in  this  area. 

The  motivation  for  the  reseiurch  effort  reported  here  came  from  an 
involvement  of  the  principal  investigator  with  a study  of  the  stability 
of  chemical  reactions  using  Catastrophe  Theory  [12,13]  . It  was  con- 
jectured that  the  discontinuous  and  limit  cycle  phenomena  at  high  angles- 
of-attack  could  be  analyzed  using  results  from  Topology,  Biftircatlon 
Analysis  and  Catastrophe  Theory  [14,15],  Recent  analysis  of  the  roll- 
coupling problem  by  Schy  and  Hannah  [16]  provided  further  support  to 
this  conjecture.  The  results  presented  here  have  not  only  confirmed 
our  original  conjecture,  but  have  also  provided  additional  insightc 
and  revealed  new  dynamic  phenomena  that  can  occur  at  high  a> 

(1.3)  New  Methodology 

The  application  of  recent  abstract  mathematical  results  from 
Topology,  Bifurcation  Analysis  and  Catastrophe  Theory  to  the  aircraft 
high  a stability  problem  has  resulted  in  the  development  of  a new  method- 
ology called  BACTM  (Bifurcation  Analysis  and  Catastrophe  Theory  Methodology) . 
BACTM  is  applicable  to  other  nonlinear  systems  such  as  submeurines,  missiles, 
aM  power  systems.  The  four  main  theorems  on  which  BACTM  is  based  are 
(i)  the  Center  Manifold  Theoremi  (ii)  Classification  Theorem  of  Elementary 
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Catastrophes;  (ill)  Hopf  Bifurcation  Theorem;  and  (iv)  Global 
Implicit  Function  Theorem.  These  theorems  are  discussed  in  Chapter  II 
emd  Appendix  A and  their  application  to  the  aircraft  stability  problem 
is  discussed  in  Chapter  III.  The  application  of  BACTM  provides  global 
information  regeurding  aircraft  dynamic  behavior  vdth  respect  to  changes 
in  initial  state  parameters,  and  control  variables.  The  use  of  BACTM 
regviires  adv2mced  ccxnputational  techniques  which  are  extremely  i]i^>ortant 
for  obtaining  quantlative  information  needed  for  stability  augmentation, 
control,  identification  2uid  design. 

(1.4)  Sumnary  of  Signific^tnt  Results 

Since  detailed  quantitative  results  of  the  study  are  presented 
in  Chapters  II  and  III,  we  briefly  summarize  here,  in  qualitative  terms, 
what  are  believed  to  be  some  of  the  significant  results.  The  current 
study  is  the  first  phase  of  a continuing  resezurch  effort  and,  therefore, 
represents  only  a first  step  in  the  direction  of  methodology  development 
emd  potential  applicability  of  BACTM.  The  contributions  of  the  present 
study  may  be  divided  into  three  general  areas; 

(i)  A New  Methodology  for  the  Analysis  of  Nonlinear  Systems 

Bifurcation  Analysis  and  Catastrophe  Theory  Methodology  (BACTM) , 
shown  schematically  in  Fig,  2.8,  which  is  based  on  several  important 
results  from  Topology,  Catastrophe  Theory  and  Bifurcation  Analysis,  is 
a powerful  new  tool  for  the  study  of  nonlinear  dynamic  systems.  It 
allows  one  to  study  the  global  behavior  of  nonlinear  systems  in  an 
(n  + m)  dimensional  space  vdiere  n is  the  n\miber  of  state  variables  and 


m is  the  number  of  control  veurlables  in  the  system.  A direct  outcome 
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of  applying  this  methodology  is  a characterization  of  the  stable 
«U)d  unst2d3le  regimes,  stability  boundaries,  domains  of  attraction, 
bifurcation  and  jump  surfaces  and  the  existence  of  limit  cycles.  The 
methodology  can  also  be  used  for  model  reduction  or  simplification 
(using  the  center  numifold  theorem),  feedback  control  system  design 
for  producing  almost  linear  behavior  for  a nonlinear  system,  model 
structure  determination  2md  input  design  for  nonlineeu:  system  identifi- 
cation. It  also  gives  a clear  representation  of  those  cases  where 
linearization  is  valid.  More  importantly,  it  gives  guidelines  for 
extending  the  linezurized  model  so  as  to  mzdce  it  valid  for  those  cases 
where  nonlinearities  zure  significzmt. 

(ii)  Unified  Improved  understanding  of  Aircraft  Stability 
and  Control  at  High  Angle s-of -Attack. 

The  application  of  Catastrophe  Theory  and  Bifvurcation  Analysis 
to  three  different  aircraft  models  has  provided  an  improved  understanding 
of  the  various  junp,  hysteresis  and  limit  cycle  phenomena  observed  at 
high  zmgles-of -attack. 

The  three  aircraft  zu:e  designated  as  aircraft  A,  B and  H.  Their 
physical  chzuracteristics  zure  presented  in  Table  2 of  Chapter  III. 

(See  also  Sec.  3.3).  Aircraft  A is  zm  FIOOA;  aircraft  B is  a small 
mzmeuverable  single-engined  jet,  zutd  aircraft  H is  a swept-wlng  fighter, 
the  F-80A.  Even  though  these  models  do  not  represent  modern  day  fighter 
planes,  their  ia^rtzmce  lies  in  their  extensive  use  in  previous  work 
(see  Chapter  III  for  discussion  zuid  references) . Use  of  these  models, 
therefore,  provides  a solid  basis  for  cosipzurlson  of  the  BACTM  with 
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other  techniques,  especially  analytic  results.  BACTM  has  been 
demonstrated  using  these  models,  and  it  is  currently  being  applied 
to  modern  fighter  plemes. 

It  has  been  shown  for  aircraft  B that  the  bifurcation  and 
catastrophe  surfaces  are  of  the  "butterfly"-type.  The  projection 
of  the  bifurcation  surface  in  the  control-state  space  onto  the  control 
plane  results  in  a pattern  vAiich  somewhat  resembles  the  shape  of  a 
butterfly.  See  Figs.  2.9,  3.10,  and  3.11  for  examples.  See  also 
Fig.  3.13  and  the  discussion  in  Sec.  3. 6. 1.3  relating  to  the  bifurca- 
tion to  an  autorotational  state  for  aircraft  B as  6e  increases  from  0° 
to  12®,  with  6a  = 6r  = 0*  (the  bifurcation  occurs  at  6e  = 9.3“). 

A bifurcation  surface  divides  regions  in  the  control  space  where 
different  numbers  of  equilibrium  states  are  possible.  Thus,  as  controls 
vary  in  such  a way  as  to  cross  the  bifurcation  surface,  catastrophes  in 
the  form  of  sudden  "jumps"  between  equilibrium  solutions  occur.  The 
solution  is  said  to  bifurcate  to  a new  equilibrium  branch  in  state 
space.  The  bifurcation  surface  marks  the  boundeury  between  the  stable 
cuid  unstable  equilibrium  solutions.  See  the  discussion  of  Sec.  2.2, 
cmd  particuleurly  Figs.  2.1  and  2.2  for  a better  picture  of  bifurcation 
phenomena.  Note  that  hysteresis  effects  are  prevalent  where 
bifurcations  occur;  which  introduces  the  possibility  that  control 
recovery  actions,  which  are  effective  in  stable  and/or  lineeu:  regions 
of  the  equilibrium  state  space,  may  be  ineffective  or  actually  enhance 
the  destabilizing  motions,  once  a bifurcation  has  occurred. 

The  physical  mechanism  causing  jumps  in  the  case  of  aircraft  B 
can  be  related  to  roll-coupling  first  observed  by  Phillips  in  1948  [6] . 
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C^u:eful  study  of  time-history  motions  of  aircraft  B emd  H 
indicates  that  yaw  rate  behavior  is  an  indicator  of  an  oncoming 
jump.  This  is  because  of  coupling  effects  between  roll  and  yaw  motions 
which  are  enhanced  by  the  peculieu:  geometry  of  high  performance  jet 
aircraft,  particularly  their  small  moment  of  inertia  about  the  roll 
axis.  See  the  discussion  in  Sec.  3. 6. 1.5  zmd  especially  3. 6. 1.6  for 
more  details.  Our  walysis,  furthermore,  predicts  several  other  effects 
such  as  hysteresis  2md  divergence.  It  also  indicates  coordinated  aileron- 
elevator  maneuvers  that  avoid  jimps.  See  Sec.  3.7  for  a discussion  of 
this  BACTN  application.  Even  the  critical  roll  rates  are  predicted 
more  accurately  since  the  sis^lifylng  assumptions  used  by  Phillips  [1] 
are  avoided. 

We  have  also  examined  the  physical  causes  of  jumps  more  closely 
and  developed  indicators  of  incipient  jump  phenomena  which  can  be 
useful  to  the  pilot  during  flight.  One  promising  indicator  of  catas- 
trophic, or  jxxnp,  behavior  is  total  vehicle  kinetic  energy  and  its 
time  derivative.  When  monitored  diurlng  dynamic  response  studies  of 
aircraft  B 2md  H,  the  kinetic  energy  was  observed  to  es^rlence  sudden 
changes  - in  advance  of  a similar  change  of  any  of  the  state  variables  - 
when  the  controls  had  moved  from  a stable  to  a jump-producing  value. 

This  was  cle^u:ly  the  case  for  the  aircraft  B model  (See  Figures  3.29 
and  3 . 30) ; the  presence  of  oscillatory  terms  in  the  aircraft  H case 
(Figures  3.92  and  3. 93) indicates  that  low-pass  filtering  or  a more 
sophisticated  "jump  alarm"  criterion  than  the  one  which  sufficed  for 
aircraft  B - l.e.,  change  of  sign  in  dE/dt  - may  be  developed. 
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For  aircraft  H,  time  history  analysis  has  revealed  families  of  limit 
cycles  in  a five-dimensional  space  which  have  not  been  obtained 
analytically  in  any  previous  study.  A time  history  plot,  shown 
in  Fig.  3.50,  clearly  shows  limit  cycle  behavior.  This  plot  is  seen 
to  be  the  result  of  behavior  predicted  by  inspection  of  the  related 
equilibrium  surface  plots,  shown  in  Figs.  3.48  and  3.49.  In  these 
plots,  both  limit  cycle  amplitude  and  hysteresis  effects  are  indicated, 
in  addition  to  the  specification  of  the  values  of  control  (6a,  aileron 
deflection)  which  will  cause  limit  cycle  behavior  for  aircraft  H in 
this  flight  condition.  See  also  Figs.  3.87  and  3.88,  which  show  limit 
cycles  for  aircraft  H in  the  state  space  (four  plots  of  p vs.  (q,r,a,3)); 
and  Pig.  3.90,  which  shows  how  a change  in  6e  from  2®  to  -2®  causes  the 
motion  to  leave  the  limit  cycle  eind  settle  (after  some  severe  transient 
behavior)  at  a stable  equilibrium  point.  The  existence  of  limit  cycle 
phenomena  has  been  postulated  previously,  but  the  analysis  has  generally 
been  confined  to  two-dimensional  cases.  The  existence  of  Hopf  Bifiirca- 
tions  in  aircraft  dynamics  has  been  demonstrated  for  the  first  time. 

The  bifurcation  plots  shown  in  Figs,  3,51  and  3.52  are  generated 
explicitly  by  BACTM.  They  are  also  related  geometrically  to  the  family 
of  equilibrium  surfaces  (Figs,  3.48  and  3.49  being  representative 
members) , in  that  the  bifurcation  surfaces  are  the  projection  of  the 
equilibriisn  surfaces  onto  the  control  space.  The  "folds"  in  the  j 

equilibrium  surfaces  cause  regions  of  differing  n\mibers  of  equilibrium 
solutions  to  exist;  the  bifurcation  surfaces  are  the  boundaries  to  these 
regions.  The  physical  mech€misms  causing  these  limit  cycles  2u:e  not 
completely  understood  except  in  special  cases  such  as  wing-rock  and 
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buffeting  (even  in  these  cases,  there  eure  several  physical  explanations) . 

Other  phenomena  that  are  easily  analyzed  using  this  methodology 
eure  "pitch-up",  autorotation  emd  divergence.  Most  of  these  phenomena 
are  not  described  precisely  in  the  literature  and  the  general  conditions 
under  which  they  are  observed  are  not  known  accurately.  Our  approach 
gives  an  exact  qucuititative  representation  of  these  phenomena  and  reveals 
this  dependence  on  aircraft  parameters  emd  control  systems.  It  is 
possible  therefore  to  study  techniques  for  avoiding  these  phenomena 
through  changes  in  design  or  through  control  strategies. 

(ill)  Implications  for  Model  Structure  Determination  and  Controller 


Design  at  High  Angles -of -Attack. 


A number  of  areas  in  v^lch  the  above  improved  understanding 
of  the  high  angle-of-attack  phenomena  cem  be  utilized  have  been 
investigated.  Two  common  problems  in  aircraft  model  identification 
are  model  structure  determination  and  model  validation.  It  is  easily 
Inferred  from  the  Catastrophe  Theory  generic  models  that  the  observed 
bifurcation  behavior  cm  be  used  to  Identify  and  validate  the  structxire 
of  the  model.  This  has  Important  implications  for  input  design  and 
flight  testing.  The  basic  idea  is  that  the  model  structure  should  be 
chosen  in  such  a way  that  it  is  capable  of  matching  the  jun^  and  limit 
cycle  phenomena  observed  in  flight.  At  high  angles-of -attack,  these 
phenomena  are  dcmlncmt  and  their  proper  representation  by  a mathematical 
model  is  essential.  Such  models  are  also  necessary  for  studying 
entry  into  spin  and  for  investigating  recoveries  from  spin  conditions. 

By  studying  different  aircraft  and  by  analyzing  the  effect  of 
feedback,  we  have  shown  how  bifurcation  surfaces  can  be  modified  and 
even  avoided  in  certain  cases.  This  has  important  in^llcatlons  for 
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aircraft  and  controller  design.  The  bifurcation  surfaces  may  be 
used  directly  in  the  design  of  stability  augmentation  systems.  (Refer 
to  Sec.  3.7  for  more  details).  Results  there  indicate  that  plots  of 
the  bifurcation  surfaces  cure  a useful  tool  for  defining  relationships 
among  the  control  variables  which  avoid  regions  vdiere  bifurcations 
occur. 

In  Sec.  3.7,  bifurcation  and  equilibrium  surfaces  are  used  to 
develop  an  aileron-rudder-interconnect  (ARI)  control  augmentation  system 
for  use  in  high  angle-of -attack  lateral-directional  maneuvers.  This 
study  was  done  with  the  aircraft  H model,  and  time  history  comparisons 
of  ARI  vs.  non-7^  (6r  = 0*)  systems  show  clearly  that  ARI  not  only 
enhances  stability  but  improves  performance  characteristics.  The  ARI 
gains  developed  by  BACTM  are  effective  because  of  their  direct  reliance 
on  the  bifurcation  surfaces,  which  contain  important  stability  information. 
A linear  ARI  model  was  used  in  Sec.  3.7  for  illustrative  purposes;  however 
BACTM  allows  for  a more  general,  global  ARI -type  relationship  which 
explicitly  accounts  for  elevon  ($e)  deflections  as  well.  It  is  possible 
to  expand  the  criteria  for  the  bifurcation  loci  themselves  to  develop 
further  criteria  for  the  ARI  surface  in  the  three  dimensional  control 
space.  Results  presented  in  this  report  show  that  the  ARI  gains  plotted 
as  a function  of  equilibrium  angle-of -attack,  as  derived  frcxn  the 
bifurcation  surfaces,  behave  in  much  the  same  way,  with  roughly  equi- 
valent values,  as  the  ARI  gains  derived  by  other,  less-global,  methods 
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(1.5)  Orqemization  of  the  Report 

The  report  is  organized  into  two  main  chapters  (II  and  III) 
supported  by  two  appendices  (A  and  B) . Chapter  II  contains  a 
discussion  of  the  Bifurcation  Analysis  emd  Catastrophe  Theory 
Methodology.  Chapter  III  contains  all  the  aircraft  global  stability 
and  control  results  for  three  example  aircrafts  (A,B  and  H) . The 
study  of  aircraft  H involving  nonlineeur  )cinematic  emd  aerodynamic  terms 
is  most  comprehensive  representing  a complete  description  of  the 
equilibrium  and  bifurcation  surfaces  in  an  eight-dimensional  space 
(p,q,r,a,$,2a,6e,6r).  Conclusions  and  recommendations  are  stated 
in  Chapter  IV.  Appendix  A contains  exact  statements  of  BACTM  theorems 
and  Appendix  B contains  detailed  algebraic  calculations  for  aircraft 
B amd  H.  A list  of  symbols  and  nomenclature  is  included  in  Appendix  C. 
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II.  Bifurcation  Analysis  and  Catastrophe 
Theory  Methodology  (BACTM) * 

OVERVIEW 

In  this  chapter,  we  describe  from  a general  system  theoretic 
viewpoint,  the  basic  principles  and  applications  of  the  Bifurcation 
Analysis  and  Catastrophe  Theory  Methodology  (BACTM) . After  giving 
a brief  historical  survey,  the  terminology  of  BACTM  is  introduced 
through  an  illustrative  example.  The  four  main  theorems  of  BACTM 
cure  discussed  thereafter  in  Section  2.3.  The  exact  statements  of  the 
main  theorems  are  contained  in  Appendix  A.  The  application  of  BACTM 
to  nonlinear  dynamic  systems  is  discussed  in  Section  2.4.  Other  methods 
for  the  analysis  of  nonlinear  systems  are  outlined  in  Section  2.5  and 
their  relationship  to  BACTM  is  duscussed.  The  specific  applications 
of  BACTM  for  Global  Stability  Analysis,  Control  System  Design  and 
Nonlinear  System  Identification,  are  described  in  Section  2.6. 

(2.1)  Historical  Background 

Nonlinear  dynamic  systems  are  extremely  rich  in  their  behavior 
and  are  known  to  exhibit  a wide  variety  of  bifurcation  and  catastrophe 
jvimp  phenomena  [37,  38].  In  recent  yeeurs,  several  developments  in 
differential  topology  have  led  to  a fairly  general  classification  of 
bifurcation  phenomena.  Following  the  publication  of  the  thought -provoking 
book  by  Them  [1] , veurlous  applications  of  Catastrophe  and  Bifurcation 
Theory  by  Zeeman  have  appeared  in  the  literature  [2,3,4].  A number  of 
* 

Certain  readers  may  prefer  to  go  through  Chapter  3 before  reading  this 
chapter.  This  process  would  help  in  associating  various  abstract 
quantities  and  concepts  defined  here  with  specific  aircraft  physical 
vjuriables  and  phenomena. 
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these  applications  have  been  to  social  and  biological  systems, 
tfhere  the  lack  of  a mathematical  model  causes  serious  difficulties  in 
the  application  of  the  theory.  This  has  led  to  seme  recent  criticism 
of  Zeeman's  work  on  Applied  Catastrophe  Theory  [34].  The  engineering 
applications  of  Catastrophe  Theory  and  Bifiircation  Analysis,  on  the 
other  hauid,  are  well  founded  since  they  are  based  on  physical  isodels 
[5,6,7,28,29].  Perhaps  the  BK>st  common  engineering  applications  of 
Bifurcation  Analysis  are  in  the  fields  of  Elastic  Stability  of  Struc- 
tures [39],  Nonlinear  Circuit  Analysis  [37]  and  Chemical  Kinetics 
[28,30,35,40]. 

A study  of  the  bifurcation  phenosiena  in  ordinary  differential 
equations  was  started  Poincare  around  1892  [8]  and  was  extended 
fxirther  by  Andronov  and  Pontryagin  [10  ].  A new  impetus  «nis  provided 
in  recent  years  by  the  theory  of  singularities  of  asnoth  maps  [11,12,13,14]. 
A full  understanding  of  these  results  requires  use  of  advanced  concepts 
from  differential  topology.  Our  approach  here  will  be  to  illustrate 
important  concepts  and  results  through  specific  examples,  relevant  to 
engineering  applications.  We  shall  discuss  four  main  theorems  that 
are  the  cornerstones  of  BACTMt  (1)  the  Center  Manifold  or  Reduction 
Theorem  [7,15];  (ii)  the  main  theorem  of  Catastrophe  Theory  due  to  Thom 
[1]  and  Mather  [14]  j (lii)  the  Hopf  BifurcatiMt  Theorem  [15]  j and 
(iv)  the  Global  Implicit  Function  Theorem  of  Palais  [17]  and  its  exten- 
sions. For  illustrative  purposes,  a simple  example  of  bifurcations  and 
catastrophes  is  presented  first. 


1] 

n 

ij 
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(2.2)  Illustrative  Exeanple;  Cusp  Catastrophe 


Consider  a scalar  diffemetial  equation  model  of  the  form 


- X + c^x  + 


(2.1) 


where  x denotes  state  euid  (c  ,c  ) denote  control  veiriables.  Since 

it  is  a scaleu:  system,  the  stable  equilibrium  points  may  be  represented 

4 2 

as  the  minima  of  the  potential  function  (j)  = x /4  + Cj^  (x  /2)  + 

These  are  shown  in  Fig.  [2.1]  for  different  values  of  C2  and  a fixed 
value  of  c = -3.  It  is  clear  from  Fig.  [2.1]  that  for  |c  | > 2,  only 
one  minimum  exists  whereas  for  [c^l  < 2,  two  minima  and  one  maximum  exist. 
The  maximum  represents  an  unstable  equilibrium  point  so  that  the  system 
can  never  stay  in  that  state.  Suppose  c^  is  vauried  slowly  from  -3  to  3 
in  such  a manner  that  the  system  state  x reaches  an  equilibrium  value 
for  each  c^.  Then  the  state  will  stay  in  the  right  hand  minimum  corre- 
sponding to  Xj^  till  this  minimum  coalesces  with  the  meucimum  at  c^  = 2. 

At  this  point  with  a small  change  in  c^,  the  state  will  jump  catastro- 
phically from  x^  to  x^.  If  c^  is  now  reversed,  the  system  will  stay  in 
state  x^  until  c^  = -2  and  then  jump  to  x^  . Thus  there  is  a hysteresis 
effect  in  that  the  location  of  the  jump’ depends  on  the  past  history  of  c^. 

We  now  draw  a three-dimensional  picture  in  the  space 

(see  Fig.  [2.2  ] ).  The  surface  of  equilibrium  points  is  called  M and 
its  projection  on  the  control  space  is  denoted  by  C.  Notice  that  the 
region  in  control  space  over  which  M is  triple -sheeted  (i.e.,  the  function 
f(x,c)  = 0 has  three  roots)  is  enclosed  by  a cusp.  The  boundary  between 
regions  with  different  number  of  equilibrium  points  is  called  the 


P»JJIP 
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Fig»  2,3;  Cusp  Behavior  for  Cj^<  o (No  Bifurcation) 
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Fig.  2.^;  Cusp  Catastrophe  for  c,>  o 
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Bifurcation  Surface  since  the  equilibriun  solutions  separate  (or 
bifurcate)  frcn  each  other  at  this  surface.  The  equations  for  the 
Bifurcation  Surface  are  easily  derived  by  setting  both  f(x,c)  and 
3f(x,c)/3x  equal  to  zero,  i.e. 


W- 


f(x,c)  = X + Cj^x  + c^  “ 0 


(2.2) 


(2.3) 


In  other  words,  the  coalescing  of  a minimun  emd  a maximum  creates 
and  inflection  point  at  which  the  first  euid  the  second  derivatives 
of  <p(x,c)  vanish  simultaneously. 

Elminating  x from  Eqs.  (2.2)  and  (2.3),  we  obtain  the  equation 
of  the  bifurcation  surface  as 


1 3^12  ^ 


(2.4) 


The  vari^d>le  c^  is  called  a splitting  factor  since  for  c^  < 0,  there 
is  no  catastrophe  (Fig.  [2.3  1 ) and  for  c^  > 0,  there  is  a catastrophe 
(Fig.  [2.41). 

Having  defined  some  of  the  terminology,  we  now  discuss  the  four 
main  theorems  of  BACTM. 

* 

(2.3)  Four  Main  Theorans  of  BACTM 

(2.3.1)  Wm  Center  Manifold  or  Tleduction  Theorem  [18,19,201; 
Consider  a system  described  by  a differential  equation 


X - f(x,c) 


(2.5) 


where  x is  n x 1 vector  of  state  variables  and  c is  m x 1 vector  of 


*The  exact  statesients  of  these  theorems  are  contained  in  Appendix  A. 
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control  variables  and  paraaeters.  Suppose  that  for  c « 0,  the 


systssi  has  an  equilibriiai  point  at  x 0 (this  can  be  achieved  by  a 

df 

translation  in  state  mace).  Assume  that  P > ^ x>0,  representing 

***  c-0 

the  dynamics  of  the  linearized  system,  has  k eigenvalues  on  the 
imaginary  axis  (zeros  and  purely  imaginary  eigenvalues  being  included) . 
Out  of  the  restaining  (n-k)  eigenvalues,  let  j be  in  the  left  half 
plane  md  i be  in  the  right  half  plane.  Also,  divide  the  eigenvectors 
of  F into  three  groups  corresponding  to  the  imaginary  <txis,  left  half 
plane  and  right  half  plane  eigenvalues  respectively.  Arrange  the  eigen- 
vectors of  F corresponding  to  k eigenvalues  on  the  Imaginary  itxis  as 

columns  of  an  n x k laatrix  E . Similarly  define  matrices  of  eigen- 

o 

vectors  E_(nxj)'and  E^(nx£)  corresponding  to  the  left-half  and  right- 
half  eigenvalues  of  F.  Consider,  now,  coordinate  transformations  or 


projections  of  the  state  vector  x onto  the  space  vpaimed  by  the  different 
eigenvectors  as  follows < 


u ■ E^x  (kx  1) 


V - E_x  (j  xl) 


w - E^x  (1  X 1) 


Then  the  Center  Manifold  Theorem  states  that  Eq.  (2.5)  is  locally  (i.e., 
for  small  deviations)  equivalent  to  the  following  system  of  equations: 

u ■ g(u,c) 


V « A V 


(2.6) 


w ■ A^w 

T T 

where  A_  ■ " E^FE^. 


A 

I 
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The  iaportant  thing  to  notice  in  Eq.  (2.6)  is  that  the  equations 
for  V azxl  w are  lineeu:  ^md  that  the  dimension  k is  generally  much 
smaller  than  n.  In  fact,  the  maximum  value  of  k depends  on  the 
dimension  m of  the  control  vector  c.  For  example,  if  m = 1,  only 
the  cases  k « 1 (a  single  zero  eigenvalue)  and  k » 2 (a  pair  or  purely 
imaginary  eigenvalues)  are  obtained  in  a "structurally  stable"  way. 

Other  values  of  k are  not  structurally  stable  since  a small  cheuige  in 

the  dynamics  of  the  system  cem  move  the  extra  eigenvalues  away  from 
the  imagin2u:y  axis.  Notice  further  that  a minimum  of  three  control 

variables  are  required  to  place  two  eigenvalues  of  F simulteuneously 
at  the  origin  (alternatively  to  reduce  the  remk  of  F by  two) . General 
relationships  of  this  type  between  m and  k have  been  given  by  Arnold  [7]. 

A direct  consequence  of  Eq.  (2.6)  is  that  the  bifurcational  study 
of  the  original  nonlinear  system  is  reduced  to  a much  smaller  system 
of  dimension  k.  The  theorem  also  implies  that  any  results  that  are 
obtained  on  the  u-system  may  be  generalized  to  higher  dimensional  systems 
by  the  process  of  building  on  or  "suspension"  [7]  with  linear  systems 
having  eigenvalues  with  non-zero  real  parts.  For  applications,  the  most 
significant  consequence  of  the  theorem  is  perhaps  the  central  Importance 
of  the  eigenvalues  on  the  imaginary  axis.  The  nonlinear  character  of  the 
system  is  essentially  unfolded  through  these  eigenvalues.  Therefore, 
the  determination  of  surfaces  in  the  control  and  parameter  space  where 
one  or  more  eigenvalues  cross  the  imaginary  cixis  is  of  paramount  impor- 
tance in  studying  the  behavior  of  nonlinear  systems.  These  surfaces 
may  be  obtained  ^ulalytically  from  models  or  inferred  e}q>erimentally 
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for  real  systems.  In  BACTM  (see  Fig.  2.8),  the  given  equations 
of  antion  for  a system  (viz.  Eq.  (2.5))  are  lineeurized  euround  each 
equilibrium  point,  either  analytically  or  numerically.  The  eigenvalues 
of  each  lineeurized  system  eure  confuted  to  determine  the  control 
conditions  at  which  the  eigenvalues  cross  the  imaginary  axis.  There 
are  cases  (e.g.  aircraft  B)  for  which  only  one  eigenvalue  crosses  the 
imagineury  axis.  Using  the  Center  Manifold  Theorem,  we  can,  therefore, 
study  their  bifurcation  Isehavior  in  terms  of  a single  state  variable. 

The  Center  Mwifold  Theorem  essentially  generalizes  and  extends  the 
llneulzation  procedure  used  so  camnonly  in  control  applications. 

More  importemtly,  it  reveals  vAien  linearization  would  break  down  and 
how  the  model  may  be  extended  to  capture  the  "essential"  nonlinearities 
of  the  system. 

We  now  discuss  theorems  by  Thom  [1]  and  Hopf  [16]  that  chjuracterize 
the  bifurcatlonal  behavior  for  k « 1,2.  The  classification  for  higher 
values  of  k is  still  lnoanQ>lete,  but  the  phenomena  exhibited  by  the 
above  two  cases  is  so  rich  that  a large  number  of  applications  ue 
already  possible. 

(2.3.2)  Main  Theorem  of  Elementary  Catastrophe  Theory  [1,14]. 

This  theorem  provides  a complete  classification  of  the  bifurca- 
tlonal behavior  of  finite  dimensional  systems  of  gradient-type  for 
m 5.  Gradient-type  systems  are  such  that  their  dynamics  locally 

minimize  a potential  function.  For  such  systems  f(x,c)  in  Eq.  (2.5) 

-3 

can  be  expressed  as  '^^(x,c)  where  ^(x,c)  is  a scaletr  potential 

function.  An  in^rtant  consequence  of  this  assumption  is  that  the 

2 2 

llne^u:  part  of  the  dynamic  viz.  F * 3 ^/3x  is  a syanetrlc  matrix 
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and  therefore  has  only  real  eigenvalues.  It  can  be  shovm  [1,7]  that 
for  m <_  5,  only  the  cases  k = 1 euid  2 ceui  be  obtained  in  a structurally 
stcdsle  way  (or  occur  in  "general  position"  [7]  i.e.,  a slight  change  in 
the  function  (j)(x,c)  or  the  parameters  would  not  eliminate  the  above 
singularities)  whereas  higher  order  singularities  of  F,  i.e.,  k ^3 
do  not  appear  in  control  space  dimension  of  less  than  6 in  a structur- 
ally stable  way.  This  CeUi  also  be  seen  from  the  fact  that  F may  be 
decomposed  into  the  product  of  a non-singular  lower  triangular  matrix 
W with  an  upper  triaingular  matrix  G via  Gaussiaui  elimination.  Then, 
for  a stable  rank  deficiency  of  3 in  F,  a minimum  of  six  elements  in  the 
the  last  three  rows  of  G must  be  zero  [35] , which  requires  having  six 
independently  varying  control  parameters.  In  some  practical  applica- 
tions, the  special  structure  of  F may  lead  to  stronger  results,  allowing 
one  to  narrow  down  the  classification  of  singularities  and  the  resulting 
bifurcations . 

Based  on  the  above  discussion  euid  the  Center  Manifold  Theorem, 
the  local  classification  of  bifurcations  for  all  finite  dimensional 
gradient  type  systems  for  m £ 5 ccUi  be  reduced  to  systems  of  dimension 
one  (k  = 1)  and  two  (k  = 2) . However,  a complete  classification  of 
singularities  depends  fxirther  on  the  nature  of  higher  order  terms 
in  the  Taylor  series  expansion  of  the  function  g(u,c)  (cf.  Eq.  (2.6)). 

The  general  theorem  by  Thom  [1]  is  given  in  Appendix  A.  Vfe  consider 
here  the  case  of  four  control  varied^les. 

Table  1 lists  the  seven  catastrophes  given  by  Thom  [1]  for  the 
case  of  m = 4.  The  first  four  are  called  cuspoids  and  correspond  to 
the  case  k « 1,  whereas  the  last  three  are  called  umbilics  and  correspond 


I 

I 


to  the  case  k * 2 
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Reduced 

Control  Space  State  Space 
Dimension , m Dimension , k 


Generic 

g(x,c) 


-(x^-c) 


-(x  +Cj^x-C2) 

-(X  +Cj^x  +C2X-c^) 

/ 5^  K 2^  , 

(x  +C,X  +C^X  +C-X-C.) 
12  3 4 


’‘l‘*’®l*2’®2 


:x2+CiXi-C3 


2 2^ 

*rv°i*i"°2 

^*1*2^°1*2’°3 


2X1X2+2C1X-C2 

*2'^°2*2"°4 


Description 
Name  of  the 
Catastrophe 


FOLD  (2.7) 


CUSP  (2.8) 


SHALLOW  TAIL  (2.9) 
BUTTERFLY  (2.10) 


HYPERBOLIC 

UMBILIC 


ELLIPTIC  . 

UMBILIC  <2.12) 


PARABOLIC 
UMBILIC  ‘2.13) 


Table  1:  Seven  Elementary  Catastrophes  for  Control  Space  Dimension,  m > 4. 
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The  classification  within  each  of  the  above  two  categories  is 
based  on  the  codimension  of  the  singuleurity,  which  qualitatively  is  an 
indication  of  the  degeneracy  of  the  singularity.  The  exact  definition 
of  the  codimension  of  a singularity  may  be  found  in  ref.  [1],  but 
roughly  speaking,  the  codimension  may  be  inferred  from  the  number  of 
Taylor  series  coefficients  that  vanish  in  the  expansion  of  g(x,c) 
around  a given  equilibrium  point.  In  this  sense,  all  equilibrium 
points  [equivalently  critical  points  of  ^(x,c)]  around  which  the  lin- 
earized system  has  no  zero  eigenvalues  are  singularities  of  codimension 
zero.  For  k = 1,  it  is  clear  that  at  least  one  (linear)  Taylor  coeffi- 
cient must  vanish  in  the  expansion  of  g(u,c)  [cf . Eq.  (2.6)}  but  higher 
order  coefficients  may  also  vanish.  The  dimension  of  the  control  space 
m again  limits  the  number  of  cases  that  are  to  be  considered  since  for 
m « 4 and  k = 1,  only  singularities  of  codimension  d £ 4 ceui  occur  in 
a structurally  stable  way.  Similarly,  for  k « 2,  and  m « 4,  only  three 
umbilic  types  of  singularities  occur  as  shown  in  Table  1.  As  an 
example,  consider  the  cusp  catastrophe  of  section  2.2.  The  function 


g(u,c)  takes  the  form  g(u,c) 


3 

(u  + Cj^U  + c^) . 


8u^ 


At  u * 0,  the  first  t%<o  derivatives  of  g(u,c)  viz 


vanish,  whereas 


u»0 


and 


-6.  Thom  [1]  shows  that  the  above  g (•) 


u=o  au-" 

function  is  representative  of  the  class  of  smooth  functions  for  which 

the  first  two  partial  derivatives  are  zero.  However,  to  achieve  this 

condition  in  a "structurally  stable"  way,  at  least  two  independent 

parameters  or  control  variables  (m>2)  must  be  available.  The  specific 

control  values  can  then  be  obtained  by  simply  setting  the  first  two 
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partial  derivatives  of  g(u,c)  equal  to  zero.  If  only  one  control 
variable  is  available  and  the  first  two  partials  of  g(u,c)  happen 
to  be  zero  for  a particuleu:  value  of  c,  then  a slight  perturbation 
in  the  function  g(u,c)  would  destroy  this  property. 

We  now  discuss  the  genericity  property  of  the  models  listed 
in  Tcd>le  1 for  each  catastrophe.  These  are  the  sinqplest  representations 
of  equivalence  classes  of  dynamic  systems  in  the  sense  that  all  other 
dynamic  systems  in  a given  class  can  be  obtained  by  smooth  deformations 
of  the  generic  ones.  The  exact  definition  of  equivalence  is  rather 
technical,  but  should  be  considered  carefully  in  every  application  of 
the  theorem.  (See  Appendix  A for  this  definition) . 

The  practical  importance  of  the  generic  representation  of  Table  1 
in  applications  derives  frtxn  the  following  facts: 

(i)  Their  global  bifurcation  behaviors  have  been  studied 

extensively  and  are  well  documented  in  the  literature  [1,21]. 
(li)  In  system  Identification  problems,  one  or  more  model 

structures  have  to  be  postulated  based  on  physical  con- 
siderations and  on  the  observed  qualitative  properties  of 
the  system.  The  generic  models  provide  a library  with  well- 
understood  dynamic  properties  and  are,  therefore,  ideal 
ceuididates  for  choosing  model  structures. 

(iii)  The  validity  of  the  generic  models  for  physical  systems  C2ui 
be  extended  by  use  of  the  building  on  or  the  "suspension" 
procedure  discussed  in  Section  2.3.1. 

(iv)  The  low  dimension  of  the  generic  models  m^dces  them  amenable 
for  further  study  in  problems  of  experimental  design  and 
optimal  control  (see  Section  2.6). 
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For  illustrative  purposes,  we  show  diagrauns  of  the  seven  eleaentary 
catastrophes  in  Figs.  2.5  and  2.6.  Only  selected  three  dimensional 
sections  are  sho%m.*  tfe  now  discuss  the  nature  of  equilibrium  surfaces 
and  phase  plane  trajectories  for  the  cuspoids  (cf  Fig.  2.5). 

For  a Cusp  Catastrophe,  the  nuucimum  number  of  equilibrium  points 
(stable  ^md  unstable)  for  a given  value  of  the  controls  c is  3 and 
for  a Butterfly,  the  s^une  number  is  5.  Each  equilibrium  corresponds 
to  a minimum  (stedile) , maxinmm  (unstable)  or  inflexion  point  of  the 
potential  function.  At  the  inflexion  point,  the  second  derivative  of 
the  potential  function  vanishes  and  a bifurcation  occurs.  In  higher 
dimensions  (n  > 1) , the  nature  of  nonzero  eigenvalues  aiodifies  the 
state  space  portrait  (see  Eq.  (2.6)).  For  exaaf>le,  one  can  get  nodes 
(i.e.,  all  eigenvalues  of  the  same  sign),  saddles  (i.e.,  eigenvalues 
of  opposite  signs)  juid  bifurcation  points  (one  or  sore  eigenvalues 
zero) . The  nodes  may  be  further  divided  into  sources  or  sin)(s  depending 
on  whether  they  are  unsted}le  or  stable.  The  doaolns  of  attraction  of 
the  sinics  are  defined  in  terms  of  separatrices  of  the  saddles.  At 
bifurcation  points,  the  detersiinemt  of  the  linearized  dynastic  matrix  F 
varnishes  auid  this  provides  one  method  for  locating  the  bifurcation 
surfaces  without  solving  explicitly  for  eigenvalues  of  F. 

The  discussion  in  this  subsection  (2.3.2)  has  concentrated  on  the 
case  of  real  eigenvalues  crossing  the  imaginary  axis,  tfe  now  consider 
the  case  in  which  a pair  of  complex  eigenvalues  of  F crosses  the 
isuginatry  atxis.  It  is  shown  in  the  next  section  Chat  this  results  in 
bifurcation  to  a limit  cycle  )cnown  as  the  Hopf  Bifurcation,  provided 
certain  conditions  aure  satisfied. 


STATE 
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(2.3.3)  Hopf  Bifurcation 

Qualitatively,  Hopf  Bifurcation  is  very  different  froB  the 
elementary  bifurcations  of  equilibrium  points  considered  so  far. 

For  general  autonomous  nonlineeu:  systems  of  the  nongradient  type, 
equilibrium  solutions  which  are  closed  orbits  can  exist.  Such 
solutions  are  also  called  Limit  Cycles  or  periodic  orbits.  It 
was  shown  by  Hopf  [16]  that  a family  of  periodic  orbits  bifurcates 
from  a simple  equilibrium  point  as  a p2u:ameter  of  the  system  is 
vauried  in  such  a way  that  a complex  pair  of  eigenvalues  of  the 
linearized  system  crosses  the  imaginary  axis  with  nonzero  speed. 
Before  discussing  the  Hopf  Bifurcation  Theorem,  we  present  a simple 
exan^le.  The  reader  is  referred  to  Ref.  [5]  for  further  examples 
of  Hopf  Bifurcations. 

Consider  a two-dimensional  system  in  polar  coordinates  (r,6), 
222  -1  *2 

where  r » x + x_  and  0 * taui  — and  the  system  is  described  by 
X 2 


(2.14) 


0 » -1 
r « r(c-r  ) 

2 

The  equilibrium  solutions  must  satisfy  the  equation  r(c-r  ) > 0 
which  lollies  r * 0 and  r ■ provided  c > 0.  For  c < 0,  only  the 

equilibrium  solution  r « 0 exists  and  the  sign  of  the  linear  term 
coefficient  i.e.  c < 0 indicates  that  r » 0 is  a spiral  sink.  But 
for  c > 0,  r * 0 solution  becomes  unstable  (spiral  source)  and  a new 
solution  r > /c  emerges.  The  latter  corresponds  to  a stable  limit 
cycle  whose  radius  grows  as  /c.  The  point  c > 0 is  weakly  attracting 
and  represents  a Hopf  Bifurcation  point.  (See  Fig.  2.7). 
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Hopf  [16]  generalized  the  above  bifurcation  to  higher  dimensions 
and  gave  conditions  for  the  limit  cycles  to  be  stable.  Further 
generalizations  of  Hopf's  results  are  contained  in  Ref.  [5,6,23]  eufid 
Appendix  A.  Examples  of  Hc^f  Bifurcation  in  five-dimensional  air- 
craft dynamics  at  high  2mgles-of-attack  will  be  presented  in  Chapter 
III. 


Hopf  Bifurcation  has  been  studied  mainly  for  the  case  of  a single 
parameter  or  control  variable  i.e.  m = 1.  Using  the  Centre-manifold 
theorem,  this  case  c2Ln  be  studied  in  two  dimensions  since  k > 2 
corresponding  to  a single  pair  of  complex  eigenvalues  crossing  the 
imaginary  axis.  For  m ^ 2,  a double  pair  of  complex  eigenvalues  or 
a single  pair  plus  a real  eigenvalue  may  cross  the  imaginary  axis 
simultaneously.  The  classification  of  Hopf  Bifurcation  for  these 
cases  is  still  incomplete.  However,  the  cases  discussed  by  Taken  [23] 
and  Arnold  [7]  for  m >=  i,2  are  quite  rich  in  their  behavior  emd  their 
uses  in  Hydrodyneunic  SteUsility  «md  flutter  have  been  reported  [7,24]. 

For  the  case  of  a scalar  parameter  u (mol) , the  Hopf  Bifurcation 
theorem  stated  in  Appendix  A requires  that  the  purely  Imaginary 
eigenvalues  cross  the  lmagin^u:y  axis  with  nonzero  velocity  at  the 


critical  value  of  V i.e.  ^ 

ap 


> 0 where  ReX(P)  denotes  the 
ji»0 


real  p2u:t  of  the  eigenvalue  of  F which  vemishes  at  p » 0.  If  this 
condition  is  not  satisfied  and  an  even  derivative  with  respect  to  p 
of  ReX()i)  is  nonvanishing,  then  depending  upon  its  sign,  either  two 
faziilles  of  limit  cycles  or  none  may  bifurcate  at  )i  « 0.  If  an  odd 
derivative  of  ReX(p)  is  nonzero,  then  at  least  locally  a one  parameter 
f sadly  of  lladt  cycles  will  exist. 


I ■ 
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Another  requiroMnt  of  the  Hopf  Bifurcation  theorea  is  that  the 


rest  of  the  eigenvalues  be  in  the  left  half  plane  bounded  away  froei 
the  ijsaginary  axis.  It  turns  out  that  for  the  aircraft  application 
(Chapter  ZIZ,  aircraft  H)  this  condition  is  satisfied  for  all  the 
cases  of  interest  since  one  is  aainly  interested  in  those  situations 
idtere  stability  is  lost  for  the  first  tisM. 


(2.3.4)  Global  Zaplicit  Function  Theoresi 

Zt  suiy  be  seen  froa  the  discussion  so  far  that  an  iaportant 
problea  in  the  application  of  Bifuraction  Analysis  and  Catastrophe 
nieory  is  the  solutiem  of  a set  of  nonlinear  algebraic  equations  for 
different  values  of  a paraMter  or  control  vector  c.  In  particular, 
the  bifurcation  and  catastrophe  behavior  is  intisuitely  connected  with 
the  fact  that  these  equations  have  sniltiple  solutions  which  bifurcate 
frost  each  other  or  coalesce  as  the  vector  c is  varied  over  its 
adadssible  dosMin.  Zt  is  also  of  interest  to  determine  the  conditions 
under  which  a set  of  nonlinear  algebraic  equations  has  a unique 
solution  since  this  isq>lies  the  lack  of  catastrophic  or  jus^  behavior. 

The  global  iac>licit  function  theorem  of  Palais  [17]  has  been 
used  extensively  In  establishing  the  uniqueness  of  equilibrium 


solutions  in  nonlinear  networks  [25].  Consider  the  equation 

f(x,c)  - 0 


(2.17) 


Palais'  theorem  (or  its  generalisation  by  Kuh  and  Hajj  [26]) 
states  two  conditions  for  uniqueness  of  solutions  vis.  (i)  non- 
vanishing of  the  determinant  of  P ■ )f/9x  for  all  (x,c),  and 


(ii)  "growth  conditions". 


f(x,c)  -*•  - as  ||x||  - 


(2. IB) 


The  importance  of  the  first  condition  is  clear  from  the  discussion 
in  Sections  2.3.1  and  2.3.2  for  systems  with  real  eigenvalues  of  F 
since  the  first  condition  implies  precisely  the  ^d>sence  of  bifur- 
cations. In  Chapter  3,  we  will  use  this  theorem  and  certain  special 
cases  of  it  to  study  the  effect  of  feedback  control  on  bifurcation 
surfaces.  In  particular,  if  one  is  interested  in  control  laws  that 
chemge  the  structure  of  a nonlineju:  system  from  bifurcational  to 
nonbifurcational,  then  the  importance  of  the  global  implicit  function 
theorems  is  obvious.  Notice,  however,  that  the  case  of  complex 
eigenvalues  leading  to  Hopf  Bifurcations  needs  further  study  since 
in  this  case  the  determinant  of  F does  not  V2mish. 

(2.4)  A Unified  Methodology  for  Applications* 

Based  on  the  above  four  theorems,  a fairly  coo^lete  methodology 
has  been  developed  for  studying  the  global  st2d>ility  and  control 
characteristics  of  nonlinear  systems.  This  methodology  con^lements 
and  extends  previous  approaches  to  nonlinear  system  stability  analysis 
such  as  Liapunov  methods  and  two-dimensional  phase  pleine  methods. 

We  describe  this  methodology  here  in  some  detail  since  it  is 
different  from  other  applications  of  Catastrophe  Theory  reported  in 
the  literature.  Fig.  (2.8)  sho«rs  the  four  basic  steps  involved  in 
BACTM  %<hich  are  discussed  below. 

(2.4.1)  Problem  Formulation 

In  this  step,  the  variables  of  interest  In  the  system  are 


divided  into  state  variables  x and  control  vari^d>les  c (this  may 
include  parameters) . The  qualitative  difference  between  the  two 
types  of  vari^d>les  is  that  the  former  are  fast-varying  and  the  latter 
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are  slow-varying.  Next,  the  equations  of  motion  for  the  state 
variables  eu:e  written  down  and  conditions  for  the  equilibrium 
state  are  developed.  If  a potential  function  exists  in  the 
problem,  it  is  identified  emd  related  to  the  equations  of  motion. 

The  admissible  ranges  of  values  for  the  state  and  control  varied^les 
2u:e  specified.  It  is  useful,  though  not  essential,  to  identify 
specific  unstable  2Lnd  jump  behavior  known  to  exist  for  the  system 
under  study. 

(2.4.2)  Cong>utation  of  Equilibrium  Surfaces 

This  step,  in  general,  requires  solution  to  a set  of  nonlinear 
algebraic  equations.  One  may  either  assume  values  for  the  control 
vauriables  and  solve  for  the  state  variables  or  vice  versa.  The 
second  procedure  turns  out  to  be  simpler  in  many  cases.  For  example, 
in  the  aircraft  case,  the  equilibrium  equations  are  generally  linear 
in  the  control  variables  but  nonlinear  in  the  state  variables.  Thus, 
it  is  easier  to  obtain  the  control  values  necessary  to  achieve  a given 
equilibrium  state.  The  numerical  techniques  used  to  solve  nonlinear 
algebraic  equations  may  be  Newton-Raphson,  Quasilinearization  or 
Conjugate  Gradient  methods.  An  efficient  numerical  procedure  is 
necessary  since  equilibrium  states  are  to  be  computed  for  a large 
number  of  control  values.  Furthermore,  for  each  control  value,  it 
is  necessary  to  compute  all  admissible  equilibrium  points.  For  high 
dimensional  systems  with  several  constraints,  fixed-point  techniques 
[36]  may  prove  to  be  quite  useful  for  the  computation  of  equilibrium 
points.  The  above  confutation  becomes  quite  easy  if  the  quilibrium 
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equations  ceui  be  reduced  to  a single  nonlinear  equation  as  will  be 
d^nonstrated  in  the  aircraft  B case  considered  in  Section  3*.  In 
addition,  if  this  equation  is  a polynomial  equation,  then  the 
reduction  to  one  of  the  canonical  models  of  Catastrophe  Theory  is 
possible.  Since  the  global  behavior  of  these  ceinonlcal  models  is 
well-understood,  the  equilibrium  surfaces  emd  the  bifurcation 
surfaces  are  obtained  easily.  In  the  more  general  case,  it  is 
possible  to  employ  Quasiline£u:ization  techniques  to  generate  all 
the  branches  of  equilibrium  solutions  by  varying  one  parameter  at 
a time. 

(2.4.3)  Computation  of  Bifurcation  Surface 

The  surface  in  the  control  parameter  space  of  singular  equilibrium 
points  at  which  the  eigenvalues  of  the  linearized  system  cross  the 
imaginary  aucis  are  called  Bifurcation  Surfaces.  At  these  surfaces, 
equilibrium  solutions  branch  out  emd  either  new  equilibrium  points  or 
periodic  orbits  (limit  cycles)  are  born  or  are  destroyed.  Therefore, 
Bifurcation  Surfaces  represent  jump  surfaces  at  which  the  system  may 
jump  from  one  equilibrium  state  to  another.  In  systems  obeying  the 
Delay  rule  [1],  the  jumps  occur  only  when  the  current  stable  equilibrium 
state  disappears  by  collision  with  an  unstable  equilibrium  state. 

The  cuialytlcal  calculation  of  the  bifurcation  surfaces  is  possible 
only  in  simple  cases.  Therefore,  efficient  numerical  techniques  have 
to  be  devised  to  solve  the  set  of  nonlinear  algebraic  equations  repre- 
senting the  Bifurcation  surfaces.  Thom's  classification  theorem  of 

*The  justification  for  this  reduction  is  provided  by  the  Centre 
Manifold  Theorem  (Section  2.3.1)  since  for  aircraft  B only  a single 
real  eigenvalue  crosses  the  origin  under  all  admissible  control 
variations . 
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Catastrophe  Theory  [1]  and  Taken* s classification  theorem  on  Hopf 
Bifurcations  [41]  can  be  used  to  describe  the  qualitative  behavior 
of  the  bifurcation  surfaces  near  the  singular  equilibrium  points. 

On  the  other  hand,  if  the  bifurcation  surfaces  are  obtained  numer- 
ically or  experimentally,  then  the  classification  theorems  can  be 
used  for  validation  of  the  results. 

The  location  of  singular  equilibrium  points  can  be  done 
simultaneously  with  the  computation  of  the  equilibrium  surface  if 
a Newton-type  algorithm  or  Quasilinearization  is  used.  The  extra 
work  involved  is  a calculation  of  the  eigenvalues  of  the  jacobian 
matrix  or  a check  for  the  location  of  eigenvalues  on  the  imaginary 
axis.  In  fact,  starting  from  one  branch  of  equilibrium  points, 
one  may  generate  other  branches  by  proper  calculation  at  the 
singular  equilibrium  points.  The  bifuraction  to  limit  cycles, 
however,  would  require  special  calculations  such  as  converting  the 
equilibrium  equations  to  polar  coordinates  or  the  construction  of 
Poincare  maps  [5]. 

<2.4.4)  Use  of  Bifurcation  and  Equilibrium  Surfaces 

An  amazing  amount  of  dynamical  information  is  compressed  in  the 
bifurcation  and  equilibrium  surfaces.  The  jump  phenomena,  hysteresis 
effects,  divergence  properties,  limit  cycles  and  domains  of  attraction 
can  all  be  inferred  from  these  diagrams.  It  is  also  possible  to  devise 
control  laws  that  will  move  the  system  from  one  equilibrium  state  to 
another  minimizing  a certain  cost  function.  The  design  of  control 
inputs  to  identify  the  model  structure  can  also  be  based  on  the  above 
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information.  The  chief  characteristic  of  these  surfaces  is  that  they 
contain  global  information  regarding  the  system.  It  may  seem  that 
for  high  dimensional  systems,  an  excessive  2UDOunt  of  conqputation  is 
required  to  obtain  these  surfaces.  On  the  other  tuuid,  the  computations 
involved  are  much  less  them  the  numerical  solution  of  the  equations 
of  motion  for  a large  number  of  control  histories.  The  latter  approach 
is  routinely  used  on  most  practical  systems,  but  can  miss  many 
significant  phenomena  since  the  number  of  cases  to  be  considered  is 
extremely  large.  The  case  of  aircraft  dynamics  at  high  ^mgles-of- 
attack  is  a good  example  of  this  situation. 

(2.5)  Relationship  to  Other  Methods  for  Nonlinear  System  Analysis 
BACTM  may  be  thought  of  as  cui  extension  of  the  classical  phase 
plane  method  from  two  to  higher  dimensions  with  special  emphasis  on 
topological  properties  of  equilibrium  ^uld  bifurcation  surfaces.  Two 
of  the  other  common  methods  for  nonlinear  system  analysis  are  (i) 
Liapunov  methods  and  (ii)  Describing  function  methods  [44].  Liapunov 
methods  (first,  second  and  extensions  by  Lure,  Popov  amd  others  [42]), 
are  useful  for  determining  whether  a given  equilibrium  point  or 
trajectory  is  stable  in  the  sense  of  Liapunov  SteLbility.  The  main 
difficulty  in  the  use  of  these  methods  lies  in  the  construction  of 
Liapunov  Functions.  Takahashi  et  al.  [43]  state  that  the  difficulty 
of  finding  Liapunov  functions  is  of  the  saune  order  as  that  of  finding 
analytical  solutions  to  differential  equations.  Even  if  a suitable 
Liapunov  function  cam  be  found,  its  use  in  computing  the  domain  of 
attraction  and  in  control  system  synthesis  may  be  very  limited.  BACTM, 
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on  the  other  hand,  does  not  suffer  from  2my  of  these  deficiencies 
since  its  use  can  be  reduced  to  straightforward  computational  steps. 
In  a sense  BACTM  is  more  general,  since  Liapunov  theorems  emd 
extensions  can  be  used  as  part  of  BACTM  to  emalyze  the  stability 
of  equilibrium  points.  In  fact,  the  First  Method  of  Liapunov  is 
used  during  step  2.4.4  2Uid  in  those  cases  where  only  a single  real 
eigenvalue  crosses  the  imaginary  etxis  (e.g.  aircraft  B,  Chapter  III), 
the  use  of  BACTM  leads  to  a potential  function  which  can  be  used  to 
construct  a suitable  Liapunov  function. 

The  describing  function  techniques  [44]  are  inherently  limited 
due  to  assumptions  of  harmonic  inputs  and  specific  nonlinearities. 
They  cannot  be  expected  to  provide  the  global  type  of  topological 
information  which  is  necesseury  for  a complete  understemding  of  the 
aircraft  nonlinear  phenomena  at  high  angles-of-attack. 


(2.6)  Applications  of  BACTM 

As  described  earlier,  BACTM  leads  directly  to  em  understanding 
of  the  global  stedsility  behavior  of  nonlinear  systems  as  control  and 
system  parameters  are  varied  globally.  Two  other  applications  of 
BACTM  are  (i)  input  design  and  model  structure  determination  for 
nonlinear  systems  and  (ii)  nonlinear  control  system  design.  These 
are  discussed  below. 


(2.6.1)  Nonlinear  System  Identification  y 

The  identification  of  nonlinear  systems  of  unknoim  structure  is  u 

of  great  practical  in^rtance,  but  woefully  little  has  been  done  in  a 

this  field.  In  Section  2.3,  we  discussed  four  different  theorems,  H 
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all  of  which  have  implications  for  nonlinear  system  identification. 

The  center  manifold  theorem  allows  us  to  reduce  the  dimensionality 
of  the  system  and  provides  guidelines  for  deciding  whether  a non- 
linear model  is  required  and  of  what  dimension.  For  excunple, 
suppose  the  unknown  system  is  perturbed  around  an  equilibrium  point 
euid  its  linearized  dyncunics  is  identified  using  standard  identification 
techniques  [33].  If  the  identified  eigenvalues  and  their  standard 
deviations  indicate  that  all  the  eigenvalues  are  far  removed  from 
the  imaginary  axis,  then  the  identified  linearized  model  can  be 
regarded  as  satisfactory.  On  the  other  hand,  if  the  95%  confidence 
limits  around  the  identified  eigenvalues  enclose  the  imaginary  axis, 
then  a need  for  the  addition  of  nonlinearities  is  indicated.  (If  the 
confidence  limits  are  too  wide,  one  may  repeat  the  experiment  and 
collect  further  data  to  narrow  the  confidence  limits) . The  eigenspace 
corresponding  to  the  eigenvalues  that  are  most  likely  to  cross  the 
imaginary  axis  would  exhibit  the  "essential"  nonlinearity. 


(2.6.2)  Nonlinear  Control  System  Design 

A standard  requirement  in  control  system  design  is  to  place  the 
poles  of  the  system  in  the  left  half  plane.  It  is  clear  from  the 
theorems  of  Section  2.3  that,  if  a control  law  achieves  this  requirement 


for  a nonlinear  system  at  all  points  in  the  state  and  control  space,  and 
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A more  difficult  problem  to  solve  is  to  design  controllers  that 
do  not  allow  the  eigenvalues  of  the  linearized  system  to  cross  the 
imaginary  axis.  The  problem  is  solved  easily  for  the  generic  models 
(2.7)-(2.13)  since  the  behavior  of  these  models  is  well-understood. 

For  excunple,  in  the  cusp  model  Eg.  (2.8),  negative  feedback  on  c^  viz. 
c_  = (-k.x  + k ) has  the  effect  of  changing  c,  to  (c,  + k, ) . If  c,  is 

2 1 O 1111 

negative,  then  by  choosing  k^  > c^  the  bifurcations  are  eliminated. 

The  sane  effect  can  be  achieved  for  the  "Butterfly**  case.  Eg.  (2.10), 

2 4 

by  choosing  c^  * - (kj^x  + k^x  + k^x  + k^)  and  choosing  the  gains 
(k^,  k^,  k^)  in  such  a way  that  the  case  shown  at  the  top  of  Fig.  2.9 
(unigue  eguilibrium  point)  is  achieved.  For  more  general  models,  the 

control  parameters  (c. ...c.)  aure  nonlinear  functions  of  the  physical 

1 4 

control  variables.  In  such  cases,  a feedback  law  on  a physical 
control  has  the  effect  of  changing  more  than  one  control  parameter 
in  the  generic  models.  The  design  of  bifurcation-free  control  laws 
is  more  complicated  in  these  cases.  However,  we  do  have  a clear  cut 
methodology  to  proceed  and  to  study  the  effect  of  different  control 
laws  on  the  bifurcation  behavior  of  the  system. 


C2=  0 

FIG  2,9:  The  butterfly  and  false  butterfly,  partitions  of  the  control 

manifold  In  the^butterf  ly"  cross  sections  (ci<o).  The  vignettes  show  the 
schematic  cross  sections  of  the  bdiavior  set  I at  various  stations 
(for  the  butterfly,  the  solid  lines  represent  stable,  observable  behavior, 
and  the  dotted  lines  represent  unstable,  unobservable  behaviour, 
and  vice  versa  for  the  false  butterfly). 
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CHAPTER  III 

Aircraft  Global  Stability  and 
Control  at  High-Angles -of -Attack 

(3.1)  Overview 

This  chapter  presents  detailed  numerical  results  of  applying 
BACTM  to  three  aircraft  configurations.  A,  B and  H.  (See  Section  3.3 
for  a description  of  the  characteristics  of  each  of  these  aircraft) . 
Section  3.2  contains  a discussion  of  specific  high  angle-of-attack 
phenomena  providing  the  main  motivation  for  the  current  study. 
Aircraft  characteristics  and  equations  of  motion  core  contained  in 
Sections  3.3  and  3.4.  Numerical  procedures  for  the  solution  of  equil- 
ibrium equations  and  for  the  computation  of  bifurcation  surfaces  are 
discussed  in  Section  3.5.  Specific  global  stability  and  control 
results  for  models  of  aircrafts  A,  B and  H in  the  form  of  two- 
dimensional  projections  of  the  equilibrium  surface,  bifurcation 
surfaces,  qualitative  dynamics,  numerical  verifications  and  physical 
explanations  for  jumps  and  limit  cycles  etre  contained  in  Section  3.6. 
(This  section  contains  the  major  numerical  results  of  the  current  study) . 
Certain  control  and  stability  augmentation  aspects  of  the  model  for 
aircraft  B are  discussed  in  Section  3.7.  Figiire  3.1  defines  and  shows 
the  various  forces,  moments,  angular  velocities  and  anguleu:  position 


for  an  aircraft. 


-46- 


SCIENTIFIC  SYSTEMS,  INC. 


(3.2)  High  Angle -of -Attack  Phenomena 

A large  number  of  the  perform2mce  requirements  expected  of  today's 
maneuverable  jet  aircraft  lead  to  flight  conditions  containing  large 
values  of  angle- of- attack  (a).  Generally,  maneuvers  requiring  a rapid 
increase  in  lift,  such  as  rolling  pullouts,  pullups  and  turns,  require 
not  only  a high  a,  but  may  also  require  a sudden,  l2u:ge  value  of  roll  rate 
(p)  to  properly  orient  the  lift  vector.  However,  as  a increases,  many 
destabilizing  effects  associated  with  flow  sep^u:ation  (such  as  stalling) 

£md  diminished  control  effectiveness  become  prominent;  and,  if  the  flight 
condition  is  severe  enough,  these  effects  can  cause  loss  of  control  of 
the  aircraft.  Furthermore,  while  the  streamlining  trends  in  aircraft 
geometry  certainly  enhance  both  the  speed  and  the  capability  to  perform 
certain  kinds  of  meuieuvers,  they  also  tend  to  compound  the  destabilizing 
tendencies  of  high-a  numeuvers.  For  some  aircraft  the  "destabilizing 
tendencies”  at  high  a are  so  severe  that  the  airframe  and/or  control 
svirfaces  cem  brecdc  apart,  partlcul2u:ly  for  such  maneuvers  as  rolling 
pullouts  where  inertial  <md  gyroscopic  coupling  is  also  prevalent. 

The  essential  feature  of  aircraft  behavior  at  high  a is  that  it  is 
nonllnezu:  in  nature.  Predominant  nonlinear  characteristics  eure  junp  emd 
limit  cycle  phenomena.  The  large-disturbance  nature  of  such  motions  implies, 
furthermore,  the  presence  of  significant  coupling  between  the  longitudinal 
and  lateral  modes.  Hence,  separation  of  modes  based  on  small-dlsturbemce 
linearized  theory  is  no  longer  valid.  Given  a significant  roll  rate  p and 
high  emgle -of- attack  a,  the  coupling  phenomenon  produces  high  sideslip 
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angles  (3)  which  may  generate  hinge  moments  surpassing  the  structural 
limits  of  the  vertical  tail.  Whether  or  not  structural  integrity  remains, 
quite  often  this  condition  induces  Jump  and  limit  cycle  behavior  which 
leads  to  departure  and  subsequent  loss  of  control. 

Prominent  features  of  flight  in  high  a regimes  include: 

• Roll  Coupling: 

The  roll  coupling  problem  is  an  old  one  dating  back  to  the  intro- 
duction of  jet  aircraft.  Phillips  [4]  produced  his  results  in  1948, 
and  the  bulk  of  the  main  results,  including  the  high  angle-of -attack 
phenomena,  were  obtained  in  the  middle  fifties  (see  Pinsker  [13]  and 
Rhoades  and  Schuler  [1],  for  example).  Certain  vertical  tail  failures 
were  caused  by  maneuvers  combining  a simultaneous  roll  and  pull-up 
(rolling  pull-out).  Large  sideslip  angles,  generated  by  the  combination 
of  roll  rate  and  high  angle-of-attack,  were  responsible. 

Previous  work  in  this  area,  while  quite  often  very  productive, 
has  been  hampered  by  the  lack  of  adequate  mathematical  tools  for 
dealing  with  the  nonlinear  models  which  are  required  to  simulate 
the  high-a  maneuvers.  Therefore,  the  studies  tend  to  be  locally 
concentrated,  so  that  at  least  quasi-linear  methods  are  applicable. 

The  risk  is,  of  course,  that  it  is  thereby  possible  to  overlook 
regions  in  the  control  space  in  which  "catastrophic"  phenomena  occur. 

For  aircraft  with  a short  wing  span  and  having  most  of  their 
weight  in  the  fuselage,  and  performing  rolling  maneuvers,  their 
"catastrophic"  phenomena  assume  the  form  of  sudden  changes  in  the 
values  of  the  djmamic,  or  state,  variables.  Such  "Jumps"  tend  to 
occur  almost  entirely  without  warning,  with  very  little,  if  any, 
change  in  control  settings.  The  roll  coupling  problem,  then,  does 
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exhiblt  many  of  the  properties  of  those  dynamic  syst«ns  vAiich 
may  be  dealt  with  using  BACTN.  Later  sections  of  this  chapter 
will  show  that  BACTM  is  a useful  tool  for  analyzing  in  a global 
manner  the  nonlinear  behavior  of  aircraft  performing  high-a, 
rolling  maneuvers. 

Low  Rudder  Effectiveness t This  is  a consequence  of  the  vertical 
tail  (hence,  rudder)  descending  further  into  the  wing's  wake 
as  a increases;  it  can  lead  to  nose  slice  (see  below).  Large 
negative  values  of  a also  inhibit  rudder  control  effectiveness. 
Result  is  weakening  of  "vieathercock  stability"  moment  • 

Reduced  Dihedral  Effect;  Related  to  the  reasons  for  reduced 
rudder  effectiveness,  in  that  the  tail  is  in  the  wing’s  wake, 
so  that  its  lift  due  to  sideslip  0 is  reduced. 

Stall;  The  stall  ^mgle-of-attack,  ttg.paT.T/  angle  of 

maximum  usable  lift  at  a given  flight  condition.  Beyond  this 
a,  while  there  may  be  greater  lift  potential,  severe  longitudinal 
and/or  lateral-directional  instability  will  most  likely  arise, 
with  ensuing  loss  of  control. 

Wing  Rock;  Uncomnanded  lateral-directional  motions,  in  the 
form  of  roll  oscillations. 

Nose  Slice;  Uncomnanded  lateral-directional  motion,  in  the 

form  of  yaw  excursions.  Also  known  as  "directional  divergence"  [2]. 


n 


k. 


1 


1 


.1 


1 


II 


'1 


1 


SCIENTIFIC  SYSTEMS,  INC. 


9- 


0 


n 

V j 


fl 


• Departure  t Aircraft  motions  iimnediately  following  loss  of  J 

control.  Many  of  these  motions  will  shortly  be  Identified 
as  "jump  phenomena". 

• Post  Stall  Gyration;  Uncontrollable  motions  2d>out  more  than 
one  body  zucis  following  departure,  e.g.,  a rapid  roll  after 
initial  yaw  divergence  ("roll  departure") . 

• Spin;  Sustained  yaw  rate  r at  emgles-of-attack  greater  than 

a . If  OL  is  positive  the  condition  is  called  "erect  spin"; 

STAXiL 

if  negative,  "inverted  spin". 

• Incipient  Spin;  In  the  initial  stages  of  spin,  there  is  often 
insufficient  balance  of  the  aerodynamic  and  inertial  moments, 
giving  rise  to  an  incipient  spin  condition.  When  these  moments 
are  in  balance,  the  condition  is  called  "developed  spin". 

• Recovery  Roll;  This  is  roll  rate  occuring  in  the  initial 
phase  of  recovery  from  spin  or  departure.  As  a gets  smaller 
during  recovery,  the  residual  values  of  roll  and  yaw  rates 
generate  a "pitch  up"  tendency,  delaying  the  reduction  in  a. 

• Pre-Stall  Buffeting;  For  Ot  less  than  buffeting  intensifies 

somewhat,  and  a "nose-rise"  tendency  develops.  This  buffeting  i 

typically  does  not  provide  adequate  warning  of  2m  oncoming  stall 
condition.  However  the  BACTM  presented  here  is  2d}le  to  indicate 
beforehand  the  flight  condition (s)  under  which  buffeting  type  of 
conditions  will  give  rise  to  "jump"  catastrophes. 

I 
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Finally,  typical  of  the  phenomena  less  directly  due  to  high 
euigle-of-attack , but  nonetheless  amplified  by  it,  is  that  of  adverse 
aileron  yaw,  in  which  the  rolling  motion  of  the  wings  generates  an  a 
differential  along  the  wing.  This  Act  produces  cm  (2d>etting)  lift 
differential  emd  a drag  differential  which  induces  yawing  motion  in 
the  adverse  sense  of  the  rolling  motion. 

Not  all  of  the  eibove  characteristics  of  high  ot  flight  are  hazardous. 
S<»te  cure  merely  bothersome,  reducing  control  effectiveness,  while  others, 
such  as  the  reduced  dihedral  effect,  may  actually  aid  the  pi  lot  in 

reestctblishing  solid  control  of  his  vehicle,  if  he  has  an  understanding 

* 

of  how  to  utilize  it. 

It  will  be  seen  later  that  many  of  these  phentanena  can  be  categorized, 
predicted,  emd,  ultimately,  controlled  by  utilizing  the  BACTM  methodology, 
presented  in  this  report. 

(3.3)  Aircraft  Characteristics 

The  eurcraft  used  in  this  study  are  identified,  respectively,  as 
aircraft  A,  B emd  H.  In  futxire  work,  it  is  hoped  to  apply  the  results  of 
this  study  to  em  aircraft  of  particuleu:  interest  to  the  Navy,  such  as 
the  F-14. 

In  general,  aircraft  A cmd  B are  distinguished  from  aircraft  H in 
that  the  latter  is  modelled  by  nonlinear  aerodynamic  coefficients,  while 
aircraft  A and  B are  modelled  only  with  line£u:  coefficients  (see  Table  2) . 

The  data  avail£d>lei  on  aircraft  A allowed  study  of  the  roll-coupling 
problem  for  zero  values  of  6e  (elevator  angle)  and  fir  (rudder  angle)  only. 

For  "pitch-down"  maneuvers  (see  Etkin,  p.  153),fie  > 0,  negative  a 
converts  to  negative  3 via  the  rolling  motion.  This  madces  the  dihedral 
desteUbilizing,  so  that  reduction  of  the  dihedral  effect  is  desirable. 
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TABLE  2:  CHARACTERISTICS  OF  EXAMPLE  AIRPLANES 

A/C 


m,  slugs 


I^,  slug  - ft 
I , slug  - ft 

y 


. slug  - ft"” 
z 

per  radian 
per  radian 


ye 

I 

Vr 


yp 

c 

niQ, 

C 

mq 


m 


n 


per  radian 

per  radian 

per  radian 

per  radian 

per  radian 

per  radian 

per  radian 

per  radian 

per  radian 

per  radian 

per  radian 

: , per  radian 

*'6a 

, per  radian 
, per  radi£m 
, per  radian 
, per  radian 


^e 


mje 


n 


‘6e 


_"6r 
qr  Ib/ft 


, per  radiem 
2 


A/C  B 


(Etkin  [5],p.447) 


745 

186.3 

10976 

1700 

57100 

12400 

64975 

13600 

-0.28 

-0.081 

0.34 

0. 

0.15 

0. 

3.85 

4.35 

-0.36 

-0.435 

-3.5 

-9.73 

-1.25 

-2.1 

0.057 

0.0218 

-0.095 

-0.0424 

-0.034 

0. 

-0.255 

-0.442 

0.09 

0.0309 

0.06 

-0.081 

-0.044 

-0.24 

0. 

0. 

0. 

0. 

0. 

-1.07 

0. 

0. 

0. 

0. 

197 

297.3 

V,  ft/sec  691  500 

s,  ft^  377  216 

c,  ft  11.3  6. 

b,  ft  36.6  36.0 


AIRCRAFT  H (Ref.  [1]).  Mach  = 0.9;  h = 20,000  ft. 
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Control  derivatives  with  respect  to  both  5a  (aileron  angle)  and  5e 

were  available  for  aircraft  B and  this  allowed  a study  of  the  roll  coupling, 

autorotation  and  divergence  problems.  Aircraft  H,  taken  from  Rhoades  and 

Schuler  [1}  and  Hacker  and  Oprisiu  [6]  contains  significant  nonlinearities 

in  the  aerodynamic  model,  and  also  gyroscopic  effects  due  to  engine  rotation 

, m , m and  m ).  Note  that  "normalized"  aero  coefficients,  such  as  , 
q P r q e 

which  are  used  by  Hacker  and  Oprisiu,  relate  to  the  more  conventional, 
"dimensionless"  coefficients,  say  C.  , via  the  factors  = q s /I  for 

Jtg  IX 

coefficients  Co  , C.  ; or  F_  = F, (b/(2V)),  for  coefficients  such  as  C, 

Jig  A 5a  2 1 ip 

and  Phenomena  such  as  bifurcations  to  limit  cycles  were  predicted  and 

verified  for  aircraft  H.  Aircraft  A,  taken  from  Haddad  [3]  , is  an  FlOOA. 
Aircraft  B is  a small  jet  airplane  taken  from  Etkin  [5].  It  has  been 
used  extensively  in  the  literature  (see  Schy  and  Hannah  [ 7) ) and  has  been 
eidiaustLvely  studied  here,  along  with  aircraft  H.  The  bifurcational  proper- 
ties of  aircraft  H are  studied  here  for  the  first  time.  The  aircraft  A and 
B are  chosen  for  analysis  here  because  of  the  amount  of  study  they  have 
already  been  subject  to,  with  regard  to  roll-coupling  and  similar  phenomena. 
Therefore,  it  is  meaningful  to  compare  the  results  presented  in  this  report 
with  earlier  work.  With  regard  to  aircraft  H,  the  above  is  true,  as  well  as 
the  fact  that  it  is  a more  comprehensive  model  - i.e.,  more  non-linear  coupling 
terms  are  present,  with  values  for  the  coefficients  of  these  terms  specified 
at  the  appropriate  flight  conditions.  These  extra  terms  for  aircraft  H 
also  introduce  limit  cycle  phenomena  which  are  known  to  exist  in 
modern  fighter  aircrafts.  For  aircraft  A and  B nonlinear 
aerodynamic  data  was  not  available  to  permit  similar  analysis. 

The  algorithms  used  have  been  kept  general  so  that  the  effect  of  additional 
nonlinearities  is  easily  studied.  In  this  report,  only  the  effect  of 
varying  the  physical  controls  (6e,5a,5r)  is  8tudle;d,  but  the  same  approach 
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can  be  used  for  studying  the  effects  of  aerodynamic  and  other  aircraft 
parameter  variations  on  the  stability  properties  of  the  aircraft. 


(3.4) 


itions  of  Notion 


The  most  general  set  of  equations  of  motion  used  in  this  study  is 
for  aircraft  H.  The  equations  for  aircraft  A and  B can  be  derived  by 
setting  some  of  the  parameters  in  the  aircraft  H model  to  zero.  Notation 
and  terminology  used  here  is  taken  from  Hacker  cuid  Qprisiu  [6]  and  is 
given  in  Appendix  C.  The  con^lete  system,  then,  is 

p = A. 6 + A X ^5a  + A q + £ r + + £ r/ta+  £ p - i.qr  + £.  6a  + £x  6r 

3 Ci6a  q r 30f  rtt  p 1 oa  Or 

(3.1) 


q = + m q + ~ 


(3.2) 


' = *^3^  ^ "a6a^^  ^ V ^ V ‘ "^3^^  ^ ^ ''6r^' 


(3.3) 


a = q - p3  + + *6e^®  * (9/Y)  (cos  9 cos  ^ - cos  0^) 


(3.4) 


3 = yo3  + p(sin  + ta)  - r cos  + ^6r^*^  ^ (g/Y)  cos  0 sin 


(3.5) 


^ = p + q tan  0 sin  <|»  + r tan  0 cos  $ 


(3.6) 


0 = q cos  ^ - r sin  ^ 


(3.7) 


In  eqs.  (3.1)  to  (3.7)  the  air  density,  speed  and  aerodynamic 
coefficients  cire  assumed  to  be  constant.  These  equations  also  include 
the  usual  assumptions  that  ta  and  6e  axe  measured  from  trim  settings  for 
steady,  straight  flight,  and  that  ussV,  o = a^  + Aas:  w/V,  and  3 ~ v/V. 
Note,  however,  that  the  complete  rotational  coupling  and  gravity  terms 
such  as  etc.,  introduce  coupling  terms  into  the  equations.  While 
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such  terms  are  nonlinear  in  the  variables,  the  coefficients  are  assumed 
to  be  constant  over  the  flight  regimes  of  interest.  It  is  also  noted 
that  the  aircraft  principal  £uces  have  been  used  here  for  convenience 
(one  benefit  being  that  inertia  products  vanish) , as  opposed  to  the 
stability  axes.  The  latter  2ure  defined  as  the  set  of  wind  cuces  at  a cer- 
tain steady-flight  reference  condition,  most  commonly  with  sideslip  angle 
S = 0.  Stability  cUces  are  generally  the  reference  axes  in  which  the  aero- 
dynamic coefficients  (stability  derivatives)  are  derived,  using  wind  tunnel  i 

data;  however,  transformations  to  body  principal  euces  of  these  coefficients 
are  straightforward.  | 
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The  complete  set  of  equations  includes  the  x-axis  force  equation 
and  the  kinematic  equation  for  yaw  angle,  ip.  However,  under  the  assump- 
tion of  constant  speed,  and  noting  that  is  not  an  element  in  any  of  the 
other  equations  save  its  own  kinematic  equation,  the  full  system  reduces  to 
the  system  of  seven  equations  presented  here,  which  hereafter  shall  be 
called  "complete".  In  order  to  successfully  explore  the  high  angle-of-attack 
phenomena  described  in  Section  3.2,  it  is  not  possible  to  introduce  the 
usual  set  of  assumptions  about  the  aircraft  dynamics  which  lead  to  a 
decoupling  of  the  above  equations  into  a longitudinal  and  lateral  set, 
each  of  which  may  be  sepcirately  analyzed.  However,  it  is  often  profitable 
to  neglect  the  gravity  terms,  as  is  done  in  the  derivation  of  the  equilibrium 
surfaces  and  in  the  examination  of  "jump"  phenomena,  which  are  not  seen  to 
depend  appreciably  on  gravity  effects.  When  this  is  done,  the  kinematic 
equations  for  roll  (<j))  and  pitch  (0)  decouple  from  the  five  force/moment 
equations  for  the  states  (p,w,r,a, 3),  leaving  five  differential  equations, 

(3.8) 
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r » ripB  + + n^r  + n^p  + n^pda  - i^pq  + n^^6r  (3.10) 


a ■ q - p3  + z^Aot  + 


a 


& * YoP  + p(sin  a + Ax)  - r cos  a + y . 6a  + y , 6r 
p o 0 6a  6r 


(3.11) 

(3.12) 


The  neglected  weight  effects  have  em  important  influence  on  peak 
responses,  so  that  if  high  accuracy  is  desired,  the  complete  set  should 
be  used  in  generating  time  histories  (Schy  and  Hannah  [7]).  The  qualitative 
properties  of  the  solution  may,  however,  be  studied  from  Eqs.  (3.8)  to 
(3.12).  For  either  set,  the  following  notation  is  used: 


a a a a 


(3.13) 


m m + m* 
q q a 


(3.14) 


- ''z  - 


(3.15) 


i - (I  - I )I 
2 z X y 


(3.16) 


(3.17) 


Also  y * Y(g/MV),  z ■ Z(g/WV),  £ ■ L/I  , m « M/I  , n » N/I  . Dotted 

X y z 


symbols  denote  time  derivatives.  The  aerodynamic  forces  are  y and  z. 


euid  the  aerodynamic  rocments  are  Jt,  m,  imd  n.  The  notation  S.^,  say. 


d£ 


is  stand£u:d  for  the  expression  , and  implies  a linearization  of 


the  contribution  of  sideslip  to  roll  rate.  Such  lineuization,  unfortun- 
ately, is  necessary  due  to  limited  data.  Finally,  the  state  vector  is 


y 
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(p,q,r,a,B,e  ,♦)  or  (p,q,r.a,6).  depending  on  the  set  used,  and  the 
control  vector  is  (6e,5a,6r),  representing  angular  deflection  commands 
for  elevator,  aileron,  and  rudder,  respectively. 

(3.5)  Solution  Procedures 

In  the  actual  implementation  of  the  BACTM  method  (see  Fig.  2.8) 
computer  progreuns  %fere  developed  tfhich  perform  four  basic  tasks: 

1)  evaluationg  the  equilibrium  solutions  for  the  equations  of  motion 
(Section  3.4);  2)  determining  the  stability  of  an  aircraft  at  these 
equilibrium  solutions;  3)  computing  the  bifurcation  surface  (control 
values  at  which  a stable  equilibrium  solution  bifurcates  from  an 
unstable  one);  and  4)  numerical  integration  of  the  equations  of  motion 
with  respect  to  time,  for  specified  control  inputs.  The  last  program 
produces  a calculated  time  history  of  the  behavior  of  an  aircraft,  and 
is  used  as  a verification  of  the  correctness  of  the  results  obtained 
from  the  first  three  steps.  The  BACTM  method,  programmed  as  described 
above,  has  been  applied  to  the  three  aircraft  (A,  B and  H)  described 
in  Section  3.3.  See  Table  2 for  their  physical  characteristics.  The 
programs  use  the  most  general  set  of  equations  of  motion  [Section  3.4 
Eqs.  (3.1-3. 7)]  for  all  the  aircraft,  typically  by  setting  some  of 
the  parameters  to  zero.  In  this  manner,  while  doing  the  time  history 
solutions,  the  influence  of  force  or  moment  terms  of  interest  may  be 
studied,  e.g.,  the  effect  of  gravity  on  the  motions. 

In  finding  the  equilibrium  or  steady-state  solutions  for  the 
aircraft  under  study,  gravity  was  ignored  (see  Section  3. 6. 1.6  and 
3. 6. 2. 6 for  a discussion  of  the  consequences  of  such  an  assumption), 
and  p,  q,  r,  a,  and  § were  set  to  zero  in  Eq.  (3.1)-(3.5).  The  last 
four  of  the  resulting  five  equations  were  solved  for  Ot,  Ot,  q,  and  r 
respectively*  and  the  expressions  qbtained  were  substituted  into 


*A  bar  over  a variable  denotes  its  equilibrium  value. 
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Eq.  (3.1),  giving  rise  to  2ui  11th  degree  polynomial  in  equilibrium  roll 
rate,  p,  with  em  additional  p ^ term.  The  resulting  equation  has  the 
following  form: 

f^ (p, fie, 6r)6a^  + f^ (p,fie,fir)fia^  + f^ <p,fie,fir)fia^  + f ^(p,fie,fir)fia 

+ fg(^,  fie,fir)  = 0.  (3.18) 

Eq.  (3.18)  is  a 4th  degree  polynomial  in  fia.  Solutions  to  this 
equation  may  be  studied  by  fixing  fie  and  fir,  emd  solving  fia  for  a wide 
range  of  values  of  p.  See  Appendix  B for  definitions  of  the  fi. 

4 

Since  the  coefficient  of  the  fia  term  is  very  small  for  the  air- 
craft considered,  a simple  cubic  equation  was  solved  using  Newton's 
* 

method.  In  solving  Eq.  (3.3)  for  a,  it  was  found  that  the  resulting 
expression  has  a denominator  vAiich  is  a function  of  fia  and  p.  Thus 

- _ N (fia, fie, fir,  F)  p.ig, 

D(p,  fia,  fie,  fir) 


i 


No  solution  will  exist  at  any  point  «diere  D is  equal  to  zero.  Such 
points  do  occur  for  some  values  of  psurameters  for  aircraft  H and 
correspond  to  more  than  one  eigenvalue  crossing  the  imagineury  sucis 
simultaneously.  The  study  of  these  singular  points  would  involve 
umbilics  or  the  Hopf  bifurcation. 

When  the  model  for  aircraft  B is  considered  (see  Etlcin  [5]  and  also 
Schy  and  Hamnah  [9]),  the  equations  are  simplified  a great  deal  since  all 
the  nonlineaor  aerodynamic  terms  are  neglected.  Eq.  (3.18)  becomes  linear 
in  (fia, fie, fir)  £md  the  Isifurcation  surfaces  are  obtained  by  solving  two 

Neglecting  the  -i^qr  term  In  Eq.  (3.8)  results  in  a cubic  equation.  On 
the  other  hand,  neglecting  yields  a quadratic  equation  in  fia. 


'] 

1 
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llnezu:  equations,  viz  (see  Schy  (9]  for  details); 


f^{p,6e,6r)  6a  + fg(p,6e,6r)  = 0 


(3.20) 


6a  + = 0 

3p  3p 


(3.21) 


It  is  further  observed  that  Eq.  (3.20)  is  a fifth-degree  polynomial 
in  p so  that  the  most  general  catastrophe  that  can  be  observed  is  of 
"butterfly"  type.  See  Section  2.3.2  and  in  particualr,  Eq.  (2.10). 

This  fact  is  verified  by  solving  Eq.  (3.2)  and  (3.21)  parametrically  in 
p (see  Section  3. 6. 1.3).  An  eigenvalue  emalysis  curound  each  equilibrium 
point  shows  that  no  stable  autonomous  limit  cycles  exist  in  this  case. 

The  above  method  for  determining  the  equilibrium  surfaces  often 
encountered  difficulties  in  that  for  sane  values  of  the  controls  and 


flight  condition  the  derivation  was  invalid,  and  care  had  to  be  taJcen  to 
recognize  these  conditions.  Further,  for  each  specific  application  of  the 
program  to  a different  set  of  equations,  the  program  had  to  be  rewritten. 


after  a long  emd  tedious  algebraic  effort.  As  a result  it  was  decided  to 
solve  the  set  of  nonlinear  equations  numerically  for  one  set  of  control 
values  and  continue  the  solution  for  the  other  set  of  control  values  in 
terms  of  a differential  equation  that  uses  the  length  of  the  arc  along  a 
solution  path  as  the  independent  variable.  The  details  of  this  technique 
will  be  provided  in  a forthcoming  report. 

The  linearized  stability  analysis  was  used  to  determine  the  stability 
of  a given  equilibrium  solution.  If  the  equations  of  motion  are  written  as. 


P ■ f (P»q»r,a,e,  6a,6e,6r) 

P 

q = f (P»q»r,oi,e,  6a,6e,6r) 
B ■ fg  {p»q»r,o,e,  6a,6e,6r) 


(3.22) 
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then  the  lineeurized  system  matrix  F is  defined  as. 


3f 

3f 

3f 

3f 

3f 

— E. 

_E 

_£ 

_J2. 

3p 

3q 

3r 

da 

3B 

(3.23) 


where  the  partials  are  evaluated  at  a given  equilibrium,  or  trim, 
setting  x = (p,  q,  7 , q 6^) . If  the  real  parts  of  all  eigenvalues  of 
F eure  negative,  the  equilibrium  is  st£Ut>le  €md  is  indicated  by  eui  "S" 
on  the  equilibrium  plots.  A real  eigenvalue  which  is  positive  indicates 
an  unstable  equilibrium  solution  and  is  denoted  on  plots  as  "U”. 

Finally,  complex  eigenvalues  with  a positive  real  component  indicate 
an  unstable  focus  with  a possibility  of  a surrounding  limit  cycle, 
indicated  by  "L".  It  should  be  noted  that  since  each  equllibriaa 
solution  has  five  corresponding  eigenvalues  it  is  possible  to  have  a 
combination  of  the  above  conditions,  leading  to  points  denoted  by  "UL”, 
"UUU",  etc. 

Of  pjurtlcular  Importance  is  the  condition  arising  when  a change 
of  control  input  causes  an  eigenvalue  to  change  from  a negative  real 
part  (stable)  to  a positive  real  part  (unstable).  The  set  of  control 
v€d.ues  at  which  such  a change  occurs  (i.e.  at  least  one  eigenvalue  equals 
zero)  constitute  the  bifurcation  surface.  This  surface  is  determined  by 
solving  the  following  set  of  six  equations  obtained  from  Eqs.  (3.1) -(3. 5) 


J 


( 


v] 


and  Eq.  (3.23): 
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f (P»qfr,a, $, 6a, fie, 6r)  = 0 
P 

f (P»q»r,a,$,  6a,6e,6r)  = 0 

q 

f^(P»q»r,a,e,6a,6e,6r)  = 0 

(3.24) 

f (p,q,r,a, 6,6a, fie, fir)  = 0 


fo(P»q»r,a,6*6a*6e,6r)  = 0 
P 

det|F(p,q,r,a*6»6a»6e»6r) I = 0 


In  the  (6a, fie)  plane,  (fir  fixed),  this  means  solving  a set  of  six 
equations  in  seven  unknowns,  producing  a one  dimensional  surface. 

Finally,  the  most  straightforward  task  is  to  numerically 
integrate  the  equations  of  motion  to  produce  time  histories.  This 
is  accomplished  using  a fourth  order  Runge-Kutta  integration  routine. 

The  program  is  written  to  integrate  either  the  set  of  five  equations 
of  motion  (gravity  neglected)  or  the  complete  set  of  seven  equations  with 
grwity  included.  The  presence  of  gravity  couples  Eq.  (3.6)  and  (3.7) 
into  the  Eq.  (3.1)  to  (3.5). 

(3.6)  Discussion  of  Results 

This  section  is  divided  into  two  main  subsections  viz.  section  3.6.1 
on  aircrafts  A and  B and  section  3.6.2  on  aircraft  H.  Each  subsection  is 
further  divided  into  discussions  of  equilibrium  surfaces  (3. 6. 1.1,  3. 6. 2.1), 
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nonlinear  root  locus  (3. 6. 1.2,  3. 6. 2.1),  bifurcation  surfaces  (3. 6. 1.3, 

3. 6. 2. 2,  3. 6. 2. 3),  qualitative  dynamics  (3. 6. 1.4,  3. 6. 2. 4),  numerical 
verifications  (3. 6. 1.5,  3. 6. 2. 5)  physical  explanations  (3. 6. 1.6, 

3. 6. 2. 6).  Figures  3.2  to  3.102,  which  pertain  to  Sections  (3,6.1)  through 
(3.7),  appear  at  the  end  of  Section  (3,7), 

(3.6.1)  Aircrafts  A emd  B 

Figures  (3.2)  through  (3.25)  pertain  to  aircraft  A 2md  B.  We 
describe  them  in  detail  below. 

(3. 6. 1.1)  Equilibrium  Surface 

Fig.  (3.2)  shows  the  equilibrium  roll  rate  p plotted  against  the 
aileron  angle  6a  for  aircraft  A (FlOOA) . The  other  control  variables 
(6e,6r)  are  zero.  A maximum  of  five  steady  states  (3  stable  euid 
2 unstable)  are  possible.  It  is  clear  from  the  figure  that  once  6a 
increases  beyond  +8*  (approximately)  a jump  in  roll  rate  will  occur  since 
the  middle  stable  equilibrium  branch  disappears.  A physical  explanation 
for  this  phenomenon  which  is  related  to  roll  cot^>ling  will  be  given  in 
section  3.6. 

Aircraft  B;  Figs.  (3.3)  and  (3.4)  show  the  behavior  of  equilibrium 
roll  rate  with  respect  to  6a  for  a pitch  elevator  (6e  * -2*)  and  a pitch 
down  elevator  (6e  ■ 2*) . A marked  difference  in  equilibrium  behavior 
is  noticed  for  the  two  cases.  For  the  pitch  up  elevator,  only  one  stable 
(S)  equilibrium  state  exists  for  -20*  ^ 6a  ^ 20®  so  that  for  all  practical 
purposes,  the  aircraft  is  stable  in  pitch  i:®)  except  for  very  high  values 

of  |6a|(>20*).  For  a pitch  down  elevator  (Fig.  3.4),  five  equilibrium 

* 

states  (3  stable,  2 unstable)  exist  in  the  approximate  range  4®<6a<4® 


I 


1 


;] 

a 


and  one  stcible  equilibrium  point  exists  for  | 6a|>  9®.  As  6a  is  increased 
The  exact  critical  values  for  6a  will  be  given  in  section  3. 6. 1.5. 
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from  zero,  the  roll  rate  Increases  almost  linearly  until  6a  % 4°  at 
which  point  a jun^  occurs  to  a value  on  the  upper  equilibrium  branch. 
Figs.  (3.5)-(3.8)  show  the  behavior  of  equilibrium  yaw  rate  r,  sideslip 
angle  angle-of-attack  cT,  jmd  pitch  rate  q versus  aileron  deflection, 
6a.  At  6a  ^ + 4 , similar  jumps  occur  in  other  states  and  the  direction 
and  size  of  the  jumps  can  be  easily  observed  from  the  above  figures.  If 
the  direction  of  6a  is  reversed,  a second  jump  would  be  observed  at 
|6a|  % 9®,  resulting  mainly  in  a reversal  of  the  direction  of  roll  rate. 
Notice  that,  without  changing  6©/  it  is  not  possible  to  bring  the  roll 
rate  back  to  zero  in  a anooth  fashion.  The  other  possibility  is  to  pro- 
duce via  some  control  action  a Icirge  change  in  the  initial  state  of  the 
aircraft  so  that  it  lies  in  the  domain  of  attraction  of  the  inner  roll 
rate  equilibrium  branch.  The  important  point  is  that  after  a jump  has 
occurred,  the  hysteresis  effect  due  to  separated  domains  of  attraction 
prevents  the  aircraft  from  returning  to  prejump  conditions  unless  a 
coordinated  control  action  is  taken  or  an  abrupt  reverse  change  in  the 
aircraft  state  is  made  by  some  other  means.  These  aspects  of  qualitative 
dyncunics  will  be  discussed  further  in  sections  3. 6. 1.4  and  3. 6. 1.5. 

(3.6. 1.2)  Nonlinear  Root  Locus 

Figure  3.9  shows  a nonlinear  root  locus  plot  obtained  by  varying 
the  equilibrium  roll  rate  from  -500  deg/sec  to  500  deg/sec  and  ccm- 
puting  the  eigenvalues  of  the  linearized  system  around  different 
equilibrium  points  for  6e  = 2®.  It  is  noticed  that  a real 
eigenvalue  crosses  the  origin  precisely  at  those  points  where 
bifurcations  occur.  Referring  to  the  Center  Manifold  Theorem,  it  can 
be  seen  that  the  bifurcation  behavior  for  this  case  can  be  studied  in 
terms  of  one  state  variable,  tfe  have  chosen  the  roll  rate  as  the  single 


variable,  but  it  is  possible  to  choose  any  other  state  variable  in  the 
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present  case.  The  reduction  of  the  bifurcation  behavior  to  one  dimension 
is  achieved  easily  for  aircraft  A and  B since  their  model  does  not  contain 
nonlinear  aerodynamic  terms.  As  discussed  in  Appendix  B,  if  the  equilibrium 
equations  are  solved  for  q,  r,  a amd  g in  terms  of  p,  then  the  following 
fifth  order  polynomial  equation  in  p is  obtained. 


(3.25) 


where,  for  aircraft  B, 


= -21.6 

f^  = -326.4  6a 
f^  •=  50.3  6e  + 358.9 

f^  - 5412.6  6a 

4 

« 11752.8  6e  - 1525.9 

5 

f . = -23015  6a 
o 

Notice  that  f is  consteuit,  (f  ,f  ,f  ) depend  linearly  on  6a  emd 
1 2 4 o 

(f_,f_)  depend  linearly  on  6e.  We  C2m  reduce  the  edxjve  equation  to  the 
3 5 

c£uionical  Butterfly  model  of  Catastrophe  Theory  [1]  by  the  following 
change  of  variables: 

b = p - 

Eq(3.26)  is  derived  by  using  the  fact  that  the  only  difference  between 
Eq(3.25)  and  Eq(2.10)  Is  that  the  coefficient  of  the  quartic  term  is 
zero  in  the  latter  case.  When  the  expression  p ■ b + f2/5fj^  is  substituted 
In  Eq(3.25)  and  binomial  series  expansions  are  used,  it  is  found  that 
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l.  i 


( 

u 


the  coefficients  Cj^,  C2>  c^  and  c^  In  Eq(2.10)  correspond  to 


-3f. 


5f, 


3f; 


2f, 


8f; 


c_  = 


25f; 


5f? 


2f . 


f 

= _li  ZL.  + _i  3 2 

°3  125  4 25  3 

1 1 


= 24 
^^4  ' 3125 


J 

,5 


6 

125 


"3^2 


6 ^4^2 


6 ^5^2 


3f 


6f, 


(3.27) 


(3.28) 


(3.29) 


(3.30) 


The  canonical  model  has  the  form: 


u 

0 

0 

I: 


^ i 


1 


b = 6b^  + 4c, b^  + 3c„b^  + 2c,b  + c 
1 2 3 4 


(3.31) 


(3. 6. 1.3)  Bifurcation  Surface 

The  bifurcation  plot  for  Eq.  (3.31)  is  well-knovm  "Butterfly"  which 
was  shown  in  Fig.  2.9  using  c^  and  c^  as  the  control  variables.  The 
bifurcation  plot  and  the  three-dimensional  Catastrophe  diagram  for 
aircraft  B using  (6a, 6e)  as  the  control  variables  is  shown  in  Figs. 

3.10  and  3.11.  Notice  the  interesting  fact  that  if  6a  = 0 and  6e  is 
increased  from  0 to  12®,  one  crosses  point  C corresponding  to  6e  = 9.3® 
at  which  the  equilibrium  state  ( 6a  = 0,  6r  = 0,  p = 0,  r = 0,3  = 0)  becomes 
unstable.  This  is  also  shown  in  Fig.  3.12  which  should  be  con®)ared  with 
Figs.  3.3  emd  3.4  for  6e  = + 2®.  It  is  clear  from  Figs.  3.11  emd 
3.12  that  autorotation  and  departure  would  occur  at  this  value  of  6e. 

Fig.  3.12  shows  that,  at  6e  = 12®,  the  central  branch  has  become  an 
unstable  one,  while  it  is  seen  to  be  stable  in  Figs.  3.3  and  3.4,  where  6e  = + 2® .>| 
An  analytical  criterion  has  been  developed  for  the  determination  of  critical 
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5e  at  which  autorotation  would  occur  and  is  given  in  Appendix  B 
(section  B.2). 

Figure  3.13  shows  the  behavior  of  equilibriiun  roll  rate  with 
respect  to  6e  for  5a  = 0®.  The  presence  of  a bifurcation  point  at 
6e  = 9.3®  is  demonstrated  beautifully , with  two  solutions  branching  off 
in  opposite  directions.  Notice  that  locally  in  the  neighborhood  of 
(p  = 0,  5e  = 9.3®),  there  are  either  3 or  1 equilibrium  solutions,  but 
globally  there  are  either  5 or  3 equilibrivun  solutions.  A local  bifurca- 
tion cuialysis  would  simply  indicate  a loss  of  stability  at  6e  = 9.3®, 
but  the  global  analysis  also  indicates  the  resulting  autorotation  equilibrium 
state  after  the  jump.  In  the  case  of  Pig.  3.13,  the  choice  between  positive 
amd  negative  roll  rates  would  be  almost  random  depending  on  the  sign  of  a 
perturbation  in  roll  rate  at  5e  = 9.3®. 


I 
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(3. 6. 1.4)  Qualitative  Dynamics  for  Aircraft  A and  B 

Having  generated  the  equilibrium  and  bifurcation  surfaces,  it 

is  now  possible  to  derive  from  them  (Figs.  3.2-3.13  , for  example) 

extremely  useful  information  about  the  aircraft  behavior.  Fig.  (3.2) 

shows  a plot  of  the  equilibrium  (or  steady  state)  roll  rate  p for 

o 

aircraft  A as  a function  of  aileron  angle,  with  oe  = 0 . Note  that 

_ * 

there  are  values  of  6a  for  which  five  solutions  of  p exist.  The 
steibility  and  insteibility  properties  of  the  equilibrium  points  are 
also  indicated  in  Fig.  3.2.  Starting  with  a state  in  the  linear  region 
of  the  figure  (i.e.  around  p = 0,  6a  = 0) , small  aileron  inputs  show 
a (linear)  variation  of  p with  6a.  Similar  plots  of  a vs.  6a  and  3 
vs.  6a  show  that  a remains  near  its  trim  setting  and  3 remains  small, 
for  small  6a.  However,  it  is  clear  from  the  figure  that  at  values  of 
about  8*  for  6a,  the  basic  solution  "disappears",  and  a "jump"  occurs 
from  one  equilibrium  point  to  another.  The  well  documented  and  much 
utilized  criterion  of  Phillips  [8]  gives  a much  higher  value  of  critical 
6a  of  <ibout  15®.  Time  histories  of  p vs.  6a  show  that  the  tendency  is 
for  the  state  to  jump  to  an  "attractor"  state  - in  this  case,  a jump  to 
the  upper  "S"  line,  (solution  branch  (2))  which  is  a stable  equilibrium 
manifold,  rather  than  the  "u"  unstable  solution.  However,  the  new 
stable  solution  represents  a large  change  in  the  value  of  p.  Furthermore, 
the  other  states  will  experience  jumps,  for  example  a from  a small  positive 
to  large  negative  value. 

*The  presence  of  the  qr  term  in  the  set  of  equilibrium  equations  would 
generate  a ninth  order  polynomial  in  p,  gravity  neglected.  See  Section  3.6.2. 
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Other  interesting  features  of  plots  like  Fig.  3.2  will  be  discussed 
below  in  conjunction  with  Etkin's  aircraft  B [5] , which  has  received  much 
more  extensive  cinalysis  here. 

For  aircraft  B,  the  elevator  is  first  set  at  either  +2®,  to 
correspond  closely  with  the  studies  done  eeurlier,  e.g.  Etkin  [5] 
eind  Schy  cuid  Hannah  [ 9 ] . When  6e  =-2®  (pitch  up)  a unic[ue  equilibrium 
state  exists  for  a large  range  of  6a  values  20®),  as  seen  in  Fig.  (3.3). 

A large  cheinge  in  the  shape  of  this  curve  is  noted  in  Figure  (3.4), 
which  corresponds  to  6e  = + 2®  (pitch-down).  Here,  five  possible  equilibrium 
values  of  p exist  around  6a  = 0®,  three  of  which  are  stable.  This  lack 
of  symmetry  about  6e  is  a result  of  the  nonlinearity  of  the  equations. 

Again,  there  are  no  gravity  terms,  and  the  (-i^  qr)  term  in  the  roll 

moment  equation,  Eq.  (3.8)  is  neglected  in  generating  the  curves  of  Fig.  (3.3) 

* 

through  Fig.  (3.4).  The  equilibrium  solutions  of  yaw  rate,  sideslip, 
euigle-of-attack , and  pitch  rate  cure  shown,  respectively,  in  Figs.  (3. 5, 3. 6, 

3. 7,3.8)*  In  all  these  figures,  note  that  there  are  two  sets  of  equilibrium 
states  at  which  jumps  may  occur,  6a  = ^4®  and  ^ 9®.  For  example,  the 
motion  could  begin  for  6a  = 0*  ( 6 e is  held  at  +2®);  as  |6a|  moves  beyond 
4®,  the  equilibrium  values  of  the  state  will  change  suddenly  to  new,  larger 
values.  The  yaw-rate  r undergoing  a large  change  in  particular  (nose  slice) 
whil0  a chcinging  sign  frcrni  negative  to  positive,  are  two  of  the  most 
drastic  effects,  although  all  five  states  e:q>erience  junps,  phenomena 
closely  related  to  departure  and  post  stall  gyration  phenomena  described 
in  Rutan,  et  al.  [10].  From  this  initial  jump,  the  equilibrium  states 

The  effect  of  the  qr  term  is  found  to  be  negligible,  as  will  be 
demonstrated  for  the  model  of  aircraft  H. 
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will  vary  smoothly  with  increasing  |6a|,  even  though  the  "equlibrium" 

state,  with  its  high  values,  is  hardly  a physically  desirable  situation, 

being  characterized  by  severe  coupling  on  all  axes.  If  |6a|  is  then 

reduced  to  0®,  it  will  be  noted  that  there  is  no  smooth  mechanism  to 

return  the  equilibrium  solutions  to  their  original  linear,  small-value 

* 

condition  by  using  6a  alone.  That  is,  one  expects  a hysteresis  condition, 
resulting  in  an  autorotational  state  (sustained  rolling  motion  for  6a  = 0®) . 
Besults  presented  in  Sec.  3. 6. 1.5  will  indeed  confrm  this  and  other 
observations . 

An  asymptotically  stable  autorotational  steady  state  is  actually 
dangerous  if  the  stationary  value  of  at  least  one  other  variable  than  p 
is  high.  This  would  then  induce  either  large  angles  of  attack  or  sideslip 
or  high  rates  of  pitch  or  yaw.  To  continue  changing  |6a|  in  the  same 
direction,  increasing  it  to  9®,  will  initiate  the  second  jump  condition. 

The  change  in  values  for  this  jvanp  are  even  more  severe  than  for  the 
inner  jump  ( |6a|  = 4®). 

Figure  (3.9)  shows  a locus  of  the  eigenvalues  of  the  linearized 
system,  plotted  as  a function  of  p.  Critical  roll  rates  (i.e,  where 
jumps  occur)  correspond  to  values  at  which  an  eigenvalue  crosses  the 
imaginary  cucis.  Note  especially  that  only  a real  eigenvalue  moves  over 
to  the  right  half  plane  in  Fig.  (3.9).  Therefore,  only  elementary 
catastrophes,  i.e.,  jumps  from  one  equilibrium  point  to  another,  are 
observed.  More  generally  (see  discussion  of  aircraft  H in  Sec.  3.6.2), 
bifurcations  to  limit  cycles  may  occur  when  a complex  pair  of  eigenvalues 
crosses  the  imaginary  cocis. 
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The  bifurcation  plot,  shown  in  Fig.  (3.10)  for  aircraft  B in 
the  (5a, 5e)  control  plwe  Indicates  that  instability  is  possible  for 
5a  = 0 cuid  5e  increased  beyond  9".  The  ph^onenon  here  is  auto- 
rotation and  departure,  and  Figs.  3.11  and  3.12  aroplify  this  prediction. 
Based  on  an  analytical  criterion,  the  5e  at  which  departure  occurs  in 
this  case  is  9.3°,  which  corresponds  to  an  0(  of  -16.3°.  The  analytical 
criterion  is  very  ea.>y  to  apply  and  involves  the  solution  of  at  most 
a quadratic  equation  in  5e.  (See  Appendix  B) 

(3.6. 1.5)  Numerical  Verification. 

Time  history  plots  obtained  by  the  integration  of  differential 

equations  (3.8)  to  (3.12)  for  aircraft  B are  shown  in  Fig.  (3.14).  Also 

given  are  the  (5a, 5e)  control  histories  (this  model  does  not  contain 

rudder  derivatives  so  that  5r  = 0 is  implied) . The  state  variaUbles  have 

initial  values  of  zero.  Over  the  first  20  seconds,  5e  was  fixed  at  2° 

(pitch  down)  and  6a  varied  in  steps  from  0°  to  10°,  and  back  to  0°.  These 

inputs  were  chosen  to  verify  the  jump  and  hysteresis  phenomena  predicted 

from  Figs.  (3.4)  to  (3.8).  Note  that  the  steady  state  variables  jump 

shortly  after  5a  changes  to  10°  and  reach  values  which  closely  correspond 

to  those  predicted  by  the  equilibrium  surface  plots  of  Figs.  (3.4)  to 

(3.8).  Furthermore,  the  return  of  5a  to  0°  at  t = 10  seconds  does  not 

return  p to  zero,  but  results  in  a rather  severe  autorotational  condition, 

a cle^u:  exan^le  of  hysteresis.  Changing  5e  properly  cam  induce  recovery 

from  this  autorotation  involving  the  roll  axis.  The  strategy  for  5e  may 
_ 

The  numbers  on  the  figure  indicate  the  number  of  possible  equilibrium 
solutions  in  that  region. 
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be  derived  from  Fig.  (3.10);  i.e.,  to  move  the  control  point  out  of  the 
Butterfly  pocket,  where  the  states  are  multi-valued,  to  the  single-valued 
region  of  small  motions,  set  6e  ^ 0°.  This  was  done,  at  t = 20  seconds, 
in  Pig.  (3.14),  and  the  varieUDles  do  indeed  jump  toward  the  initial  state 
at  t = 0.  It  should  be  obvious  that  the  (global)  surface  of  Fig.  (3.11) 
can  be  used  as  a starting  point  to  derive  control  strategies  which  will 
transfer  the  aircraft  from  one  equilibrium  state  to  another  without  jun^s. 

In  order  to  investigate  more  closely  the  physical  nature  of  the 
jvimp  phenomenon  (see  next  section),  several  time  histories  were  run  in 
the  "inner  jump"  region  ( | 6 a | 4**)  of  the  well-treated  pitch  down  case 

(6e  = 2®),  and  some  of  these  are  shown  in  Figs.  (3.15)  to  (3.24).  For 
all  of  these  cases,  the  initial  value  of  the  state  variables  is  zero. 

Figure  (3.15)  shows  how  the  state  variables  respond  to  step 
variations  in  6a,  for  a fixed  6e  (*2®).  Each  setting  of  6a  is  held  for 
3 seconds,  to  allow  the  system  to  reach  its  equilibrium  state.  As  can 
be  seen,  there  is  a definite  jump,  occurring  betwen  t = 6 and  t = 9 sec. 

The  jump  is  especially  sharp  for  the  variables  r and  a,  and  in  all  cases 
these  time  histories  reflect  exceptionally  well  the  behavior  predicted 
by  their  respective  equilibrium  curves  (Figs.  [3.4]  to  [3.  8 ]).  Note  how 
3,  q and  r decrease  in  magnitude  after  the  juirp,  and  note  how  a changes 
sign  from  negative  to  positive,  as  Fig.  (3.7)  predicts.  (This  particular 
feature  must  be  especially  disconcerting  to  the  unsuspecting  pilot,  who 
initially  comnanded  a pitch  down  mode) . There  is  no  evidence  either  in 
the  equilibrium  plots  or  in  the  time  histories  of  self-induced  buffeting 
phencanena  or  wing  rucK  (a  more  con^Iicated  model,  A/C  H,  does  exhibit  these 
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cind  other  lioiit  cycle  motions  - see  Sec.  3.6.2);  however.  Fig.  (3.15) 
does  evidence  many  of  the  phenomena  typical  of  the  high  cuigle-of-attaclc 
regime  descriJsed  in  Sec,  (3.2)  and  extensively  elsewhere,  including 
reduced  rudder  effectiveness,  nose  slice  emd  departure.  Noteworthy  is 
the  observation  that  it  is  possible  to  exploit  the  excellent  correspondence 
between  the  equilibrium  curves  cind  the  time  histories  by  using  the 
equilibrium  plots,  which  have  the  advemtage  ]x>th  of  being  relatively 
inexpensive  to  generate  emd  of  being  global  in  nature.  It  is  also 
possible  to  identify  the  flight  conditions  and  control  settings  which 
lead  to  harmful  high-a  motions. 

The  yaw  rate  curve  of  Fig.  (3.15)  seems  particularly  sensitive 
to  incipient  jump  conditions  (the  reasons  for  this  will  be  given  in 
the  next  section),  and  it  provides  a clue  that,  while  the  jump  has 
cleeurly  been  Initiated  by  the  time  6a  = 6®,  nevertheless  an  inspection 
of  smaller  values  of  6a  might  prove  fruitful. 

Figures  (3.16)  to  (3.18)  each  show  time  history  plots  for  one 
setting  of  6a  per  solution,  again  with  zero  initial  conditions.  For 
6a  3®  (Fig.  [3.16]),  there  is  no  juit^,  although  it  takes  upwards  of 
8 to  9 sec.  to  bring  the  system  to  steady  state.  Then,  as  is  quite 
accurately  predicted  by  the  corresponding  equilibrium  curves  (Figs. 

[3.4  ] to  [3.8  ]),  a jun^  occurs  for  6a  = 4®  at  about  6 sec.  into  the 
solution  (Fig.[3w  17] ) . Yaw. rate  is  the  cleanest,  eind  also  the  earliest, 
indicator  of  the  jun®)  situation.  Pitch  rate  experiences  a marked 


increase  in  anplitude  as  early  as  1.5  sec;  however,  the  seune  thing  happens 
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in  the  (stcible)  6a  = 3®  case,  except  that  q manages  to  return  to  a 
value  in  the  steady  state  somewhat  close  to  its  "first  plateau"  value 
prior  to  t = 1.5  sec.  The  mechanism,  then,  influencing  this  pre-jump 
behavior  in  q is  not  related  to  conditions  causing  the  jump,  but  is 
a result  of  the  (nonlinear)  lateral- longitudinal  coupling  in  the  equations. 
In  the  jiamp  case,  yaw  rate  changes  sign  (neglecting  the  cross-over 
at  cibout  1.5  sec.  due  to  Dutch  roll  transients),  upsetting  the  moment 
balance  in  the  pitch  moment  equation,  which  forces  q to  keep  growing. 

Since  r manages  to  stay  positive  at  5a  = 3®,  the  coupling  influence 
on  q is  a stcibilizing  one.  Similarly,  none  of  the  jump  changes  in  the 
other  varicibles  precedes  the  change  in  r. 

Finally,  the  case  6a  = 5®  is  presented  in  Fig.  (3.18)  to  show 
that  a more  extreme  junp  setting  will  cause  the  jump  to  begin  sooner. 

The  jump  does  not  come  to  con^Jletion  noticeably  quicker,  however. 

The  above  observations  indicate  strongly  that  the  region  in 
control  space  for  which  a jun®>  will  occur  is  separated  from  the  stable 
region  by  a distinct  boundary  (in  fact,  the  bifurcation  line  of  Fig. 
[3.10]).  There  is  no  intermediate  region  in  which  a juiip  may  or  may 
not  occur.  One  can  expect,  then,  the  highly  nonlinear  situation  of  one 
value  of  5a  producing  well-behaved,  stable,  motions,  yet  a value 
(6a  + e)  causing  a juit®>.  Bifurcation  Theory  implies  this,  the  eibove 
observations  support  the  contention,  and  the  plots  of  Figs.  (3.19)  - 
(3.22)  tend  to  verify  it.  For  exan^le,  6a  = 3.7®  (Fig.  [3.20])  is  clearly 
stable  while  6a  = 3.9®  is  not  (Fig.  [3.22]).  As  before,  the  initial 
value  of  all  state  varicibles  is  zero.  From  the  data  which  generated  the 
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blfurcatlon  surface  of  Fig.  (3.10),  the  critical  value  of  6a  Is 
given  as  3.828**.  This  Is  In  close  agreement  with  the  results  presented 
below,  considering  the  effect  of  nonlinearities  in  Liapunov's  First  Method 
would  slightly  reduce  the  stability  limits. 

While  the  examination  of  time  history  plots  adds  much  insight  to 
the  understanding  of  jump  phenomena,  there  are  pitfalls;  it  is  possible 
to  overlook  conditions  which  indicate  em  incipient  jump.  This  is 
illustrated  by  the  6a  = 3.8®  case,  Fig.  (3.21).  Through  t = 12  sec., 
everything  seems  to  be  well-behaved,  although  a sharp-eyed  observer  might 
detect  that,  unlike  the  6a  » 3.6,  3.7®  cases,  a steady  state  has  not  been 
reached.  Compare  in  particular  the  A and  R (a  euid  r)  plots  in  each  of 
Figs.  3.20  (6a  = 3.7®)  and  3.21  (6a  = 3.8®).  Note  that  they  have  achieved 
a steady  state  condition  in  the  3.7®  case  while  there  is  (barelyl)  noticeable 
oscillation  still  occurring  for  6a  * 3.8®.  Note  also  in  the  3.8®  case  that 
the  Q and  B (q  and  3)  plots  continue  to  move  away  from  the  origin,  which 
is  not  the  case  when  6a  = 3.7®.  Further,  r has  gone  negative  at  around 
10.5  sec.,  and  the  previous  discussion  indicates,  therefore,  that  the 
problem  may  have  reached  an  incipient  jump  condition. 

Therefore,  the  6a  = 3.8®  solution  was  run  for  24  seconds  and 
Fig.  (3.23)  shows  clearly  that  a jump  has  occurred.  Besides  showing 
that  as  little  as  a 0.1®  change  in  6a  c^m  produce  significantly  greater 
changes  in  the  aircraft's  behavior,  due  to  bifurcation  phenomena,  there 
is  also  manifested  a definite  "caveat"  against  relying  too  heavily  on  time 
histories  to  define  the  stability  behavior.  The  expense  of  generating  these 
runs  (compounded,  of  course,  as  the  model  beccxnes  more  complex)  and  the 
fact  that  only  a point-by-point  examination  over  the  control  space  is 
possible,  preclude  exclusive  use  of  the  time  histories  for  a definitive 
evaluation  of  high-a  performance,  even  though  useful  insight  may  be  derived 
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from  solutions  taken  at  selected  points  of  interest.  The  problem  with 
using  time  histories  to  define  the  globa'i.  bifurcation  surface  is 
compounded  further  by  two  facts:  (1)  the  .learch  procedure  is  iterative 
in  nature,  since  one  is  wise  to  initiate  it  with  a "coarser”  grid  of 
control  settings  (Figs.  [3.16J  to  [3.18]),  and  then  proceed  to  a "finer" 
grid  (Figs.  [3.19]  to  [3.22]);  (2)  as  is  evident  by  the  6a  = 3.8®  example, 
the  closer  one  is  to  the  bifurcation  surface,  the  longer  it  takes  for  the 
jixap  to  initiate;  which,  if  nothing  else,  adds  to  the  computational  require- 
ments. 


The  next  section  attempts  to  explain  what  is  happening  physically 
during  a jump,  and  will  report  on  initial  investigations  into  a formulation 
whose  Implementation  has  the  potential  of  providing  the  pilot  (or  auto- 
pilot) with  on-board  real-time  warning  of  incipient  jump. 

(3. 6. 1.6)  Physical  Explanation. 

In  the  previous  section,  predictions  derived  from  an  inspection  of 
the  global  equilibrivan  and  bifurcation  surfaces,  regcurdipg  the  behavior 
of  a jet  aircraft  (B)  in  high  angle-of-attack  regimes,  were  seen  to  be 
readily  verified  by  integrating  the  systen  equations  (3.8)  to  (3.12). 

The  detail  of  these  verifications  resulted  in  some  further  insight  as  to 
the  mechanism  causing  a jump  condition,  and  in  Sec.  3. 6. 1.5  some  of 
the  observed  kinematics  were  described  - specifically:  (1)  the  sensitivity 
of  yaw  rate  r to  jvmp  conditions  ("jump  conditions"  being  the  change  in 
the  control  settings  for  a particular  equilibriiim  solution  in  a direction 
which  would  cause  the  equilibrium  solution  to  "junp"  to  a different  branch. 
See,  for  example.  Fig.  3.4,  where  a jump  condition  exists  if  fia  increases 
frcci  some  value  less  than  about  4 ® to  a value  greater  than  4®).  (2)  the 

ii^iortance  of  coupling  as  a aechanisa  for  either  stabilizing  or  destabilizing 
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the  motion;  (3)  the  m^ulner  in  which  the  coupling  roectumism  and  the 

change  of  sign  in  r interact  to  bring  about  incipient  jump  (euid  nose  slice) . 

The  computer  runs  which  generated  the  time  history  plots  presented 
in  the  previous  section  were  modified  to  provide  time  history  plots  of 
roll  and  yaw  moment  tejrms  which  comprise  the  right-hand  side  of 
equations  (3.8)  to  (3.12).  Note  the  ed^sence  of  gravity  terms  in  these 
equations.  The  terms  are  clustered  into  groups  for  plotting,  as  follows 
(Notation  is  for  aircraft  H,  which  is  a more  complete  formulation) : 

Roll  (Eq.  3.8) 

Aerodynamic  terms 

A « i-B  + £,  p + J,  r 
3 P r 

Control  terms 

C - A^^fia  + A^^fir 

Inertia  term 
I - -ij^qr 

Nonlinear  terms 

Yaw  (Eq.  3.10) 

Aerodynamic  terms 

A - npB  + n^r  + n^p 

Control  terms 
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Inertia  terms 


I = -i3Pq 


Nonlinear  terms 


N = n . a6a  + N ^pa 
aoa  pot 


Each  of  these  four  terms,  together  with  its  respective  moment 

(p  or  r) , is  plotted  against  time.  ([p  or  r]  = A+C+I+N)  It  should  be 

noted  that  N=0  in  both  roll  and  yaw  for  aircraft  B,  These  terms,  however, 

do  appear  in  the  results  for  aircraft  H.  A typical  plot,  the  roll  equation, 

is  shown  in  Fig.  3.24.  This  plot  represents  the  case:  aircraft  B, 

6a  = 3.7®,  6e  = 2®,  6r  = 0®;  the  companion  state  variable  time  history 

is  given  in  Fig.  3.20.  In  Fig.  3.24,  all  terms  are  converted  to  units 
2 

of  deg/sec  . Also,  the  roll  axis  I-term  is  scaled  by  a factor  of  100  in 
order  to  detect  its  behavior.  This  term  has  little  effect  on  the  presence 
cuid  magnitude  of  jumps,  obviously,  but  is  presented  for  completeness. 

A run  made  with  produced  essentially  the  same  results,  as  discussed 

below.  The  yaw  axis  I-term,  however,  is  most  emphatically  not  insignif- 
icant. 

An  inspection  of  a set  of  these  plots  for  aircraft  B in  the  flight 

I 

conditions  depicted  by  Figs.  3.20  and  3,23  (the  6a=3.7®  and  3,8®  cases 
respectively)  leads  to  the  following  conclusion  (Figs,  3,24  to  3,27) : 

The  mechanism  of  the  jump  is  related  directly  to  the  inertia 
properties  of  the  aircraft,  the  particular  flight  condition  and  the 

1 

degree  of  control  surface  deflection.  As  Rhoades  and  Schuler  [1]  aptly 
state  it,  the  "real  villains"  are  the  (i2Pr)  and  (-i^pq)  inertia  coupling 
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terms  which  appear,  respectively,  in  the  pitch  and  yaw  moment  equations. 

To  a certain  degree,  the  relationship  between  inertia  coupling  emd  jump 
phenomena  has  also  received  widespread  recognition  among  other  investigators 
(e.g.,  Phillips  [4],  Etkin[5],  Pinsker(13],  and  Schy  and  Hannah  [9]). 

The  flight  condition  and  degree  of  control  surface  deflection  are  also 
important  factors,  however.  Along  with  the  vehicle  geometry,  these  factors 
determine  whether  the  state  of  the  vehicle  is  in  a region  of  incipient  jump, 
or  whether  it  evolves  to  a stedale  steady  state. 

Further,  for  the  specific  case  of  aircraft  B initiating  a roll 
maneuver  in  a pitch-down  elevator  setting,  the  lateral-directional 
stability  is  seen  to  be  the  most  fragile.  This  means  that  the  (-i^pq) 
term  in  the  yaw  rate  equation  is  the  chief  mechanism  for  inducing  the  jump. 
The  following  observations  with  regard  to  this  flight  situation  are  supple- 
mented by  referencing  the  relevant  equilibrium  and  bifurcation  surfaces,  as 
presented  in  Sections  3. 6. 1.1  euid  3.6.1. 3,  and  as  discussed  there  and  in 
Sections  3. 6. 1.4  and  3. 6. 1.5. 

The  yaw  rate  inertia  coupling  term  (-i^pq)  is  destabilizing  in  the 
sense  that  it  tries  to  drive  r negative,  and  this  has  been  observed,  by  the 
mechanism  of  coupling,  to  induce  rapid  changes  (i.e.,  jumps)  in  the  values 
of  the  other  state  variables.  The  factors  most  responsible  for  causing  the 
jump,  then,  are  most  ccxnpletely  present  in  the  term  (-i^Pq) . The  positive 
eleven  deflection  (+2®)  causes  pitch  down,  q negative?  positive  aileron 
similarly  generates  negative  roll  rate  p;  with  i^  > 0,  the  coupling  term, 
then,  is  also  negative.  Also,  aircraft  which  are  most  susceptible  to  jump 
phenomena,  such  as  aircraft  B,  are  members  of  the  class  of  so-called  high- 
performance  aircraft.  Their  gecHnetry  is  characterized  by  thin,  stubby  wings. 
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which  means  low  aspect  ratio,  and,  critically,  low  The  relatively 

lower  values  of  mcikes  i^^,  i^  find  i^  larger,  thereby  enhancing  greatly 
the  (destabilizing)  effect  of  the  yaw  rate  inertial  coupling  term.  Another 
effect  is  that  the  roll  damping  term,  which  involves  is  related 

directly  to  wing  size.  As  this  and  similar  terms  grow  smaller,  the 
system's  ability  to  dissipate  energy  diminishes. 

The  sequence  in  which  the  magnitudes  (and  sign)  of  the  important 
force  and  moment  terms  change  is  also  significant.  The  principal  agent 
for  this  is,  again,  coupling,  as  will  be  seen  following  a detailed  descrip- 
tion of  the  history  of  each  of  the  moment  terms  which  make  up  the  roll  (p) 
and  yaw  (f)  equations. 

Five  terms  comprise  the  roll  moment  equation:  command  (IL^JSa) , damping 

(Jl  p) , dihedral  , cross  derivative  (A  r) , and  inertia  coupling  (-i,qr). 

p P r 1 

Of  these  terms,  for  the  flight  condition  under  discussion,  the  damping, 
inertia  and  cross  derivative  terras  are  seen  as  "stabilizing";  that  is, 
balancing  the  moments  due  to  the  command  and  the  dihedral  term.  In  this 
condition,  the  dihedral  effect  is  destabilizing  - the  rolling  moment 
generated  by  the  command  deflection,  6a,  is  enhanced  by  the  dihedral  term. 

In  the  context  here,  "destabilizing"  terms  - after  transients  have  decayed, 
at  t - 2 seconds  - are  negative,  and  the  "stabilizing"  ones  are  positive. 
Refer  again  to  Figs,  3,24  to  3.27. 

For  the  6a  = 3.7“  case  (stable),  after  transient  decay  at  about 
t = 2.0  sec.,  the  stabilizing  terras,  which  contribute  the  greatest  to  the 
value  of  the  A-term,  and  Include  the  I-terra,  in  roll  (Fig,  3.24),  achieve 
a balance  with  the  destabilizing  terras  (the  dihedral  term  in  A and  the 
C-term.  In  the  yaw  eucis,  the  critical  I-term  is  destabilizing,  but,  for 
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the  steible  control  setting,  it  is  seen  to  begin  to  decay  (Fig,  3.26).  This 
decay  again  maintains  a balance  among  the  yaw  moment  terms,  so  that  r 
rapidly  damps  to  0.0. 

For  the  6a  = 3.8*  case  (xuistable) , sho%m  in  Figs.  3.25  (roll)  Md 
3.27  (yaw),  6a,  by  directly  affecting  p in  the  yaw  I-term,  causes  enough 
of  em  Increase  in  magnitude  of  this  (destabilizing)  term  so  as  to  disrupt 
the  balance.  Fig. 3. 27  shows  this  critical  I-term  failing  to  decay  at  all 
(until  after  15  sec,  after  the  jump  has  occurred).  As  a result  of  this, 
r goes  negative,  starting  just  after  the  I-term,  at  2dx>ut  9 sec.  This 
changes  the  sign  of  r,  which  is  a coupling  •* agent"  into  the  roll  axis  (and 
also  the  pitch  axis) . The  most  direct  effect  in  the  roll  axis  as  a result 
of  this  jump  condition  is  seen  in  the  I-term,  Fig.  3.25.  The  effect  on  p, 
however,  is  quite  small  (recall  that  I for  roll  is  blown  up  in  the  figures). 
Figure  3.23  shows  that  the  p history  reflects  this.  The  behavior  of  r,  as 
induced  by  the  yaw  I-term,  is  most  noticeable  in  the  yaw  and  pitch  axis, 
as  the  histories  of  q euid  r (and  a and  3)  in  Fig.  3,23  show.  Although 
the  roll  I-term  plays  no  major  role  in  inducing  the  juit\p,  it  nonetheless 
may  be  useful  as  am  indicator,  or  a detector,  of  an  incipient  jump  (A  more 
"natural"  indicator,  vehicle  kinetic  energy,  is  introduced  and  discussed 
below) . The  terms  mostly  affected  by  the  jump  are  seen  to  undergo  tram- 
sitions  to  new,  and  much  different,  steady  state  values  following  the  jun»p, 
as  Figs.  3.25  and  3.27  show. 

It  should  be  eii5)hasized  here  that  these  greatly  divergent  results  arise 
from  a difference  in  the  value  of  6a  of  only  0.1*,  One  value,  6a  * 3,7*, 
lies  barely  on  the  stable  side  of  the  bifurcetion  locus  in  the  control 
space,  while  the  value  6a  - 3.8*  lies  barely  outside  the  stable  region. 
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An  emalogy  to  this  situation  of  the  control  state  being  in  the  neighborhood 
of  the  bifurcation  line  is  illustrated  by  the  "Catastrophe  Machine"  of 
E.C.  Zeeman  (See  Reference  [4]  of  Chapter  II) . The  machine  consists  of  a 
cardboard  disk  pivoted  about  its  center,  with  two  rubber  bands  fixed  to  a 
point  on  the  disk  near  its  perimeter.  The  other  end  of  one  of  these  is 
fixed  to  the  mounting  board  emd  the  "machine"  is  operated  by  moving  the 
other  riibber  band,  whose  free  end  location  may  be  considered  as  a control 
point.  The  angle  formed  by  the  fixed  point,  the  pivot  of  the  disk  and  the 
common  point  of  attachment  to  the  disk,  is  a measure  of  the  behavior  of  the 
machine.  It  is  not  hard  to  visualize  that  many  movements  of  the  control 
point  cause  a smooth  rotation  of  the  disk.  However,  when  the  control  point 
is  at  a bifurcation  point,  the  disk  will  suddenly  swing  with  large  amplitude. 
Near  this  point,  on  the  stable  side  of  the  bifurcation  point,  there  is 
basically  a balance  between  "stabilizing"  emd  "destedsilizing"  moments 
produced  by  the  rubber  bands  attached  to  the  fixed  point  emd  to  the  control 
point,  respectively.  This  is  analogous  to  the  6a  = 3.7®  case  described 
2dx}ve.  At  this  point,  on  the  verge  of  a catastrophe,  only  a minor  distur- 
bance is  necessary  to  induce  a jump.  The  balance  between  the  two  roughly 
equal  and  opposite  moments,  near  the  bifurcation  point,  is  a precarious  one. 
The  value  of  6a  which  is  just  large  enough  to  force  r to  go  negative  (in 
this  particular  flight  regime) , as  dictated  by  the  yaw  I-term,  will  be 
that  value  which  causes  a jump. 

The  cdsove  discussion  is  not  intended  to  minimize  the  importance  of  the 
dihedral  term  on  roll  behavior.  Indeed,  a coiqmon  simplification  is  to  neglect 
totally  one  of  the  roll  rate  jump-indicator  terms,  (-i^^qr) . Runs  made  first 
with  ij^  » 0,  and  then  with  both  i^^  and  = 0,  both  show  a jump  occurring 


SCIENTIFIC  SYSTEMS,  INC. 


-82- 


near  6a  * 4*  (6  = 2®,  6r  * 0*),  seme  Initial  conditions  as  before.*  This 

e 

justifies  the  simplification.  Again;  the  principal  importeuice  of  the  r-terms 
with  regard  to  roll  behavior  lies  in  their  role  as  indicators  of  a jump 
condition. 

A study  was  also  done  on  the  effect  of  adding  the  gravity  terms  to  the 
system  equations.  A run  made  with  the  same  control  history  and  initial  con- 
ditions as  the  one  depicted  in  Fig.  (3.14),  except  for  the  addition  of  the 
gravity  terms,  shows  very  similar  behavior,  as  Fig.  (3.28)  indicates.  This 
result  agrees  with  observations  made  by  Schy  2md  Hannah  [9]. 

A question  now  arises  as  to  whether  or  not  there  may  exist  a more  general 
dynamic  Indicator  of  incipient  jump  conditions.  This  is  an  Importemt  concern. 
Insofar  as  the  time  history  of  a jixnp  situation  manifests  only  very  subtle 
differences,  prior  to  the  jump,  compared  to  a neighboring  non- jump  control 
setting.  If  the  pilot  had  a reliable,  real-time,  on-board  dynamic  indicator 
of  incipient  jump  conditions,  he  would  then  have  a much  better  chance  of 
taking  the  proper  action  to  return  to  stable  control  settings.  Providing 
such  information  to  eui  autopilot  obviously  has  similar  benefits.  The  studies 
done  with  aircraft  B seem  to  Indicate  that  monitoring  the  rate  of  change  of 
total  vehicle  kinematic  energy  might  very  well  play  an  important  role  in 
creating  a general  control  strategy  for  avoiding  jump  phenomena.  In  any 
event,  it  adds  useful  insight. 

Consider 

E * j m + Y I U • U (3.32) 

where  E is  the  vehicle  total  kinetic  energy,  m is  the  vehicle  mass,  V the 
velocity,  1 the  inertia  tensor,  and  w is  the  vehicle  wgular  velocity, 
coordlnatlzed  in  a body-fixed  frame  consistent  with  1. 

*A  run  with  only  ij^  ■ 0 still  produces  a jump  for  6e  ■ 2®,  6a  neeu:  4®, 


E » mV  (oa  + $6)  + 


(3. 


where  E is  the  time  rate  of  change  of  E and  the  other  terms  have  been 
identified  elsewhere. 

Eq.  (3.33)  was  confuted  for  the  6a  = 3.7“  zmd  3.8“  snins  with 
aircraft  B,  6e  = 2“,  and  the  preliminary  results  (see  Figs.  3.29  and 
3.30)  show  that,  for  stzdale  6a  (3.7“),  E is  positive  and  increases  during 
the  transient  period,  but  then  decreases  emd  changes  sign  shortly  there- 
after, then  decays  negatively  towards  zero.  For  the  jump  situation, 
however,  there  is  an  inflection  point  in  E as  it  decreases  towards 
zero,  and  it  again  starts  increasing.  This  inflection  point  occurs  in 
time  just  before  the  point  that  yaw  rate  goes  negative  (about  4,5  sec,). 
The  resultant  buildup  in  E because  of  E's  failure  to  go  negative  in  a 
jump  situation,  until  the  jump  is  completed,  is  interpreted  as  imparting 
the  vehicle  with  enough  energy  to  jump  the  energy  barrier  separating 
the  current  state  from  a new  equilibrium  state  of  higher  energy.  Once 
the  jump  has  occurred,  there  is  a higher  level  of  energy,  and  E does 
decay  to  the  origin.  Fig.  3.30  shows  this  for  the  6a  =»  3.8“  jump  setting 
The  jump  is  well  underway  at  about  t = 15  sec,  (see  also  Fig,  (3.23) 
for  a time  history  of  this  case) , where  E is  a local  maximum.  There- 
after, E falls  off  in  magnitude  towards  zero. 

The  connection  between  the  behavior  of  the  vehicle's  kinetic 
energy  and  jump  phenomena  is  not  mere  coincidence.  The  equilibrium 
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values  of  emy  of  the  state  variables,  for  aircraft  B,  may  be  interpreted 
21S  loc2d.  extrema  of  a potential,  or  Liapunov,  function  ^ . For 
exanple,  if  equations  (3.8-3.12)  are  solved  for  equilibrium  values  of  the 
state  variables,  the  algebraic  system  reduces  to  a scalar,  polynomial 
equation  of  ninth  order  for  p,  in  the  form  (see  Appendix  B) 


f (p,  6a,6e,6  r)  * 0 (3.34) 

P ® 

after  q,  1r,a,  euid  have  been  eliminated  using  equations  (3.8-3.11) 

[A  fifth  order  polynomial  in  p results  if  the  (-i^qr)  term  in  Eq.  (3.12) 
is  neglected] . From  this  equation,  the  Liapunov  or  potential  function 
is  defined  eis 


Thus,  ^ may  be  eeisily  generated  numerically  by  integrating  (~fp)  over  p, 
with  the  control  veuriables  fixed.  J(p)  will  typically  have  maxima  and 
minima.  The  extrema  of  T correspond  to  the  equilibrium  points,  emd  the 
e]q>res8ion  will  determine  whether  the  extremum  in  question  is  a 

«jp 

maximum  or  minimum.  Since 
that  the  minima  of  J (stable  points)  occur  where  < 0,  and  conversely 


this  is  equivalent  to  saying 


for  the  meucima  (unstable  points) . This  is  consistent  with  the  results 
of  linearization,  which  require  the  eigenvalues  of  the  Jacobian  matrix 
to  have  negative  real  parts  for  stable  behavior  in  the  neighborhood  of 
an  equilibrium  solution.  J changes  shape  as  the  controls  vary,  and  the 
point (s)  in  the  control  space  idiere  the  number  of  extrema  of J changes 


A potential  function  can  be  defined  for  the  case  in  which  the  equilibrium 
equations  can  be  reduced  to  a single  equation,  as  is  shown  in  Appendix  B. 
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(as  one  or  more  extrema  beccxne  inflection  points)  are  bifurcation 
points,  2Uid  their  locus  generates  the  bifurcation  surface.  A system 
state  ^Ich  is  at  a local  minimum  of  ^ which  disappears  in  this  way 
(as,  say,  6 a Increases)  then  "jtxmps"  to  a surviving  local  minimum. 

The  family  of  local  minima  of  in  the  state  space  are  called  "attractors". 
For  fixed  controls,  the  initial  conditions  determine  which  equilibrium 
state  is  attained  by  the  system  as  is  seen  in  Fig.  3.31  for  a representa- 
tive situation.  As  initial  conditions  vary  over  the  p-q  space  (controls 
fixed)  in  Fig.  3.31,  the  solution  may  or  may  not  diverge.  The  x's  are 
equilibrium  points,  [extrema  of  a Liapunov  (potential)  function^  (Pr<l)]> 

S meems  the  point  is  a stable  node,  U,  unstable  saddles. 

From  this  viewpoint,  J is  seen  to  generate  a potential  energy  field, 
with  "barriers"  and  "sinks".  Changing  the  controls,  say  da,  is  then  seen 
as  changing  the  potential  energy  field  of  the  system.  Referring  to 
Fig.  (3.31),  if  the  solution  is  in  the  vicinity  of  and  6a  changes  so 
that  moved  toward  Sj^,  then  a jump  occurs  for  that  6a  at  which  Uj^  merges 
with  S^,  cemcelling  the  extrema. 

Notice  that  the  sum  of  potential  and  kinetic  energies  (or  Hamiltonian) 
of  the  system  stays  constant.  When  the  equilibrium  state  is  unstable 
corresponding  to  a meucimum  of  the  potential  function,  emy  disturbance  will 
generate  an  exchange  of  kinetic  2tnd  potential  energies,  viz.,  em  Increase 
in  the  kinetic  energy  and  a decrease  in  the  potential  energy  of  the  system. 
Fig.  3.30  shows  this  increase  in  kinetic  energy  during  a jump  condition. 
After  the  jump,  the  exchange  between  kinetic  and  potential  energies  is 
completed  and  the  system  settles  to  a new  stable  equilibrium  point. 


i 


SCIENTIFIC  SYSTEMS,  INC. 


-86- 


f 


< I 


(3.6.2)  Aircraft  H. 

(3. 6. 2.1)  Equilibrium  Surface  Nonlinear  Root  locus. 

The  equilibrium  surface  for  aircraft  H is  a 3-dimensional  surface  in 
eight  dimensional  space  consisting  of  five  state  varied>les  (Pfq,r,a,8)  and 
three  control  V£uriables  (6a,6e,6r).  We  present  here  projections  of  the 
equilibrium  surface  on  different  plues  formed  by  pairs  of  state  emd 
control  veu:ied}les.  Following  the  approach  tedcen  for  aircraft  A emd  B, 
first  the  equilibrium  roll  rate  p versus  6a  curves  are  shown  for  different 
values  of  (6e,6r).  (Figs.  3.32-3.42).  a compeurison  with  similar  figures 
for  aircraft  B (Figs.  3.3  and  3.4)  reveals  that  for  aircraft  H,  up  to  9 
equilibrium  points  may  exist.  A linearized  stability  analysis  around  each 
equilibrium  point  reveals  that  unstable  equilibrium  points  may  have  one 
or  more  complex  eigenvalues  with  real  parts  in  addition  to  positive  real 
eigenvalues.  Points  denoted  by  letter  L have  one  pair  of  complex  eigen- 
values with  positive  real  part  whereas  points  denoted  by  E have  in  auldi- 
tion  a real  eigenvalue  in  the  right  half  pleme.  Interestingly  enough,  in 
most  of  the  cases,  only  one  equilibrium  point  is  str  /a  "sae  Figs.  3.32  to 
3.41)  for  a fixed  value  of  the  control  variable.  -v^ve'. , em  ex2unination 
of  the  root  locus  plot  (Fig.  3.47)  shows  stzd>ility  is  ccnerally  lost  by 
a pair  of  conplex  eigenvalues  crossing  the  imaginary  zucis.  Based  on 
the  discussions  in  Chapter  II,  this  indicates  possibility  of  Hopf  Bifurca- 
tions to  limit  cycles.  A conplete  representation  of  limit  cycles  requires 
a minimum  of  two  state  dimensions  chosen  to  span  the  eigenspace  of  the 
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^ L purely  imaginary  eigenvalues.  For  this  reason,  limit  cycles,  which 

> ■ 

^ I are  nothing  but  equilibrium  orbits  cannot  be  represented  on  figures  such 

V'  • 

i ^ 

} as  3.33.  In  the  next  four  sections,  we  would  show  limit  cycles  in  terms 
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of  time  histories  and  phase  plane  plots.  He  conplete  our  discussion  of 
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the  equilibrium  plots  here  by  referring  to  Figs.  3.43-3.46,  which 
show  the  remaining  state  variables  (q,r,a,3)  versus  6a  for  6e  » 6r  ■ 0® 
and  to  Figs.  3.59-3.63  which  show  equilibrium  states  with  respect  to  6e 
cuid  6r.  The  effect  of  neglecting  different  terms  in  Eqs. (3.8)- (3.12)  for 
aircraft  H on  the  equilibrium  surface  cam  be  seen  from  Figs.  3.64-3.66.  A 
cotiparison  of  Figs.  3.69,  3.65  and  the  corresponding  figure  with  (-i^qr) 
neglected  (not  shown  here)  reveals  that  there  is  virtually  no  difference 
between  these  three  cases.  However,  Fig.  3.66  shows  that  neglecting  all 
the  nonlinear  aerodynamic  terms  for  aircraft  H has  a serious  effect, 
particularly  on  the  minor  branches  of  the  equilibrium  surface.  It  is 
interesting  to  note  that  for  |5a|  < 10®  the  shaipe  of  the  p vs.  6a 
curve  is  quite  similar  in  all  the  cases. 

(3. 6. 2. 2)  Hopf  Bifurcation  to  Limit  Cycles 

Figures  3.48  amd  3.49  show  the  projection  of  limit  cycles  on  the 
p-aixis  for  different  values  of  6a  obtained  by  solutions  to  Eqs. (3.8)-(3.12) . 
A con^arison  of  the  two  figures  shows  the  hysteresis  phenomena,  which  is  a 
consequence  of  the  finite  domains  of  attraction  surrounding  limit  cycles. 

It  is  seen  frcrni  Fig.  3.48  that  when  6a  is  decreased  from  zero  to  negative 
values,  limit  cycles  appear  at  6a  -12®  when  the  equilibrium  point  chamges 
from  an  S to  Ir-type  (representing  the  crossing  of  a ccnnplex  pair  of  eigen- 
values into  the  right  half  plaine) . The  limit  cycle  first  grows  amd  then 
decays  in  amplitude,  finally  disappeauring  around  6a  -30®.  A stad}le 
limit  cycle  generally  disappears  by  collision  with  an  unstable  limit 
cycle.  Figs.  3.48  and  3.49  only  show  the  stable  limit  cycles  obtained 
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by  solving  the  equations  of  motion. 


SCIENTIFIC  SYSTEMS,  INC. 


-88- 


r 

i 


i 

1 


1 


I 


When  6a  is  Increased  from  negative  to  positive  values,  Figs.  3.49 
shows  that  limit  cycles  exist  for  6a  values  up  to  -2®.  In  fact  for 
6a  B -2®,  6e  « 2®,  the  state  portrait  consists  of  at  least  5 equilibrium 
points  (one  stable,  4 unsteible)  and  2 equilibrium  orbits  or  limit  cycles 
(1  st2U3le,  1 unstable) . The  stable  equilibrium  point  Is  surrounded  by 
em  unstcdile  limit  cycle  which  in  trun  is  surrounded  by  a sted}le  limit 
cycle.  Therefore,  depending  on  the  Initial  state,  the  aircraft  trajec- 
tories may  either  converge  to  the  st2d3le  equilibrium  point  (Fig.  3.48) 
or  to  the  limit  cycle  (Fig.  3.49).  Actual  phase  plane  plots  showing 
this  phenomena  will  be  presented  in  section  3. 6. 2. 5.  The  existence  of  two 
steQsle  limit  cycles  for  a particular  value  of  6a  will  also  be  demonstrated. 
To  con^lete  our  discussion  of  Hc^f  Bifurcation,  we  show  in  Fig.  3.30  time 
history  plots  where  6a  is  Increased  in  steps  from  -10®  to  -20®.  The 
development  amd  growth  of  limit  cycles  is  seen  clearly  in  this  figure. 

(3. 6. 2. 3)  Bifurcation  Surface 

Figs.  3.51,  3.52  and  3.53  show  the  bifurcation  surface  in  the 
control  space.  Compared  to  the  Butterfly  bifurcation  surface  for 
aircraft  B (Fig.  3.10),  the  bifurcation  surface  for  aircraft  H is  much 
more  complicated.  A three  or  higher  dimensional  visualization  of  the 
equilibrium  surface  c£m  be  obtained  by  combining  Figs.  3.51-3.53  with 
Figs.  3.32-3.42.  The  presence  of  eui  umbilic  (possibly  parediollc  type 
[14])  is  indicated  by  the  shape  of  curves  around  N emd  N'.  These  regions 
correspond  to  the  case  where  two  distinct  eigenvalues  of  the  linearized 
system  go  to  zero  simultaneously.  The  dotted  bremches  (1,1',  J,J',  K,K', 

L and  L')  on  Fig.  3. 51a  represent  (6a,  6e)  points  at  which  a pair  of 
conplex  eigenvalues  crosses  the  imaginary  axis,  resulting  in  Hopf 
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Bifurcations  to  limit  cycles. 

The  bifurcation  breuiches  (solid  lines)  separate  each  pleme  into 
a group  of  individual  regions.  Each  point  (corresponding  to  a specific 
control  setting)  in  euiy  specific  region  has  the  same  number  of  equilibrium 
solutions  as  any  other  point  in  the  same  region,  and  the  number  of  such 
solutions  is  indicated  by  the  number  which  appears  in  the  region.  Catas- 
trophic behavior  can  occur  when  the  controls  are  varied  so  as  to  cross 
one  of  the  bifurcation  branches  to  a region  with  a fewer  numlser  of 
equilibrium  solutions,  as  some  of  the  steible  solutions  disappear. 

The  bifurcation  behavior  arotuid  M eind  M'  (see  also  Fig.  3.42)  in  Fig.  3.51a 
requires  further  study.  A proper  representation  of  this  laehavior  would 
require  the  use  of  two  behavior  or  state  variables  corresponding  to  the 
case  of  an  umbilic  catastrophe.  Further  con4>utations  are  also  needed 
to  obtain  the  limit  cycle  bifurcation  boundaries  in  the  (6e,6r)  and 
(6a, 6r)  planes.  Both  the  adjove  cementations,  however,  are  straightforward 
and  are  not  expected  to  present  any  additional  conelications . 
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(3. 6. 2. 4)  Qualitative  Dynamics  for  Aircraft  H 

Dynamically,  the  added  complexity  of  the  aircraft  H model  vis-a-vis 
aircraft  A amd  B tends  to  produce  more  periodic  motions  and  even  tighter 
coupling  them  the  simpler  models  which  describe  these  aircraft.  As  a 
consequence  of  this,  one  expects  to  observe  some  of  the  phenomena  which  are 
oscillatory  in  nature,  and  which  have  not  been  observed  in  aircraft  A 
and  B.  These  would  Include  limit  cycle  phenomena,  buffeting  and  wing 
roclc,  in  addition  to  the  phenomena  which  have  characterized  the  high-a 
behavior  of  the  aircraft  discussed  in  Sec.  3.6.1. 
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The  basic  numeuver  to  be  investigated  is  similar  to  the  meuieuver 
applied  to  aircraft  B;  i.e.,  a pitch  dovm  (6e  = 2**),  with  a commanded 
aileron  deflection  near  the  stability  surface.  For  this  condition,  the 
equilibrium  surface  plot  of  (steady  state)  roll  rate  p vs  6a  (Fig.  [3.33]) 
shows  the  usual  linear  behavior  in  the  small  control  deflection  region 
which  is  one  of  five  possible  solutions  for  small  6a.  While  this  plot  has 
most  of  the  traits  of  the  corresponding  one  for  aircraft  B (Fig . [3.4]), 
there  are  extra  curves  for  high  |6a{  which  reflect  the  increased  order 
of  the  equilibrium  equation  for  p.  The  6e  = 0*  plot  shows  this  even  more 
convincingly,  as  these  "extra  curves"  have  moved  to  the  small  |6al  region, 
providing  nine  possible  solutions  for  p (Fig.  [3.32]). 

Returning  to  Fig.  [3.33] » a linearized  stability  2unalysis  at  selected 
points  along  the  equilibrium  reveals  whether  the  equilibrium  point  has 
stable  (S) , unstable  because  of  real  positive  roots  (U) , or  unstable 
because  of  right-half  complex  roots  (L) , eigenvalues.  The  L surfaces  were 
not  present  in  eeurlier  studies  (Sec.  3.6.1),  emd  they  indicate  regions 
of  possible  limit  cycles.  Refer  also  to  the  root  locus  plot  of  the 
eigenvalues  with  p as  the  parameter,  given  in  Sec.  3.6. 2. 2 (Fig.  [3.47]). 
In  Fig.  [3.33],  the  main  equilibrium  bremch  starting  from  (p  = 0, 

6a  s 0)  becomes  unstable  around  |6a{  « 14",  when  a pair  of  complex 
roots  crosses  the  imaginary  eocis  (see  also  Fig.  [3.47]).  If  the  aircraft 
is  excited  at  6a  = -12",  there  are  five  possible  equilibrium  points,  but 
!!ll  of  them  are  unsted>le.  The  Hopf  Bifurcation  Theorem  predicts  the 
possible  existence  of  a family  of  stable  limit  cycles.  Results  which 
confirm  these  predictions  will  be  presented  in  following  sections. 

If  5a  starts  from  0"  emd  increases,  the  solution  will  begin  to 
approach  the  limits  of  the  linear  region  of  steady  state  solutions.  Just 
prior  to  the  juiqp,  at  stable  6a  % ^ 12",  p % + 120"/sec.  (Fig.  [3.33]) » 
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there  is  observed  a small  region  where  stedile  limit  cycle  behavior  is 
experienced.  At  this  point,  one  expects  oscillations  in  all  axes  to 
develop,  somewhat  like  pre-stall  buffeting.  As  6a  increases  to  a jump 
setting,  the  equilibrium  states  jun^  to  limit  cycle  solutions.  Here, 
if  |6a|  were  decreased,  there  would  again  result  a hysteresis  effect, 
culminating  in  ein  autorotational  limit  cycle  whe  i |6a|  decreases  to  0.0®. 

On  the  other  hcind,  a further  increase  in  |6a|  would  result  in  the  equilibrium 
solutions  entering  a region  characterized  by  stable  equilibrium  points. 

(3. 6. 2. 5)  Numerical  Verification 

The  verification  of  the  limit  cycles,  predicted  in  the  previous 
section  as  a result  of  analysis  of  the  equilibrium  surface  plots  (e.g., 

Fig.  [3.33]),  is  done  first  by  solving  the  equations  of  motion 
nvimerically  for  different  values  of  6a  from  -12®  to  -40® , and  6e  = 2® . 

The  results  shown  in  Fig.  [3.50]  are  again  in  general  agreement  with 
the  equilibrium  plot  for  this  case  (Fig.  [3.32]).  The  significance  of 
these  results  is  that  the  family  of  limit  cycle  oscillations  shown  are 
believed  to  be  the  first  case  of  five-dimensional  limit  cycles  shown 
to  exist  for  an  aircraft.  The  appearemce,  growth,  decay,  and  extinction 
of  these  limit  cycles  is  typical  of  Hopf  Bifurcations;  and  this  has 
also  been  conjectured  to  play  an  important  role  in  the  onset  of  tur- 
bulence (see  Ruelle  and  Tcikens  [12]. 

As  in  the  case  of  aircraft  B (Sec.  3.6.1),  time  history  runs  were 
made  to  study  the  behavior  of  aircraft  H in  the  vicinity  of  its  "inner" 
jump  control  setting  for  6e  = 2®  viz.  6a  12®.  Refer  to  the  roll  rate 
equilibrium  plot  (Fig.  [3.33])  for  a global  perspective  of  behavior  in 
this  region.  Time  history  plots  for  aileron  values  equal  to  8,10,11,12, 

14  and  16®  are  shown,  respectively,  in  Figs.  [3.67  - 3.72]).  For  6a  = 8,10 
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and  11°,  the  time  response  is  linear,  and  the  steady  state  values  for  each 
value  of  £a  behave  as  indicated  by  the  corresponding  equlllbrlun  surface 
plots.  There  is,  however,  a noticeable  tendency  towards  increasing  oscilla- 
tory response  as  5a  increases  in  this  linear  region.  At  5a  > 12°,  just 
beyond  the  bifurcation  value , the  jump  is  quite  obvious  at  t >>  3 sec . , 
and  the  post  jump  trajectory  grows  to  a limit  cycle.  This  behavior  is 
quite  accurately  predicted  by  Fig.  [3.33],  since  Hopf  Bifurcations  to 
limit  cycles  occur  when  a pair  of  complex  eigenvalues  crosses  the  imaginary 
axis  from  the  left  half  plane.  The  final  two  figures  in  the  sequence 
(Fig.  [3.71,3.72])  show  that  the  limit  cycle  amplitude  decreases  and  finally 
vanishes  with  increasing  5a.  This  behavior  is  also  in  agreement  with  Fig. 
[3.33],  as  there  one  notices  the  stzdsle  (S)  region  of  the  upper  equilibrium 
bremch  which  exists  for  high  values  of  |5a|.  Except  where  noted,  gravity 
is  neglected  in  these  and  subsequent  runs. 

The  hysteresis  phencxnenon  was  adequately  discussed  and  verified 
for  the  simpler  models  of  Sec.  3.6.1,  and  thus  will  not  be  repeated  here. 
There  is  a clear  hysteresis  situation  for  aircraft  H as  well,  setting 
up  an  autorotational  state.  Figs.  [3.73]  emd  [3.74]  show  the  hysteresis 
behavior  of  aircraft  H,  with  multiple  solutions  for  a = 11.2°  zmd 
5e  = 2°.  The  initial  conditions  for  Fig.  [3.73]  are  p » -80  deg/sec, 
q = -8.1  deg/sec,  r = 2.5  deg/sec,  O*  -2.5°,  $ » 3.80°;  and  for 
Fig.  [3.74],  p « -310  deg/sec,  q *•  -23.69  deg/sec,  r “ -38.56  deg/sec, 
a = 6.92°,  3 > 6.14°.  For  both  cases  5a  was  set  to  11.2  °at  t » 2 sec., 
from  an  initial  setting  of  6«  (Fig.  [3.73])  and  13.5  (Fig.  I3.74J). 

This  is  not  an  autorotational  situation,  but  the  figures  do  show  great 
differences  in  behavior.  There  is  nearly  limit  cycle  behavior  in  Fig. 
I3.73],  with  much  less  severe  amplitude  and  frequency  than  is  seen  in 
Fig.  [3.74],  where  limit  cycle  behavior  is  very  much  evident. 


The  distinguishing  characteristic  of  the  aircraft  H model,  with  its  non- 
linear aerodynamic  coefficients,  is  the  existence  of  limit  cycles  in  the 
motion,  over  large  portions  of  the  flight  regime.  BACTM  provides  a sys- 
tematic approach  for  predicting  the  existence  of  limit  cycle  behavior, 
exhibited  by  time  histories  such  as  shown  in  Figs.  [3.50]  and  [3.74]. 

As  fvirther  validation,  a "finer  grid"  study  of  the  response  of  air- 
craft H was  conducted  for  6e  = 2®  eu>d  6a  in  the  neighborhood  of  the  bifurca- 
tion point  6a  = 12®.  A family  of  plots  for  6a  = 11.4,  11.46,  11.48,  11.5 
and  11.6®  is  presented,  respectively,  in  Figs.  [3.75]  to  [3.79].  The 
initial  conditions  for  these  runs  are  p(0)  = -80.  deg. , q(0)  = -8.1{deg/sec. , 
r(0)  = 2.5  deg/sec,  a(0)  = -2.5®,  3(0)  = 3.8®.  From  these  plots,  it  can 
be  deduced  that  the  critical  value  of  6a  at  which  danping  of  an  oscillatory 
mode  is  zero  lies  between  11.46  and  11.48®.  For  values  of  6a  just  short  of 
the  critical  value  (Figs.  13.69,  3.75  3.76]),  there  is  excitation  of  steady 
amplitude  and  frequency  oscillations  of  very  low  dainping.  The  case  6a  = 
11.46  (Fig.  [3.76])  was  run  for  64  sec.,  much  longer  them  the  aircraft 
characteristic  times,  and  there  is  only  a mild  decay  in  anplitude.  There- 
fore this  value  for  6a  is  the  critical  value,  for  all  practical  purposes. 

Once  the  junp  has  occurred,  for  6a  * 11.48®  and  11.50®  cases,  there 
is  no  question  that  growth  to  a limit  cycle  is  the  predominant  phenomenon. 
Again,  the  state  space  plots  of  Sec.  3. 6. 2. 6 (Figs.  [3.87]  and  [3.88]) 
amplify  this  result.  Observe  also  that  the  limit  cycle  magnitude  is  greater 
for  the  smaller  value  of  6a.  Fig.  [3.80]  also  shows  this  tendency,  as  it 
presents  a time  history  for  4 different  values  of  6a,  13.5,  12.,  11.4  and 
11.2®,  with  6e  ■ 2®  and  6r  * 0®,  as  usual  and  the  state  initial  conditions 
on  the  post-jump  equilibrium  branch.  Figs. [3.87]  and  [3.88],  plotting  p vs. 
the  other  state  variables,  runs  for  72  sec.,  by  which  time  the  limit  cycle 
amplitude  is  cdeaTly  established.  The  time  history  shown  in  Fig.  [3.80] 
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nms  for  t * 32  sec.,  stopping  shortly  before  this  steady  state  amplitude 
Is  reached. 


(3. 6. 2. 6)  Physical  Explanation 

The  distinguishing  feature  of  the  model  for  aircraft  H Is  that  the 
added  nonlinear  aerodjrnamlc  terms  enhance  the  overall  nonlinearity  and, 
specifically,  cause  limit  cycle  behavior  over  large  regions  of  the  control 
space.  Even  though,  physically,  limit  cycles  are  very  different  from 
equilibrium  points  compared  to  the  non-osclllatory  jump  phenomena  observed 
for  aircraft  B,  their  sudden  appearance  may  be  thought  of  as  an  oscillatory 
jump. 

As  before,  the  roll  and  yaw  moment  equation  will  be  examined  as  being 
typical  of  the  remaining  state  variables.  For  aircraft  H,  there  are  four  extra 
moment-generating  terms:  the  pitch  cross  derivative,  and  three 

terms  linear  In  B,  r and  6a,  with  the  coefficient  being  proportional  to 

a,  (l  ra),  oB),  and  (1  , 6a).  The  Inclusion  of  such  terms  obviously 

ra  Bn  aoa 

makes  the  aircraft  H model  more  realistic  In  hlgh-a  regions.  The  added 
terms  are  all  nonlinear,  and  they  tend  to  make  the  behavior  generally 
more  oscillatory.  However,  the  basic  conclusions  of  Sec.  3,6.1  remain 
valid  here  also.  As  before,  yaw  rate  Is  the  variable  most  sensitive  to 
the  mechanisms  which  create  a jump  situation;  l.e.,  1^,  which  Is  large 
when  Is  small  (as  It  Is  for  high  performance  aircraft)  and  p and  q, 
which  respond  directly  to  the  control  settings  6a  and  6e.  For  6e  fixed 
at  a pitch  down  setting  (2**),  and  for  Increasing  6a,  the  yaw  moment  terra 
(-IjPq) , will  eventually  force  the  secular  part  of  r to  remain  below  zero 
(see  Fig.  [3.77]).  Thus,  as  before,  an  Imbalance  between  the  aero- 
dynamic and  Inertia  moment  terms,  produced  by  the  large  values  of  6a 
and  Ij  In  the  latter,  causes  the  jump  to  occur.  Figs,  3.83  to  3.86  show  the 
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roll  and  yaw  moment  time  histories  for  a non- jump  (6a  = 10**)  and  jump 
(6a  = 11. 6**)  control  setting.  See  Section  (3. 6. 1.6)  for  definitions 
of  the  A,  C,  I and  N terms;  and  Figs.  3.68  and  3.79  for  the  state 
variable  time  histories.  The  time  history  initial  conditions  for  these  runs 
are  p = -80**/sec,  q = -8.1**/sec,  r = 2.5‘’/sec,  ot  = -2.5**,  3 = 3.8**. 

Fig.  3.79  shows  that,  for  6a  = 11.6**,  the  jump  is  underway  by  about 
t = 6.0  sec.  However,  the  roll  I-term's  behavior  (Fig,  3.84)  is  seen 
to  begin  to  diverge  from  the  stable  (6a  = 10**,  Fig.  3.83)  I-term  by 
5.0  sec.  As  explained  for  the  aircraft  B case,  this  term  is  not  involved 
in  producing  the  jump,  but  seems  to  indicate  well,  if  magnified  enough, 
the  presence  of  an  incipient  jump.  Note  also,  particularly  in  the 
stable  r plot.  Fig.  3.85,  that  the  nonlinear  (N)  terms  play  a notice- 
able role  in  the  motion.  Further,  note  that  the  yaw  I-term  is  greater 
in  magnitude  than  the  control  (C)  -term.  The  yaw  I-term  for  the  unstable 
6a,  11.6**,  Fig.  3.86,  begins  to  diverge  in  Isehavior  from  the  stable  I-term, 
Fig.  3.85,  at  about  5.0  sec.  As  was  observed  with  aircraft  B,  this 
term's  behavior  for  aircraft  H is  similarly  the  main  catalyst  for 
inducing  the  jump.  This  is  clear  in  Fig.  3.86  as  the  I-term's  growth 
proceeds  that  of  the  other  terms.  As  before,  this  causes  r to  assume 
consistently  negative  values,  and  the  coupling  effect  in  pitch  and  the 
other  variables  is  obvious  from  Fig.  3.79. 
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h study  was  made  to  compare  the  effect  of  gravity  terms  on  the 
dynamics  of  aircraft  H.  Time  history  plots  for  dynetmlcs  without,  and 
with,  gravity  terms  are  presented,  respectively,  in  Figs.  [3.81],  emd 
[3.82].  The  initial  conditions  for  both  runs  at  6a  « 2*  are  p = 0.40 
deg/sec.,  q = -2.80  deg/ sec. , r * -.16  deg/sec.,  o * -2.35*,  g = 0.18*. 

The  elevator  setting  is  6e  * 2*  throughout.  As  for  the  aircraft  B study, 
the  results  are  quite  comparable,  the  most  noticeable  difference  here 
being  the  smoothing  of  peedcs  near  the  juBqp  region  %ihen  gravity  is 
included.  The  jump  consequently  does  not  occur  until  6a  1^,  slightly 
higher  than  the  critical  value  of  6a  when  gravity  is  neglected, 

A limit  cycle  is  a self-sustained  oscillation  in  a dyn2uaic 
system.  The  amplitude  of  this  oscillation  is  not  dependent  on  the  initial 
conditions  of  the  ;vroblem,  as  is  the  case  for  linear  systems,  '"it  it 
depends  on  the  nonlinear  nature  of  the  system  model  itself.  See  Figs. 
[3.74]  and  p.77]  for  exan^les  of  limit  cycles. 

Hopf  has  shown  that  a llsilt  cycle  always  exists  under  certain 
mild  conditions  when  a pair  of  roots  of  the  linearized  system  is  purely 
imaginary.  The  Liapunov  Center  Theorem  is  one  consequence  of  Hopf's 
results,  and  the  Liapunov  stability  analysis  of  the  equilibrium  solutions 
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for  aircraft  H have  shown  that  the  conditions  are  present  for  bifurcations 
to  limit  cycles.  Time  histories  of  several  flight  conditions  (see, 
for  example.  Figs.  [3.711,  [3.78],  [3.81],  and  [3.82’)  have 
verified  this  behavior. 

Phase  plame  plots  of  the  state  vzu:i2d>le  p (roll  rate)  vs.  the 
remaining  vauriables  (q,r,a,6)  demonstrate  the  evolution  of  limit  cycles. 

An  entire  set  of  such  plots  is  presented  in  Figs.  [3.87]  and  [3.88]. 

These  plots  represent  a rtm  made  with  aircraft  H with  6e  ~ 2" , 5r  = 0" , 

6a  ■ 11.46®,  and  initial  conditions  p = -310  deg/sec.,  q = -23.7  deg/ 
sec.,  r * -38.6  deg/sec.,  a = 6.92®,  emd  3 = 6.14®.  Due  to  scaling 
within  the  plot  algorithm,  portions  of  the  large  zuiplltude  limit  cycles 
have  been  cut.  The  initial  conditions  place  the  equilibrium  state  on 
the  "outer"  solution  branch  of  Fig.  [3.33].  A very  similar  time  history, 
also  on  this  equilibrium  brzmch,  but  with  time  as  the  independent  variable 
is  given  in  Sec.  3. 6. 2. 5,  Fig.  [3.74  ].  The  major  differences  are  that 
in  Fig.  [3.74  ],  the  control  settings  are  6a  = 13.5®,  emd  then  11.2® 
for  30  sec.,  whereas  in  Figs.  [3.87]  and  [3.88],  6a  » 11.46®  is  used  for 
72  seconds  to  allow  the  limit  cycle  amplitude  to  reach  a steady  state 
value . 

If  the  above  run  is  repeated,  except  that  the  initial  conditions 
are  on  the  "inner"  equlibrium  solution  branch  - i.e.,  p > -80  deg/sec., 
q ■ -8.1  deg/sec.,  r ■ 2.5  deg/sec.,  a “ -.25®  emd  3 “ 3.8®,  the  solu- 
tion is  "attracted"  to  a completely  different  region,  as  Fig.  [3.89] 

(p  vs.  q)  shows,  in  caii%>arison  to  Fig.  [3.87].  Figure  [3.76]  is  the 
compemion  plot  to  Fig.  [3.89],  showing  time  as  the  independent  variable. 

A more  con^lete  series  of  time  history  runs,  varying  the  Intial  conditions. 
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would  define  the  barriers  between  attractor  regions  in,  for  exanple, 
the  p-q  plane.  However,  the  equilibrium  surfaces  themselves  perform 

this  function  far  more  efficiently. 

Fig.  (3.901  shows  how  a change  in  elevator  setting  to  a value 
outside  of  the  "Butterfly"  region  of  the  bifurcation  surface  affects 
the  subsequent  motion.  The  run  is  the  same  one  as  that  which  generated 
Fig.  [3.88],  except  that  the  initial  conditions  were  right  on  the 
Fig.  (3.  8]  limit  cycle,  and,  at  t - 8 sec.,  6e  was  changed  to  -2“. 

As  can  be  seen,  there  is  very  erratic  behavior  (plots  of  the  other  state 
variables  are  similar,  due  to  co^pling)  in  the  course  of  transition 
from  the  limit  cycle  to  a stcUsle,  steady  state  solution.  The 
equilibrium  surface  predicts  the  results , but  not  the  transient  behavior 

very  much  in  evidence  in  Fig . (3 . 90l • 

Finally,  runs  were  made  from  the  same  initial  conditions  as 
the  runs  started  from  the  "inner"  equilibrium  branch  6^-2*,  r = 0®; 
however,  «a  was  set  to  10.6",  a non-junp  value,  for  one  run,  and  then 
to  13.5",  which  forces  a jui^p.  Plots  of  p vs.  3 for  both  cases  are  shown 
in  Fig.  [3.91],  which  shows  the  effect  of  the  "disappearance"  of  an 
"attractor"  equilibrium  solution,  as  it  merges  with  an  unstable  equilibrium 
solution  for  clumglng  6a. 

As  with  aircraft  B,  the  behavior  of  kinetic  energy  of  aircraft  H under 
jump  and  non- jump  conditions  was  studied,  and  Figs.  3.92  and  3.93  show  E and  E 
vs.  time  for  6a  » 10"  (no  jump  ) and  6a  = 11.6"  (jump)  settings  of  a 
rvin  with  6e  - 2",  6r  - 0",  and  zero  initial  conditions.  The  results 
indicate  large  exhanges  of  kinetic  and  potential  energy  typical  of  limit 
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cycle  motions.  Note  the  smoothness  of  the  (decaying)  oscillations  of  E 
in  the  stedjle  case.  Fig. [3.92],  as  opposed  to  the  behavior  of  E in  Fig. [3.93]. 
Table  6 and  Fig.  [3.79]  also  show  non-sinusoidal  inflections  in  E at  6.5 
and  7.5  seconds,  hinting  of  a jump,  which  does  occur  at  ^dx>ut  8 seconds. 
Unlike  aircraft  B,  E,  the  energy,  does  not  asymptotically  approach  a new, 
higher  value.  This  is  due  to  the  limit  cycle  nature  of  the  jump  state.  The 
time  history  of  the  6a  = 10*  run  is  shown  in  Fig.  [3.68]. 


(3.7)  Control  and  Stability  Augmentation  (A/C  B) 


The  effect  of  feedback  controls  on  the  bifurcation  behavior  may 
be  studied  by  using  Eqs.  (3.26)-(3.31) . It  is  seen  from  Fig.  2.9  that 


l2u:ge  positive  values  of  c^  produce  unique  stable  equilibrium  solutions. 
Consider  a negative  feedback  on  6a  from  roll  rate  viz.  6a  * “Kj^P  + 
where  > 0.  This  has  the  effect  of  changing  all  the  f's,  but  for 
simplicity,  let  us  concentrate  on  the  change  in  f^  and  f'.  The  new 
values,  denoted  by  f£  and  f^  would  be: 


f'  = 326. 7Kj^  - 21.64 


f'  = 11752.78  6e  + 23015Kj^  - 1525.9 


(3.34) 


(3.35) 


Thus,  negative  feedback  (K^  > 0)  increases  the  values  of  f^  euid 
and  by  a proper  choice  of  both  cm  be  made  positive.  The  effect 
on  c^  is  seen  from  the  last  term  in  Eq.  (3.16).  It  is  seen  that  for 
* 0 and  6e  > .13  rad,  this  term  is  negative  emd  increases  with 
6e.  The  net-effect  is  a decrease  in  c^  which  eventually  results  in 
unstedjle  behavior  for  (6a  » 0,  p = 0) . With  > 0,  this  problem  cm 
be  avoided  and  c^  can  be  increased  to  a value  for  which  the  stzdsle  region 
around  (6a  ■ 0,  p * 0)  is  conslder2d}ly  enlarged.  The  effect  of  more 
complicated  feedback  laws  can  also  be  studied  in  similar  fashion. 
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Control  emd  Stability  Augmentation  (A/C  H). 

Bifurcation  emd  equlllbrltn  surfaces  can  be  of  great  help  in 
deriving  stability  augmentation  systems  or  feedback  control  laws, 

«dd.ch  need  not  necessarily  be  linear.  An  example  involving  the 
aircraft  H model  will  now  be  considered.  More  detailed  analysis  will 
be  presented  in  forthcoming  reports. 

A combat-type  maneuver  of  basic  interest  is  the  rolling  pullout. 

This  meuieuver  typically  calls  for  coordinated  deflections  of  all  three 
controls  (Sa, fie, fir).  Lateral  response  is  often  enhanced  by  the  addition 

of  the  alleron-rudder-lnterconnect  (ARI)  and  stability-axis  yaw  danger 
command  augmentation  systems  [11,15].  The  ARI  system  causes  the  rudder 
(fir)  to  deflect  in  conjtinctlon  with  roll  control  inputs  (fia)  so  as  to 
eliminate  adverse  yaw  due  to  these  surfaces,  and  to  improve  roll  response 
in  general.  Typically,  the  ARI  gain(8)  are  scheduled  with  angle-of-attack 
In  implementation,  more  and  more  rxidder  deflection  results  from  a lateral 
control  stick  movement,  and  less  and  less  aileron  deflection,  as  angle- 
of-attack  Increases.  The  stability  axis  yaw  damper  deflects  the  rudder 
(fir)  in  response  to  a signal  proportional  to  (r-pa) . The  purpose  of  this 
feedback  is  to  reduce  sideslip  excursions  during  rolling  numeuvers  at 
high  angles-of-attack,  and  to  improve  lateral-directional  damping  in 
generiUl . 

Equivalent  lateral  stability  augmentors  can  be  derived  from  the 
bifurcation  surfaces  for  aircraft  H.  The  example  given  here  shows  results 
relatively  small  deflections  of  fie.  The  ARI  law  may  be  tiritten  as 


fir  ■ k(fie)fia 


(3.36) 
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In  Eq.  (3.36),  while  6r  is  linear  in  5a,  the  gain  k is  actually 
a parameter  veurying  with  5e  setting,  as  will  be  seen  below.  This 
law  can  be  rather  straightforwardly  extended  to  a more  general  form 
of  the  type 

6r  = f(5a,5e)  (3.37) 

Details  of  this  extension  will  appear  in  future  reports. 

The  form  of  Eq.  (3.36)  may  be  ascertained  directly  frcm  the 
bifurcation  surface,  of  which  Fig.  [3.53]  is  a typical  example. 

Fig.  [3.53]  shows  the  bifurcation  surface  in  the  (5a, 5r)  control  plane, 
with  5e  * 0*.  In  this  figure,  bifurcations  in  the  fom  of  jumps  occur 
when  a (5a, 5r)  locus  crosses  one  of  the  branch  lines.  The  branches 
A and  B ^u:e  of  particular  interest  because  they  axe  associated  with 
lower  values  of  equilibrium  roll  rate,  p;  hence,  these  branches  would 
mark  the  first  jump  encountered  as  (5a, 5r)  move  away  from  the  lineeu: 
(trim)  region,  which  is  centered  at  the  origin. 

If  5r  is  fixed  during  a roll  maneuver  (5a  varies) , then  the 
(5a, 5r)  locus  in  Fig.  [3.53]  is  a horizontal  line.  If  5r  « 0”,  this 
line  intersects  breuiches  A and  B at  1 5a  | =11.5".  The  corresponding 
equilibrium  plot,  p vs.  5a  for  5e  * 5r  * 0*(Fig.  [3.32]),  confirms 
that  jumps  will  occur  at  these  settings,  when  the  problon  starts  near 
5a  » 0*.  Thus,  the  "non-j«np"  region  for  5e  = 0“  and  5r  * 0“  is  |6a|  < 
11.5*.  However,  this  region  can  obviously  be  expemded  by  relating  5a 
and  5r  so  that  their  locus  passes  through  the  points  N and  N*  in 
Fig.  [3.53].  The  simplest  locus  is  a line  of  form  (3.36),  with  k ■ 
-0.152.  Fig.  [3.94]  shows  the  ^ vs.  5a  plot  for  this  case.  Notice 
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that  the  trim  (basic)  stability  region  has  been  esipanded  to  |5a|  £17* 
frcn  + 11.5*  for  the  6r  * 0*  case.  Time  histories,  for  the  6r  = 0®  case  and 
the  hRI  case  eure  compared  in  Figs.  [3.95  and  3.96]  , respectively. 

It  is  obvious  that  ARI  aids  not  only  in  stability,  but  in  perfomumce  as 
well.  Note  that  ARI  causes  decoupling,  even  for  the  large  6a  values. 

In  order  to  excite  the  instsibillty  at  6a  “ 18.5*,  it  was  necesscury  to 
hold  that  setting  for  much  longer  than  the  8 sec . sho%m  in  Fig . [ 3 . 96al . 
Fig.  [3.96b]  shows  that  when  6a  ■ 18.5*  for  an  extended  time,  a very  sdld 
divergence  ensues.  Thus,  for  ARI  when  6e  * 0* , the  effective  range  of  6a 
is  increased  upwards  to  |6a|  -36*.  It  is  therefore  obvious,  from  the  time 
histories  sho%m,  that  stability  augmentation  control  laws  c^m  rather  easily 
be  derived  from  the  bifurcation  plots. 

If  k is  not  the  value  which  intersects  the  points  N and  N',  the 
results  are  seen  to  be  less  than  optimal.  Fig.  [3.97]  shows  the  6e  0* 
case  with  k * +0.152  (the  "optionxn"  value  of  k for  this  case  is  -0.152). 

The  plot  is  p vs.  6a.  Note  how  the  stable  region  for  6a  has  shrunk  to 
+ 8”.  (Other  equilibrium  branches  are  not  shown  for  the  purposes  of 


clarity) . 

When  eleven  is  nonzero,  the  use  of  the  equilibrium  plots  becesws 
more  iitportemt,  as  vrell  as  the  bifurcation  plots.  This  is  because  limit 
cycle  behavior  enters  the  major  branch.  Fig.  [3.98]  desKmstrates  this, 
showing  a pitch  down  case,  6e  - 5*.  Thus,  for  coBa)ined  longitudinal- 
lateral  maneuvers,  where  eleven  as  well  as  rudder-aileron  deflections 
become  prominent,  there  are  two  kinds  of  stability  criteria  that  must  be 
considered:  the  first,  seen  already  for  the  6e  ■ 0*  case,  is  the  jui^  in 


There  is  one  unstable  root  at  this  value  of  6a,  which  has  a time  constant 
of  14.6  sec. 


I 


SCIENTIFIC  SYSTEMS.  INC. 


-103- 


iJ 

J 


stability,  indicated  by  the  bifurcation  loci;  the  second  is  the  entry  into 
a lialt  cycle  region,  tdilch  is  seen  to  occur  for  magnitudes  of  5a  which 
axe  less  than  values  at  vdiich  jumps  occur.  In  Fig.  [3.98  ] , the  jumps 
occur  at  |5a|  - 35*,  tdille  limit  cycle  behavior  may  be  expected  to  begin 
at  around  + 20*.  This  figure,  for  the  p2u:ticul2ur  case  5e  5*,  was  plotted 
using  the  ARI  relationship  (3.36),  with  the  "optiminn"  k being  used  (k  = 

-0. 250) . As  for  the  5e  « 0*  case,  the  value  of  k is  determined  as  the  slope 
of  the  line  tdiich  conne<^ts  the  points  N anfl  N'  on  the  relevant  blfiircation 
surface.  Fig.  [3.99].  It  is  worthwhile  to  consider  the  "optimum"  gain  k 
as  the  slope  of  the  line  which  connects  two  points  \diich  are  on  loci 
marking  the  onset  of  limit  cycle  behavior.  The  two  points  would  be  selected 
so  as  to  maximize  the  stable  range  of  <Sa. 

In  Fig.  [3.99],  only  the  major  bifurcation  branches  eure  plotted. 

Note  that,  for  5e  <■  5*,  the  non-jus^  region  of  5a  is  generally  quite 
insensitive  to  variations  in  k (the  ARI  gain) . But  it  turns  out  for  this 
case  that  the  non-limit-cycle  range  of  5a  is  greatest  for  the  value  of 
k (-0.250  here)  which  connects  the  points  N and  N*.  More  work  is  being 
completed  in  this  area,  viz.,  the  effect  of  the  ARI  gain  k on  both  sted>llity 
ranges  for  5a:  limit  cycle  and  jump  (roll  departure).  The  results  will 
appear  in  forthcoming  reports. 

Tisw  history  results  for  the  pitch-down  roll  numeuver  case  are 
presented  in  Fig.  [LOOa] . It  can  be  seen  that  there  is  limit  cycle  behavior 
both  at  the  pre-jus^  aileron  setting  of  24*,  as  well  as  the  post-jim^  value 
of  40*.  The  jump  occurs  at  around  17  sec.,  just  after  5a  is  changed  to  40*. 
These  results  are  consistent  with  the  predictions  of  Fig.  [3.98]. 
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The  studies  conducted  thus  far  show  that  definite  Improvement  In 
performance  results  from  doing  control  system  design  using  the  bifurcation 
surfaces.  Even  for  the  most  simple  situation,  a control  law  of  the  form 
(3.36),  there  are  significant  Improvements.  For  the  above  pitch-down  case, 
the  ARI  p vs.  6a  plot.  Fig.  [3.98],  shows  that  limit  cycle  behavior  will 
be  avoided  If  |6a|  <18*.  If  there  Is  no  ARI  (6r  - 0*;  see  Fig.  [3.34]),  this 
region  Is  only  + 13*.  In  this  case  (and  also  for  the  6e  - -5*  case,  pre- 
sented below),  the  jump  stability  region  Is  not  particularly  sensitive  to 
k,  as  Fig.  [3.99]  Indicates.  However,  for  smaller  values  of  6e,  as  noted 
earlier,  jump  stability  Is  particularly  sensitive  to  k (for  smaller  values 
of  fie,  l.e.,  |6e|£4*,  limit  cycle  behavior  Is  not  predominant  on  the  major 
equilibrium  branches).  Time  histories  are  compared  In  Fig.  [3.100)for  the 
fie  “ 5*  case.  For  the  same  fia  history,  there  Is  a noticeable  improvement 
In  response  when  ARI  Is  active  (Fig.  [3.1003])»as  opposed  to  no  ARI 
(fir  - 0*  Fig.  [3.100b]).  Note  especially  the  reduced  amplitudes  for  all 
of  the  fia  settings,  and  the  lack  of  oscillation,  after  transients,  at 
5a  ■ 12*.  The  rolling  pullout  maneuver  Is  a combined  pltch-up  (fie  <0") 
and  roll,  executed  together.  For  the  case  fie  ■ -5®,  the  "optimum” 

ARI  gain  k - 0.118.  This  gain  Is  optimum  In  the  sense  that  the  rudder 
control  Is  specified  by  Eq.  3.36,  and  In  the  sense  that  k has  a value  which 
Is  the  slope  of  the  line  connecting  the  N-N'  points  In  the  relative  bifurca- 
tion surface. 

The  equilibrium  plot,  pT  vs.  fia,  for  aircraft  H In  a rolling  pullout 
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maneuver  Is  given  for  the  ARI  (k  0.118)  and  non-ARI  Ot  ■ 0"; fir  “0") 
cases  In  Figs.  [3.101  and  3.40],  respectively.  As  for  the  fie  - 5*  case. 
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a comparison  of  these  figures  cleeurly  shows  that  limit  cycles  are 
avoided  over  a much  l2u:ger  6a  range  (;f32®)  with  ARI  than  without  (+20® ). 
Additionally,  time  history  comparisons  of  the  ARI  and  non-ARI  systems 
(Fig.  [3. 102]), as  for  the  0 eUid  5®  settings  of  6e,  show  improved  response 
over  all  of  the  values  of  6a.  Note  paurticularly  that  a does  not  change 
sign  when  6a  moves  to  30®  from  12®,  amd  longitudinal  buffeting  (a  oscilla- 
tions) is  greatly  reduced  at  this  setting  of  6a.  Furthermore,  roll  rate 
p has  a lower  average  value  at  6a  = 30®  for  the  ARI  case,  as  can  be  veri- 
fied by  comparing  the  respective  histories  of  roll  angle,  (P.  Note  further, 
for  6a  = 30“,  that  the  frequencies  of  oscillation  au:e  reduced  with  ARI, 
and  that  there  is  a tendency  for  the  oscillations  to  dampen  out.  Fig. 

[3. 102b] shows  the  damping  out  of  the  variables  more  clearly,  as  6a  is 
held  to  30®  for  16  sec.  This  tendency  is  not  evident  in  the  non-ARI  case; 
Indeed,  6a  =»  30®  is  a setting  for  possible  limit  cycle  activity,  according 
to  Fig.  13.40].  The  time  history  Fig. [3.102c] at  6a  30®  certainly 
supports  the  results  of  Fig.  [3.40]. 

Plots  of  the  ARI  gain  k versus  equilibrium  angle-of -attack  a show 
nonlinear  behavior  in  the  a < 10®  region;  however,  the  values  of  k as 
derived  from  the  bifurcation  surfaces  2u:e  in  the  neighborhood  of  the  ARI 
gain  values  presented  in  [ll] . Results  of  a comparison  of  the  BACTM  ARI 
gain  with  others  [ll,15]  will  be  presented  in  a later  report.  At  this 
time,  it  seems  that  BACTM  is  fully  capable  of  deriving  satisfactory 
command  augmentation  schemes,  as  the  above  example  indicates. 


Equilibrium  Roll  Rate  p vs.  Aileron  Angle  6a  for 
Aircraft  A (6e  = 0®) 
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■ ■ 0*.  Equilibrium  pitch  rate  vs.  6a  - Deg/Sec. 
Curve  Denote  Equilibrium  Roll-Rate) . 
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Aircraft  B 


-500  deg/sec  < p < 500  deg/sec. 


■E 


BIFURCATION 

POINT 


Nonlinear  root  locus  plot  (Numbers  along  the  locus  denote 
equilibrium  roll-rate,  fie  - 2*.  Different  syiri>ols  refer 
to  different  eigenvalues.  For  complex  eigenvalues,  only 
one  l•eld>er  of  the  pair  Is  shown. 
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Flg.  3.10:  A/C  B,  Bifurcation  Surface  In 

(5a,  5e)  plane.  The  numbera  In 
each  region  Indicate  number  of 
equilibrium  points;  6 ■ 0. 
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FI6  3ill:  Efftct  of  Incrtoting  tltvotor  angit  6 9 from  A to B on  aircraft 
rol I -ratfj  autorotation  and  doparturt  occur  at  C. 
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Equilibrium  Roll-Rate  vs.  6a. 
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Aircraft  B,  Equilibrium 
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Fig.  3.31: 


Relation  between  stable, 
unst2U)le  equilibrium  points  ^Lnd 
attractor  regions  in  p-q  space. 
(Qiialitative  Representation) . 


roll  rate  (deg  x 10/ sec)  vs.  aileron  (deg) 

Boxed  letters  are  reference  for  bifurcation  plots 
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roll  rate  (deg  x 10/sec)  vs.  aileron  (deg). 
Boxed  letters  cure  reference  for  bifurcation 
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Fig.  3.49:  Equilibrium  Surface  and  Limit  Cycles,  A/C  H. 

Letters  indicate  eigenvalues  at  equilibrium: 

U = unstable;  S = stable;  L = unstable  cyclic. 
Vertical  bdrs  indicate  maximum  magnitude  of  roll 
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Bifurcation  locus 

in  (5a, 5e)  plane,  5r  = 2";  A/C  H; 

Limit  Cycle  branches  not  shown. 


to  letters  on  Figs.  3.26  and  3.36.  (Breinches  1,1',  J,J',K,K',  L and  L'  corresponding 
to  purely  imaginary  eigenvalues  are  not  shown) . 


SCIENTIFIC  SYSTEMS.  INC. 


-165- 


I 

f 

i 


OOCSSf 

66Z89i 


LLZ99i 

LLZmi 

^,9K09-{- 

^9J09+ 

99«05+ 


3 CT.  3 

0*3  3*0 

w a 


THETfltXlO^  ) PHHXIO 


SCIENTIFIC  SYSTEMS,  INC, 


SCIENTIFIC  SYSTEMS,  INC. 


-177- 


0 

0 

0 

a 

u 

0 

(j 

Q 

0 

u 

D. 

0 


0 

0 

0 

D 

0 

0 

fl 


ITO.OO  180.00  190.00  200.00  2*0.00  220.00 


terms  time  history  for  5a 


-JMj  M^llipjll,,  IJ  H III  mill-lM 


SCIENTIFIC  SYSTEMS,  INC. 


-194- 


00- 021 


00-06 


00' or  00- 0 

2*»*a39/0ia 


00* or- 


00- oe- 
loay 


00-021- 


li 

V 


t i' 


u 


SCIENTIFIC  system; 


u 

u 

0 


DED)(X10 


SCIENTIFIC  SYSTEMS,  INC. 


u 
u 

Fig.  3.69;  A/C  H fie  <=  2°,  fir  - 0°,  fia  - 11.46° 

U Phase  plot,  p vs  q,  for  Initial  conditions 

on  linear  equilibrium  bremch  (t^  «=  25  sec) 
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Fig.  3.90;  a/C  H 6r  - 0°,  6a  « 11.46° 

Phase  plot,  p vs  8,  initial  condi- 
tions on  outer  equilibrium  branch 
t^  = 64  sec;  6e  ■ 2 , then  -2  at 

t « 8 sec. 
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Fig.  3.99;  Bifurcation  locus  in  (fia.fir)  plane  (units,  deg.). 

A/C  H;  6e  * 5*.  Only  the  major  branches  are  shown. 
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Fig.  3.101t  Equilibrium  surface,  A/C  H; 

6e  - -5*;  ARI  gain  k « 0.118. 
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Fig.  3.102b:  Ti»«  history,  X/C  H; 

Part  (b):  ARX  gain  K - 0.118;  transient 
re^wnse  at  5a  > 30*. 
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IV.  Conclusions  and  Recommendations  for  Further  Research 

(4.1)  Conclusions 

Based  on  the  emalysis  of  nonlinear  models  for  aircraft  A,  B and 
H,  we  Ccm  conclude  that 

(i)  A large  number  of  unsted>le,  jump  cuid  limit  cycle  phenomena 
occurring  at  high  angles-of-attack , can  be  analyzed  in 
a unified  fashion,  by  using  BACTM  (Bifurcation  Analysis  and 
Catastrophe  Theory  Methodology) . 

(ii)  BACTM  provides  a global  representation  of  the  equilibrium 

and  bifurcation  surfaces  for  nonlinear  dynamic  systems.  The 
qualitative  dynamics  of  the  system  for  different  initial 
conditions,  controls  and  system  par6uneters  can  be  obtained 
easily  from  the  equilibrium  and  bifurcation  surfaces.  In 
particular,  the  control  emd  parameter  values  for  which  jumps 
and  limit  cycles  appear  as  well  as  domains  of  attraction  are 
obtained  directly  from  the  equilibrivun  auid  bifurcation  surfaces. 
BACTM  also  provides  insights  into  the  control  emd  identification 
problems  for  nonlinear  dynamic  systems. 

(iii)  Five  degree-of- freedom  nonlinear  aircraft  models  in  terms 

of  roll  rate  p,  yaw  rate  r,  pitch  rate  q,  emgle-of-attack  a 
and  sldellp  angle  3 with  three  control  surface  deflections 
(6a,5e,5r),  when  2malyzed  using  BACTM,  reveal  the  presence  of  a 
large  number  of  limit  cycle  auid  jump  phenomena  which  enconpass 
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roll  co^pllng,  directional  divergence,  wing  rock,  buffeting 


and  other  phenomena  known  to  occur  at  high  euigles-of-attack . 


(iv)  For  aircraft  A and  B represented  using  linear  aerodynamics 


effects  only,  the  five  degree-of-freedom  model  can  be  reduced 


to  a scaleur  fifth  order  polynomial  model  which  is  diffeomorphic 


to  a "Butterfly"  catastrophe.  The  stability  and  control  properties 


of  this  model  reveal  two  sets  of  jun^s,  hysteresis  effects  emd 


autorotation  for  high  values  of  elevator  deflection.  An  analytical 


criterion  is  developed  to  predict  the  critical  elevator  deflection 


or  amgle-of- attack  for  autorotation.  The  effect  of  feedback 


control  on  bifurcation  behavior  is  studied  and  it  is  shown  that 


a control  law  which  keeps  the  eigenvalues  of  the  linearized 


system  in  the  left  half  pleme  leads  to  bifurcation- free  behavior. 


The  results  predicted  by  BACTM  are  in  close  agreement  with  the 


results  obtained  by  numerical  integration  of  the  equations  of 


motion,  including  the  effects  of  gravity.  Physical  explematlons 


axe  given  for  the  bifurcations  and  jumps  emd  it  is  shown  that 


chwges  in  yaw  rate  are  an  indicator  of  6m  incis^lent  jump  condition. 


(v)  The  five  degree-of-freedom  model  for  aircraft  H involving 


several  nonlinear  aerodynamic  effects  exhibits  Hopf  Bifurcations 


to  limit  cycles.  Several  families  of  limit  cycles  are  shorn  to 


exist  and  depending  on  the  initial  state  and  the  control  settings. 


the  system  trajectories  converge  either  to  one  of  the  stable  limit 


cycles  or  to  a stable  equilibrium  point.  The  limit  cycles  them- 


selves increase  and  decrease  in  magnitude  and  may  correspond  to 


phenomena  such  as  wing  rock,  pre-stall  buffeting  and  divergent 


oscillations  towards  spin  conditions.  A physical  e^qplanatlon  in 


SCiENTiFiC  SYSTEMS,  INC. 


-217- 


0 


[ 

i A 

L 

( 


terms  of  the  exchange  of  kinetic  and  potential  energies 
is  developed  for  jump  and  limit  cycle  behavior.  The  extreme 
sensitivity  observed  at  high  angles -of -attack  to  minor  design 
changes  is  also  easily  demonstrated  in  terms  of  the  bifurcation 
and  equilibrium  surfaces. 

(vi)  Results  on  kileron-Rudder  Interconnect:  In  Sec.  3.7,  the 
use  of  EACm  to  synthesize  a conmand  augmentation  system 
which  schedules  the  rudder  as  a lineeur  function  of  the 
aileron,  the  gain  being  a nonlinear  function  of  elevon, 
has  been  demonstrated.  The  values  of  the  ARI  (aileron- 
rudder-interconnect)  gains  vexe  derived  directly  from  the 
plots  of  the  bifurcation  surfaces.  These  gains  were  seen 
to  be  nonlinear  functions  of  angle -of -attack,  when  computed 
under  the  criterion  of  maximizing  the  range  of  aileron  values 
which  avoid  jtnp  or  limit  cycle  behavior.  Time  history  com- 
parisons show  significant  improvement  in  dynamic  response  and 
performance  using  the  BhCTM  ARI  gain,  as  opposed  to  nulling 
the  rudder.  The  results  indicate  that  scheduling  rudder  as 
a general  nonlinear  function  of  both  aileron  and  elevon  would 
enhance  globally  the  aircraft  performance  at  high  emgles-of -attack. 

(4.2)  Recomnendations  for  Future  Research 

The  next  step  in  tlw  application  of  BACm  to  high  angle -of -attack 
phenomena  is  to  use  six  degree -of -freedom  equations  and  study  spin 
entry,  prevention  and  recovery. 

The  aircraft  departure  into  spin  followed  by  aircraft  loss 
in  several  cases  is  a very  serious  problem.  The  methodology  developed 
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under  the  current  project  provides  a suit^lble  tool  for  a con^lete 
analysis  of  this  phenomenon.  Certain  dynamical  behavior  for 

aircraft  B and  H has  already  been  observed  which  has  similarities  to 
departure  into  spin.  An  exact  verification  of  this  fact  would  require 
solving  the  spin  equations  of  motion  and  relating  the  two  sets  of 
equilibiriura  points.  A first  step  %«ould  be  to  use  the  spin  equations 
used  by  previous  investigators  [1,2]  and  solve  them  for  equilibrium 
surfaces  and  bifiurcations  using  BACTH.  By  considering  the  most  general 
set  of  equations,  it  would  be  possible  to  obtain  ginning  as  well  as 
nonspinning  equilibrium  states,  including  others  that  may  not  have 
been  identified  before.  This  will  lead  to  a study  of  the  post-stall 
departure  phnomenon  and  entry  into  spin.  Since  the  equilibritsa  and 
bifurcation  surfaces  will  be  obtained  for  different  control  values, 
techniques  for  spin  prevention  and  recovery  would  be  investigated  easily. 
It  tpould  also  be  possible  to  explore  the  use  of  feedback  laws, 
stability  augmentation  techniques  and  different  aircreift  configurations 
on  the  bifurcation  surfaces.  Notice  that  the  characterization  of 
bifurcation  surfaces  would  iiuJce  the  concepts  of  incipient  spin,  stability 
margins,  etc.  much  more  exact  compared  to  current  heuristic  methods  for 
defining  these  variables. 

Another  direction  for  future  research  would  be  to  consider 
applications  of  BACIM  to  other  dynamic  vehicles  such  as  ships,  sub- 
marines and  missiles.  Similar  types  of  junps  and  instabilities  are 
known  to  occur  in  these  vehicles.  The  use  of  BACTM  for  nonlinear  control 
design,  model  structure  determination  and  experimental  design  also 
needs  further  investigation. 
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APPENDIX  A:  MAIN  THEOREMS  OF  BACTM 


I 

I 

i 


A.l  CentT  Mzmifold  Theorem  or  Reduction  Theorem, 

[18,19,20,241 ;* 

Suppose  that  a set  of  n Ordinary  Differential  Equations  can  be  written  as: 

u « Au  + £,  (u,v,w)  (A.l) 

A 

V = Bv  + fg  (u,v,w) 
w = Cw  + f ^ (u,v,w) 

v^ere  A,B,C  are  square  matrices  of  orders  j,k,£  respectively 
( j + k + £ ==-n)  whose  eigenvalues  have  their  real  parts 

r 

<0,  = 0,  >0  resp.  Hie  nonlinear  functions  f , f , f are  C 

ABC 

differentizdile  (r^2)  euid  vemish  along  with  their  derivatives, 

r-1 

at  (u,v,w)  * 0.  Then  there  is  a c suboamifold  M,  of 
dimension  k,  which  passes  through  (u,v,w)  * 0»  is  tangent 

there  tou«0,  w=0,  euid  which  is  invariauit  under  the  adiove  , 

equations  (i.e.,  is  composed  of  solution  curves  of  those 
equations) . 

Proposition  [20]  i Suppose  there  is  a curve  (u(li),  v(M),  w(Vl))  depending  on  a | 
parameter  or  control  vector  p such  that  (u(0),  v(0),  w(0))  =0  and  for  each 
pe  R™ , [u(p),  v(p),  w(p)]  is  a fixed  point  for  Eq.  (A.l).  Then  for  all  ^ 

sufficiently  small  p,  [u(p),  v(p),  w(p)J  lies  in  any  center  manifold  M. 

A. 2 Main  Theorem  of  Element€ury  Castrophe  Theory  [1,2,3]; 

Let  xEr”  €uid  ceR™.  Then  f(x,c)  » 0 is  eui  m dimensional  1 

n+m  ^ 

nuuiifold  M in  R . Define  the  projection  map  X;  M C where  C 

is  the  space  of  control  variables. 

* Reference  numbers  same  as  for  Chapter  II  I 
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Let  F denote  the  space  of  c«  -functions  on  ^n 


with  the  Nhltney  c««>-t<^logy. 


Theorem:  If  nKS,  there  is  an  open  dense  set  F^F  which  is  the 
set  of  generic  functions.  If  4*  is  generic,  then 

1)  The  mzmifold  is  2m  m-manifold; 

2)  Any  singularity  of  the  projection  map  is  equivalent' to  one 
of  a finite  number  of  types  called  elementary  catastrophes; 

3)  Xo^  is  stable  under  small  perturbations  of  4>.  The  number  of 
elementary  catastrophes  depends  only  t^on  m,  the  dimension  of 
control  space,  as  follows: 


Elementary 

Catastrophes 


I 1 I 2 I 3 I 4 1 5 16 


7 11 


1 2 


Here  equivalence  Implies:  two  me^s  and 

2u:e  equivalent  if  there  exist  diffeonorphisms  h 2md  k such 
that  the  following  diagram  is  coamutative. 


If  and  X^  have  singularities  at  *g^g' 

respectively,  then  the  singularities  are  equivalent  if  the 

above  definition  holds  locally  with  hx^'^„.  Stable  aeams  that 

t 9 

X^  is  equivalent  to  x^  for  all  g in  a neighborhood  of  ^ in  F. 


11 


Ml 
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A.3  The  Hopf  Theorem  in  R t [S) 

The  following  theorem  is  essentially  due  to  Andronov 
(1930)  cuid  Hopf  (1942)  euid  was  suggested  in  the  work  of 
Poincare  (1892).  For  detailed  discussion  see  Ref.  [5]. 

k 2 

Theorem;  Let  be  a C (kM)  vector  field  on  R such  that 

x^(0)  = 0 for  all  y and  X = (Xp»0)  is  also  Let  dX^  (0,0)  denoting  the 

linear  part  of  X^  have  two  distinct,  complex  conjugate  eigenvalues  X(y)  and 
X(vi)  such  that  for  y>0.  Re  X(y)>0.  Also  let  | 


Then: 


(a)  There  is  a function  y;(-e,e)-^R  such  that  (xj^,0,y  (Xj^) ) 

is  on  a closed  orbit  of  period  ^21t/|X(0)  | emd  radius  growing  like 
^ of  X for  Xj^-^O  and  such  that  y(0)«=0. 

3 

(b)  There  is  a neighborhood  U of  (0,0,0)  in  R such  that  any 
closed  orbit  in  U is  one  of  those  aQ»ve.  Fvurthermore,  if  0 is 
a "vague  attractor"  for  X^,  then 

(c)  y(Xj^)>0  for  all  Xj^#0  and  the  orbits  are  attracting. 

The  Hopf  Theorem  in  R^ 

Theorem;  Let  be  a kM,  vector  field  on  r” , with  all  the 

assumptions  of  the  cibove  Theorem  holding  except  that  we  assume  that 
the  rest  of  the  spectrum  is  distinct  from  the  two  assumed  single 
eigenvalues  X(y),X(y).  Then  conclusion  (?i)  is  true.  Conclusion 
(b)  is  true  if  the  rest  of  the  spectrum  remains  in  the  left  half 
pleme  as  y crosses  zero.  Conclusion  (c)  is  true  if,  relative 
to  X (y)  ,X(y) , 0 is  a "vague  attractor"  and  if  when  coordinates  are 


SCIENTIFIC  SYSTEMS,  INC. 


-230- 


chosen  so  that: 


dXQ(O) 


0 |X(0)| 

-|A(0)|  0 

0 0 


(0) 
d3X^ (0) 


X(0)^a(d3X  (0)) 


[i.e.X  (0)  does  not  belong  to  the  spectnim  of  d^X  (0) ] . 


A. 4 Global  Tmpllcit  Function  Theorans  [26,27]; 

The  following  are  extensions  of  Palais  Theorem  [17] ; 


Theorem  1;  Let  f be  a continuously  differentiable  mapping 
_n  . . _m 


of  e”  X e"  into  E®.  There  then  exists  a unique  continuously 
differentiable  mapping  g of  E™  x e”  into  E®  with  the  property 
that  g(y»u)“x  for  all  x,yeE”  satisfying  f(x,u)“y,  provided 
that  i)  det[3f/3xj“0  at  each  point  in  E®  x e",  and  ii) 
for  each  uCe”,  | (f  (x^u)  1 1"**  as  ( |x|  (■+«». 


Theorem  2;  Let  f be  a continuously  differentiable  mapping  of 


e"  into  e”.  For  k=l,  . . n,  let  J,  denote  the  matrix 


consisting  of  the  elements  in  the  first  k rows  and  first  k 
columns  of  the  Jacobian  matrix.  Sty>pose  there  exist  a positive 
constant  e such  that 


|detJ. I |>e. 


I det'^2 


jdetJx 


f m • • f 


(det'^n 


detJ 


n-1 


>e. 


for  all  xCe”.  Then  there  exist  a unique  solution  of  xEe"  for 


each  yes  . 


1 


.i 


:] 


I 
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Appendix  B 


DerivationB  of  Equlllbriui  Equations  and  Analytic 
Stability  Criterion  for  Autorotation 

In  this  appendix,  further  details  of  the  derivations  of 
the  equilibrium  equations  in  terms  of  the  roll  rate  and  the  critical 
elevator  angle  for  autorotation  are  given  (cf.  Chapter  III) . 


B.l  Derivation  of  Equlibrium  Equation  in  Terms  of  Roll  Rate,  p. 

Starting  with  the  five  equations  of  motion  (gravity  ignored) : (a 
denotes  deviation  from  trim) . 

6 * + p(sina^  + a)  - r cos  + y ♦ y ■ 

fg(6,a,p,r^a,6r)  (B»l) 

a - q - p3+  z^a  + Zjg^e  - f^(B,a,p,q,6e)  (B.2) 

q » m^a  + m^q  + i2P>^  “ 

fq(B,a,p,q,r,5e)  (B.3) 

^ "a6a“^*  "r'  V "po^  ’ ^3*^ 

n^^Sr  » f^{B,a,p,q,r,6a,6r)  (B.4) 

i-V*  W*’  * V * V * V * V * ‘ra™ 

-ij^qr  + “ fp(B,oi,p,q,z;5a,6r)  (B.5) 

Setting  all  time  derivatives  to  zero  (i.e.,  fo,  f , f ,f  »f  “ 0), 

^ B a q r p 

five  equations  in  five  unicnowns  are  obtained.  Equation  B.2,  correspondingly 
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is  then  solved  for  q: 

q ■ p0  - s^o- *^5e  ■ fO,a,p,6e)  (B.6) 

Substituting  into  the  restaining  four  equations  gives 

0 « fg(6ja»P»r,6a,£e,5r)  (from  B.l)  (B.7) 

0 » f^(3fa,p,r,6a,5e,6r)  (from  B.3)  (B.8) 

2 

0 ■ f^(3,a,p,p  ,r,6a,6e,5r)  (from  B.4)  (B.9) 

0 ■ f^(3«a,p,r,6a,5e,5r)  (frost  B.5)  (B.IO) 

B.7  is  then  solved  for  r and  substituted  into  the  others,  leaving: 

0 » f^(3.a.p.p^,ja.3e.5r)  (frost  B.3)  (B.ll) 

0 ■ f^(3.a,p,p^,6a,6e,6r)  (frost  B.4)  (B.12) 

0»f  (3»  3 ^#0l,a^,p,p^,3a,6e,6r)  (from  b.5)  (B.13) 

P 

Next,  B.ll  is  solved  for  3 end  substituted  into  B.12  emd  B.13. 

0 - f^(a,p,p^,p^,p  ”^  ,6a,6e,6r)  (frost  B.4)  (B.14) 

0 - f 

q 


IB.15) 
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Flnally.  B.14  Is  solved  for  a: 


f ri (P  f *P^  »P”^ , 6a, 6e , 6r ) 
fr2<P»P^'P”^*^a*^®»*>^) 


This  Is  substituted  into  B.15,  leaving: 


0 - fp(p,p^,p^...p^,p  ^,6a,6e,6r) 


0 * f (p.p^. . .p^^,6a,6e,6r) 
P 


(B.16) 


(B.17) 


(B.18) 


Equation  B.18  is  a quartic  in  6a  [of.  Eq.  (B.3)].  It  sin^lifies  to 
a cubic  if  -i^^qr  tern  is  negelected  in  Eq.  (B.5)  and  to  a quadratic 
if  is  set  to  zero  in  Eq.  (B.4).  Finally,  if  all  the  nonlinear 

dynamic  terms  and  -i^qr  are  neglected,  a fifth  order  polynomial  in  p 
is  obtained  (cf.  Aircraft  B,  section  3.6.1). 


B.2  Analytic  Stability  Criterion  for  Autorotation. 

In  this  section,  we  derive  an  analytic  e:q>ression  for  critical  6e 
at  which  the  equilibrium  point  (p  >=  0,  r « 0,  3 “ 0,  6a  * 0,  6r  = 0) 
becomes  vinsteUsle  in  the  sense  that  a real  eigenvalue  crosses  the  origin 
from  the  left  half  plane.  First,  the  equilibrium  equations  corresponding 
to  q - 0 emd  a » 0 [Eqs.  (B.2)  and  (B.3)]  are  solved  for  d and  q in  terms 


of  5e. 


m z»  — m* 
q oe  oe 


m - m z_ 
a q a 


6e  K 6e 

a 


(B.19) 
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I q . ^ 6«  (B.20) 

m - ■ « ^ 

a q a 

N*  now  linaariz*  tha  aquations  of  aotlon  (B.1-B.5)  around 
(p.r-6-0,  a , q)  and  sat  tha  dataxainant  of  tha  jaocd>ian  to 
zaro.  This  laads  to  tha  following  aquation  for  critical  6a  (naglacting 
Ti^qr) : , 


2 

a5e  4-  b6a  4-  c > 0 

where 

• “ *^q*o**$q  *6a  * *rq  ^ra^  * *0o  * *ra  ^ra^ 

S ^ *rq  *r>  ^ *$  * **«  *r 


(B.21) 

(B.22) 

(B.23) 


c 


*6*6*^  ‘r*‘p 


(B.24) 


and 


0q  ng  + n^g 


+ ** 

_ go  g 

^ ng  + n^g 


-n 


B ng+n^g 


rq 


Vb, 

°B  ' V 


ra  n^  4-  n 


B * Vb 


^r-^B^B 
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IJ 

u 

LI 

IJ 

U 

U 

IJ 

D 


0 

D 

G 

0 

D 

0 


II 


For  th*  aircraft  B model.  Eq.  (B.21)  slnqplifles  further  since  a » o. 

# c 

Thus  the  critical  elevator  emgle.  ~ ^ deg.  See  Fig.  3.13. 

For  aircraft  H,  two  flight  conditions  from  Rhoades  and  Schuler  [1, 
Chapter  III]  %fere  examined  emd  for  both  these  conditions,  Eq.  (B.21) 
was  found  to  have  imaginary  roots.  Thus,  aircraft  H model  indicates 
that  autorotation  is  not  possible  for  6a  = 6r  0.  See  Fig.  3.48  for 
confirmation  of  this  observation. 

The  above  euialytical  results  are  confirmed  by  the  numerical 
simulation  results  given  in  section  3.6  of  this  report. 


i 
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Appendix  C 

Mcxaenclature  for  Chapter  3 

The  orientation  of  the  aircraft  body  axes  *’^9* 

with  respect  to  inertial  space,  or  eeurth-fixed  eoces,  is  defined  by  the 

following  Euler  angles,  in  sequence  from  reference  to  body  axes: 

i|>(yaw),  6 (pitch)  and  The  orientation  of  wind  €uces  (X^,Y^,Z^), 

stability  2uces  (X  ,Y  ,Z  ),  and  bo^  euces  and  the  conventions  for  linear 
s s s 

velocities  (u,v,w),  anguleu:  velocities  (p,q,r) , forces  (X,Y,Z),  moments 
(L,M,N) , and  control  surface  deflections  (6a, fie, fir)  are  all  shown  in 
Fig  .3.1.  The  symbols  and  definitions  used  in  this  report  2u:e  in  as  much 
agreement  with  convention  as  cleurity  allows.  The  body  axes  are  oriented 
so  as  to  align  with  the  stability  axes  for  straight  and  level  flight  in 
the  reference  flight  condition. 


B 

a 

“o 

fia 

fie 

fir 


u,v,w 

(I) 


X,Y,Z 


L,M,M 


orientation  (Euler)  angles  of  reference  to  body  axes. 

reference  value  of  6 

angle  of  sideslip 

angle-of-attack 

reference  value  of  a 

aileron  co—mnd  deflection 

elevator  command  deflection 

rudder  cosmund  deflection 

body  axis  components  of  linear  velocity 

body  axis  components  of  angular  velocity 

vector  of  angular  velocity  components,  (p,q,r) 

body  axis  components  of  force 

body  axis  components  of  moment 


1 


J 

Q 

u 

tJ 
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mass  (context  distinguishes  this  from  pitching 
moment  per  Iy#m) . 


I ,I  ,I  moments  of  inertia 
X y z 


matrix  (or  tensor)  of  moments  of  inertia 


il  - <1,  - 


^2  ” ^^2  " > Non-dimensional  Inertia  Coefficients 

^3  “ <^y  - V/"z  I 

E vehicle  kinetic  energy  jL 

2 2 2 ^ 

V total  aircraft  velocity  = (u  + v + w ) 

c wing  mean  aerodynamic  chord 

b wing  spem 

S lifting  surface  area 

g acceleration  of  gravity  («9.8,  metric,  32.2,  Eng.) 

p air  density 

— 12 

q dynamic  pressure,  ; equilibrium  pitch  rate  (distinguished  by 

context  from  dynamic  pressure) 

A-  rolling  moment  per 

m pitching  moment  per 

n yawing  moment  per  I 

z 

y side  force  (over  aircraft  mass  and  speed) 

z aerodynamic  force  along  axis  (over  mass  and  speed) 

Cy  - .y/(qs) 

C,  ,C  ,C  ■ L/(qSb),  N/('i5b) 

Ann 
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Superscrlpts 

equilibrium,  or  steady  state,  value 
• time  derivative  in  inertial,  or  earth-fixed,  coordinates 

0 reference  value 


Subscripts 

oi,S,a,P,q»r»  6a,6e,6r  denote  partial  derivatives  with  respect 

to  the  relevant  quantity  (viz,  C , C*,  C , C for  aircraft  A amd  B; 

y X n n 

emd  y,z,£,m,n  for  aircraft  H) . Thus,  for  ex2uiiple. 
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Special  Definitions 
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